Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,9250TA68.7% 855.3 / inf 659intergenic (‑1008/‑160)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (82/122);  new base A (282/165);  total (369/290)
Fisher's exact test for biased strand distribution p-value = 6.93e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 698 total aligned reads displayed.

GTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCT  >  CP000730/1929777‑1930073
                                                                                                                                                         |                                                                                                                                                                
gcgACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCATGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                                                                                                  <  2:119259/147‑1 (MQ=255)
    atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatc                                                                                                                                                             >  1:734310/4‑146 (MQ=255)
          ctGTATATTCACAGATTGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                                       <  1:1492841/148‑1 (MQ=255)
            gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                                     <  1:1229962/138‑1 (MQ=255)
             tattgtggCCGTTGATATAATAGTAGAGTCCCCTATCGCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                                     <  1:144010/142‑1 (MQ=255)
             tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCCATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    >  2:1313381/11‑149 (MQ=255)
             tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  1:75659/139‑1 (MQ=255)
             tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  1:288505/139‑1 (MQ=255)
             tatatccacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  1:927666/139‑1 (MQ=255)
             tatagtcacagatggtATAATAATAGGGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  2:734026/139‑1 (MQ=255)
             tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  1:1637308/149‑1 (MQ=255)
             tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  1:1874051/149‑1 (MQ=255)
             tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  2:1143094/149‑1 (MQ=255)
             tataTTCACAGACGATATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                    <  2:52544/149‑1 (MQ=255)
              atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                                                    <  1:84006/149‑1 (MQ=255)
               tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                                  <  1:1838573/141‑2 (MQ=255)
               tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                                  <  1:1264212/141‑2 (MQ=255)
                attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTGCTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                                                 >  1:146746/8‑147 (MQ=255)
                attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                                                 >  2:699971/8‑147 (MQ=255)
                atacacagagagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                                                 >  2:1735552/13‑147 (MQ=255)
                 ttcactgatggtATAATAATAATATCCCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                <  2:409447/143‑1 (MQ=255)
                 ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                                                >  1:307216/7‑146 (MQ=255)
                 ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                >  1:641925/7‑149 (MQ=255)
                 ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                >  2:48930/7‑149 (MQ=255)
                   cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGAAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGAGGCTGGGCTGTTGCGTATACTACTTATTGGCTTAGTa                                                                                                                                              >  2:602524/5‑149 (MQ=255)
                   cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                              >  2:1905952/5‑149 (MQ=255)
                    acagatggtATAATAATAGAGTCGCCTATCTCTCCGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                             >  2:429292/4‑148 (MQ=255)
                    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                             >  2:1261095/4‑148 (MQ=255)
                    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                             >  2:1467004/4‑148 (MQ=255)
                    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTAATATACTTATTTGCTAAGTaa                                                                                                                                             >  2:1374125/4‑148 (MQ=255)
                     cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                            >  2:1553932/3‑147 (MQ=255)
                      agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                           >  1:1381242/2‑149 (MQ=255)
                      agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                           >  1:1654069/2‑149 (MQ=255)
                        atggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCCTTACCTGGCTGATCATTTTTGTTCACATCATAGCACCAGTCATCAATGGCTGTGCTGTTGCTTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                         >  1:1974918/1‑149 (MQ=255)
                         tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a                                                                                                                                    <  1:793301/146‑1 (MQ=255)
                          ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑a                                                                                                                                 <  2:1594984/147‑1 (MQ=255)
                           gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑ac                                                                                                                                >  1:1090730/2‑148 (MQ=255)
                            tataATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑aca                                                                                                                               <  1:688313/149‑3 (MQ=255)
                             ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaa                                                                                                                              <  2:971706/149‑4 (MQ=255)
                             ataataataGAGTCGCCTATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaa                                                                                                                              >  1:1024161/1‑146 (MQ=255)
                             ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATa                                                                                                                               <  2:92787/149‑1 (MQ=255)
                              taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaa                                                                                                                             <  2:233650/149‑5 (MQ=255)
                              taataaGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAATACTGGGATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaa                                                                                                                             <  1:950476/149‑5 (MQ=255)
                              aaataatAGAATCGCATATCTTCTATGCGTCAA‑‑TCGACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCATTCATCAGTTGCTGTGCTGTTGATTATTATACTTATTGGATTAGTTTACGC‑‑‑‑a‑‑acaaa                                                                                                                             <  2:1771345/148‑5 (MQ=255)
                               aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaat                                                                                                                            >  2:324231/1‑144 (MQ=255)
                                ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                           >  2:856857/5‑143 (MQ=255)
                                ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                           <  1:1768451/149‑7 (MQ=255)
                                ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                           >  1:417695/1‑143 (MQ=255)
                                 tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACaa‑‑a                                                                                                                         >  1:1009797/1‑149 (MQ=255)
                                 tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACaa‑‑a                                                                                                                         <  1:1729725/149‑1 (MQ=255)
                                 taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataa                                                                                                                          <  2:891231/146‑8 (MQ=255)
                                 taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataa                                                                                                                          <  1:494378/146‑8 (MQ=255)
                                   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatt                                                                                                                        >  1:1421663/2‑140 (MQ=255)
                                   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTAGACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatt                                                                                                                        >  1:1340367/2‑140 (MQ=255)
                                   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatt                                                                                                                        >  1:539425/2‑140 (MQ=255)
                                   aaagaaTCGCATATCCGTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatt                                                                                                                        >  2:1227964/3‑140 (MQ=255)
                                    tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatta                                                                                                                       >  1:1926062/1‑139 (MQ=255)
                                    tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatta                                                                                                                       >  1:1048257/1‑139 (MQ=255)
                                    tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTACCGGCTTCGTAAACGC‑‑‑‑a‑‑acaaataatta                                                                                                                       >  2:1353694/1‑139 (MQ=255)
                                    tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatta                                                                                                                       >  2:1536006/1‑139 (MQ=255)
                                    tagAGTCGCCTATCTCTCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatta                                                                                                                       >  1:1198743/1‑139 (MQ=255)
                                     agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAg                                                                                                                     <  1:1790225/149‑1 (MQ=255)
                                      gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaat                                                                                                                     >  2:1861085/1‑137 (MQ=255)
                                      gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaat                                                                                                                     >  2:811380/1‑137 (MQ=255)
                                       agTCGCCTATCTCTCAGGCGTCAATT‑‑GATGAAGAGAGGAGGTGCATTACTTTCCTAATCATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACGTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                    >  2:1097273/1‑136 (MQ=255)
                                       agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑agcaagtaactaaaa                                                                                                                      >  2:67518/1‑135 (MQ=255)
                                       agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                    >  1:1902671/1‑136 (MQ=255)
                                       agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                    >  2:2007947/1‑136 (MQ=255)
                                       agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑aaaaataattaata                                                                                                                    >  2:141328/1‑136 (MQ=255)
                                        gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTg                                                                                                                  <  1:896826/149‑1 (MQ=255)
                                        gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaatag                                                                                                                   >  2:375065/1‑135 (MQ=255)
                                        gTCGCCTATCTCTCAGGCGTCAATGTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTg                                                                                                                  <  1:258630/149‑1 (MQ=255)
                                                  tctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGcc                                                                                                     <  1:1315418/149‑1 (MQ=255)
                                                     cagtcttCAATGTAGACGCAGAGAGGAGGTGTATAAGGTGAGTCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCa                                                                                                 <  2:1003319/143‑1 (MQ=255)
                                                        gCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCGAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGGGACATATAGCCGCACCa                                                                                              >  1:1517381/1‑149 (MQ=255)
                                                         cGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCaa                                                                                             <  2:26369/149‑1 (MQ=255)
                                                         cGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCaa                                                                                             <  2:1697115/149‑1 (MQ=255)
                                                                    aCGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATcccc                                                                                  <  1:170445/149‑1 (MQ=255)
                                                                    aCGCAGAGAGGAGGGGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTCGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATcccc                                                                                  <  2:1828959/149‑1 (MQ=255)
                                                                      gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc                                                                                <  1:216808/149‑1 (MQ=255)
                                                                             ggaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg                                                                         <  1:1165489/149‑1 (MQ=255)
                                                                             ggaggTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg                                                                         <  1:1590388/149‑1 (MQ=255)
                                                                                     ataAGCTGTTGCTTCTTTTCTTTCACTGCCTAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGt                                                                 <  2:105986/149‑1 (MQ=255)
                                                                                           ggTGATTCTTAATTTCGTTCACATCATACCACCAGTCATCAGCGGCTGTGCCATTGCGTTTCTTTCTTATTGGCTAAGCACACCC‑‑‑‑A‑‑ATACAA‑‑AAGAGGTGACATATAGCCGCACCACTAAAACTCCCCTCACTACCGCAAATAGTGAgg                                                             <  1:37878/149‑1 (MQ=255)
                                                                                                      ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                                     <  1:883357/149‑1 (MQ=255)
                                                                                                      ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                                     >  1:340288/1‑149 (MQ=255)
                                                                                                         ttttCGTTCACATCATAGCACCGGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGta                                                  <  1:160513/149‑1 (MQ=255)
                                                                                                         tttgCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTAACATATAGCCGCTCCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGta                                                  <  2:1504103/149‑1 (MQ=255)
                                                                                                             cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                              <  1:1550079/149‑1 (MQ=255)
                                                                                                              gTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑AACAAATAATTAATAGGCGACTTTATAAGCCAGCAACGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGgtgca                                                  <  2:117251/149‑3 (MQ=255)
                                                                                                                  acaTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑AACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGcattaa                                              >  1:597261/1‑146 (MQ=255)
                                                                                                                      catAGGACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                                     <  1:1135057/149‑1 (MQ=255)
                                                                                                                          gCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAA‑‑ATAATT‑‑AATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGgtgcattgattgttaat                                      >  2:148588/1‑135 (MQ=255)
                                                                                                                          gCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑AACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGcattaattgttaaa                                      >  1:570530/1‑138 (MQ=255)
                                                                                                                                gTCATTAGTGGCTGTGCCATTGCGTGTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg                           <  2:778076/149‑1 (MQ=255)
                                                                                                                                  gatcaGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTc                         <  2:63612/148‑1 (MQ=255)
                                                                                                                                   ttcaGTGACTGGCTCATTGCGTTGTTTTCTCATTGGCTAAGTGGAACC‑‑‑‑A‑‑ATACAA‑‑AACAGGTGACATATAGCTGCGATAATAAAAGGCCCCTCACTACCGCCAATGGTGGGGGGGTTGGTGTATAAGTAAATACTTATTTTCGTTGTCt                        <  2:1650725/148‑1 (MQ=255)
                                                                                                                                              tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATACAGCCCCGCCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg             <  1:464095/149‑1 (MQ=255)
                                                                                                                                                         ttttttCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCt  >  1:66448/1‑149 (MQ=255)
                                                                                                                                                         ttttttCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCt  >  2:1963189/1‑149 (MQ=255)
                                                                                                                                                         |                                                                                                                                                                
GTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCT  >  CP000730/1929777‑1930073

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