Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,970 T→C 59.8% E194E (GAA→GAG)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,9700TC59.8% 526.9 / inf 828E194E (GAA→GAG) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (176/155);  new base C (234/258);  total (411/417)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 831 total aligned reads displayed.

AACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGC  >  CP000730/35848‑36110
                                                                                                                          |                                                                                                                                             
aatcACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTa                                                                                                                            <  2:699055/138‑1 (MQ=0)
                   cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATAc                                                                                                         <  2:433822/141‑1 (MQ=255)
                   cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATAc                                                                                                         <  1:612302/141‑1 (MQ=2)
                   cAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCAGGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATAc                                                                                                         >  1:828045/1‑141 (MQ=11)
                    aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTGTTTTTATATAGAGCGTAAATACATTCAATACTTT‑TTAAAGTAGTGTTTGCCGTATCGATAcg                                                                                                        >  2:205711/1‑140 (MQ=11)
                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACtt                                                                                                       >  1:945615/1‑141 (MQ=255)
                     gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                       >  2:534153/1‑141 (MQ=12)
                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                       <  2:598443/141‑1 (MQ=17)
                     gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                       >  1:577214/1‑141 (MQ=17)
                      tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                      >  1:129961/1‑141 (MQ=32)
                        ttcttttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                                                    >  1:764059/6‑141 (MQ=12)
                        ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                                                    <  1:725291/141‑1 (MQ=17)
                         tCGCTTAACTGTCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                                                                                   <  1:463973/141‑1 (MQ=12)
                          cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                                                  >  2:93772/4‑141 (MQ=17)
                            tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTGACGGCTGTATTGATACTTTGATAcc                                                                                                <  2:752790/140‑1 (MQ=11)
                              tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATAAACTCAATACCTT‑TTAAAGTACTCTTTGCTGTATTGATACTTTGATACCtt                                                                                              >  1:932078/1‑141 (MQ=12)
                               aaCTTGCTAGCATGATGCTCATTTAGTAGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGatactttgatactttg                                                                                             >  2:228124/1‑141 (MQ=11)
                               aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTg                                                                                             >  2:471446/1‑141 (MQ=255)
                               aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                                                                             <  2:865667/141‑1 (MQ=17)
                                aCTTGCTAGCATGATGCTAATTTCGTGGGATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACTCTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                                                                            <  2:279873/141‑1 (MQ=12)
                                 cTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtc                                                                                           >  2:403710/1‑141 (MQ=255)
                                    gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                                                        <  1:441723/141‑1 (MQ=17)
                                     ttAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                       <  1:973948/140‑1 (MQ=17)
                                     cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                       >  1:853514/1‑141 (MQ=32)
                                      tAGAATGATGCTAATTTCGTAGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTATATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                      <  1:49640/141‑1 (MQ=255)
                                       aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                                                                                     <  1:475996/141‑1 (MQ=17)
                                        gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTa                                                                                    <  1:906494/141‑1 (MQ=17)
                                        gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATAGATACTTTGATACCTTGTCTTTCTTa                                                                                    >  1:1007847/1‑141 (MQ=17)
                                          atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                                                                                  >  2:341772/1‑141 (MQ=18)
                                          atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                                                                                  >  1:301526/1‑141 (MQ=21)
                                            gatgCTATTTTCGTGGCATGGCGAAAATTCGTAGATCTGAACTGAGCCGCGGTTTTTTTTATCGACAGCGTAAATACACTCAAGAGCTG‑TTAAAGTATTCTTTGCTGTGCTGATTCTTTGAGACCTTGTCTTTCATACttt                                                                                <  2:266550/141‑1 (MQ=2)
                                                cTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAata                                                                            >  1:955381/1‑141 (MQ=25)
                                                 tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTTTTGATACTTTGATTTCTTGTCTTTCTTACTTTAatat                                                                           >  1:265813/1‑141 (MQ=11)
                                                 tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                                                                           >  1:887549/1‑141 (MQ=21)
                                                 tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGGAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatct                                                                           >  2:316379/1‑139 (MQ=255)
                                                  aaTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAatata                                                                          >  2:143429/1‑140 (MQ=255)
                                                   aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTGTTCTTTGCTGTATTTATGCTTTTATACCTTGTCTTTTTTACTTGAATATGa                                                                         >  1:73965/1‑141 (MQ=11)
                                                    tttCGTGGCATGGCGACAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATGGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATgga                                                                        >  2:461125/1‑139 (MQ=255)
                                                    tttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGAc                                                                        >  1:563066/1‑141 (MQ=25)
                                                     ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatatcccg                                                                       >  2:522173/1‑137 (MQ=255)
                                                     ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTATTTTTTACTTTAatatttag                                                                       >  2:1048918/1‑137 (MQ=255)
                                                      tcttGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                                                      <  1:522844/138‑1 (MQ=25)
                                                       cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt                                                                     <  2:980696/141‑1 (MQ=255)
                                                       cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAACATGACGGt                                                                     >  2:299882/1‑141 (MQ=255)
                                                        gcggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg                                                                    <  2:1038028/139‑1 (MQ=255)
                                                          ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt                                                                  >  2:18479/1‑141 (MQ=37)
                                                           gCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATc                                                                 <  2:484975/141‑1 (MQ=38)
                                                            cacagccaaaaCCCGCAGATCCGAAGAGACCTGCGGTTCGTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt                                                                <  1:12571/134‑1 (MQ=14)
                                                              tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                              >  2:387484/1‑141 (MQ=37)
                                                              tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                              >  1:640559/1‑141 (MQ=38)
                                                               ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                                                             <  1:1018132/141‑1 (MQ=255)
                                                               ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTATAATATGACGGTTATCTTGc                                                             <  1:668424/141‑1 (MQ=255)
                                                                 cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAAGACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTGTTTCTTACTTTAATATGACGGTGATCTTGCta                                                           >  1:710373/1‑140 (MQ=255)
                                                                 cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGCAAATACACTCAACACCTT‑TTAAAGTATTCTTTGCTGTATTGTTACTTTGGGACCTTGTCTTTCTTACTTTACTATGCCGGTGATCTTGCTc                                                           >  2:963442/1‑141 (MQ=11)
                                                                  gAAAATCCGTAGATCTGTAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCa                                                          <  2:678134/141‑1 (MQ=255)
                                                                    aaaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTACCGTATCCATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAAt                                                        <  1:201349/141‑1 (MQ=11)
                                                                     aaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                       <  1:807486/141‑1 (MQ=255)
                                                                     aaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                       >  2:1012570/1‑141 (MQ=255)
                                                                      aTCCGTAGATCTGACGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGGATTGATACTTTGATGTCTTGTATTTCTTACTTTAAGATGACGGTgggcgttcacgatga                                                      >  2:719978/1‑127 (MQ=11)
                                                                      aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGa                                                      <  1:749865/141‑1 (MQ=255)
                                                                      aTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGa                                                      <  1:539822/141‑1 (MQ=255)
                                                                        ccGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAgg                                                    <  1:626824/141‑1 (MQ=255)
                                                                           tAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCCATGAGGttt                                                 >  2:360361/1‑140 (MQ=255)
                                                                           tAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACACAATACCTT‑TTAAAGTCTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGGCGGTGATCTTGCTCAATGAGGTTa                                                 >  1:998322/1‑141 (MQ=255)
                                                                            aggggtGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGATCAATGAGGTTAt                                                >  2:651143/6‑141 (MQ=11)
                                                                             gATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                                               <  2:839886/141‑1 (MQ=255)
                                                                                cTGAAGCGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                            <  1:405487/141‑1 (MQ=255)
                                                                                cTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                            >  1:864279/1‑141 (MQ=38)
                                                                                cTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                            <  2:880915/141‑1 (MQ=255)
                                                                                    agagACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGGTGTATTGATTCTTTGATTTGTTGTCTTTTTTACTTTAATTTGGCGGTGGTCTTGCTCAATGAGGTTATTCTGatat                                        >  2:115234/1‑141 (MQ=11)
                                                                                    agagACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTACAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGatat                                        <  2:568003/141‑1 (MQ=255)
                                                                                      agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGATGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAttt                                      >  2:452076/1‑141 (MQ=255)
                                                                                       gtccTGCGGTTCTTTTTATATAGAGCCTAAATACACTAAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTc                                     <  2:806716/139‑1 (MQ=255)
                                                                                        accgacGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCg                                    <  2:247154/136‑1 (MQ=255)
                                                                                        aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCATTAAATTATTTAACT‑TTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCg                                    <  1:1037636/141‑1 (MQ=14)
                                                                                         ccTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGa                                   <  1:216666/141‑1 (MQ=255)
                                                                                           tGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGTTCTTGTTCAATGAGGTTATTCAGATATTTCGAtt                                 >  1:324408/1‑140 (MQ=255)
                                                                                           tGCGCTTCTTTTTATATAGAGCGGGTTGACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATg                                 <  2:740811/141‑1 (MQ=11)
                                                                                            acGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGt                                >  1:277519/2‑141 (MQ=255)
                                                                                               gTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACa                             <  1:425298/141‑1 (MQ=38)
                                                                                                 tCTTTTTATATAGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                           >  2:911925/1‑141 (MQ=37)
                                                                                                   tttttATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGa                         >  1:360333/1‑141 (MQ=40)
                                                                                                          ataGAGCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAg                  <  1:178545/141‑1 (MQ=255)
                                                                                                             gagCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtt               <  1:587217/141‑1 (MQ=255)
                                                                                                             gagCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtt               >  1:790833/1‑141 (MQ=255)
                                                                                                              agCGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGCTATTTCGATGTACAATGGCAGTCAGgttt              >  1:920417/1‑141 (MQ=255)
                                                                                                               gCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttta             >  2:633626/1‑141 (MQ=255)
                                                                                                                cGTAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATAAAGAAGTACAATGGCAGTCAGgtttaa            <  1:1053084/141‑1 (MQ=14)
                                                                                                                  tAAATACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagt          <  1:382807/141‑1 (MQ=255)
                                                                                                                   aaaTACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATACAGATATTTCGATGTACAATGACAGTCAGgtttaagtt         <  1:681599/141‑1 (MQ=255)
                                                                                                                   aaaTACATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt         <  1:236474/141‑1 (MQ=255)
                                                                                                                   aaaTACATGCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt         <  1:16000/141‑1 (MQ=255)
                                                                                                                   aaaTACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtt         >  2:778331/1‑141 (MQ=255)
                                                                                                                    aaTACATTCAATACCTT‑TTACAGTATTCTTTGCCGTATTGATACATTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTAGCTCAATGAGGTTATTCAGATCTTTCGATGTACACTGACAGTCAGgtttaagttt        >  2:582142/1‑141 (MQ=14)
                                                                                                                    aaTACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCCGTCAGgtttaagttt        >  1:69187/1‑141 (MQ=255)
                                                                                                                    aaTACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt        <  2:618368/141‑1 (MQ=255)
                                                                                                                     aTACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta       >  2:162303/1‑141 (MQ=255)
                                                                                                                     aTACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta       >  1:352336/1‑141 (MQ=255)
                                                                                                                      tACACTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTGCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa      >  2:667602/1‑141 (MQ=255)
                                                                                                                        cATTCAATACCTT‑TTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTaaaa    <  2:552189/141‑1 (MQ=255)
                                                                                                                          |                                                                                                                                             
AACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTT‑TTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGC  >  CP000730/35848‑36110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑