Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,9280TA60.9% 95.5 / 168.3 139intergenic (‑1011/‑157)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base A (49/35);  total (79/60)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 143 total aligned reads displayed.

CTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTT  >  CP000730/1929781‑1930076
                                                                                                                                                       |                                                                                                                                                                
cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                                                                                                                                                                 >  1:1121636/1‑149 (MQ=255)
 tgggCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                                                                                                <  1:277153/146‑1 (MQ=255)
    gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                                                            >  1:183987/1‑149 (MQ=255)
    gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                                                            >  1:49257/1‑149 (MQ=255)
    gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCGCTCAAGCGTCA‑‑CTCGCCGAATAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                                                                                            >  1:439010/1‑149 (MQ=255)
     cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                                           >  1:754182/1‑149 (MQ=255)
     cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                                           >  1:132276/1‑149 (MQ=255)
      ttGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCTTTGCGTTTTTTTCTTAt                                                                                                                                                           >  1:981883/1‑149 (MQ=255)
      ttGTATATTCACAGATGGTATAATAACAGAGTCGCCTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                                          <  1:355960/149‑1 (MQ=255)
       tgtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                                         >  1:895637/13‑149 (MQ=255)
       tgtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                                         >  1:602427/13‑149 (MQ=255)
       tgtatattcacagatggtATAAAAATAGAGTCGTCTATCTCTCAGGCGTCAAAT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                                         >  1:616745/13‑149 (MQ=255)
       tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGGGTTTTTTTCTTAtt                                                                                                                                                          >  1:393240/1‑149 (MQ=255)
       tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt                                                                                                                                                          >  1:807893/1‑149 (MQ=255)
        ggagatgcacagatggtATAATAATAGAGGCGCCTATCTCTTAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                                        <  1:517121/138‑1 (MQ=255)
        gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                                        <  1:896506/149‑1 (MQ=255)
         tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                       >  1:839872/11‑149 (MQ=255)
         tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTTTTTTCGTTCACATCATAGCACCAGTCATCAGTGTATGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                                        <  1:729361/149‑1 (MQ=255)
         tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                                        <  1:237891/149‑1 (MQ=255)
         tataTTCACAGATGGTATAATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                       <  1:386217/149‑1 (MQ=255)
         gatattcacagatggtATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                                       <  1:52371/139‑1 (MQ=255)
          atattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                                      <  1:1114300/140‑1 (MQ=255)
          atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGTTGCGTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                                                       <  1:491541/149‑1 (MQ=255)
            attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                                                                                                    <  1:905206/142‑3 (MQ=255)
              tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                                  >  1:992980/6‑149 (MQ=255)
              tcacagatggtATAATAATAGAGTCGCCTATCTCACAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGACCCAGACATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttcgt                                                                                                                                                  >  1:487263/6‑145 (MQ=255)
                gCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                                                 >  1:322561/1‑149 (MQ=255)
                acagatggtATAATAATAGAGTCGTCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                                >  1:563005/4‑148 (MQ=255)
                acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                                >  1:76401/4‑148 (MQ=255)
                 cagatggtAAAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                               <  1:424423/147‑3 (MQ=255)
                 aagatagaataataAAAGAGCCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACCTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                               <  1:303495/141‑3 (MQ=255)
                  cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                                                               <  1:139008/149‑1 (MQ=255)
                  agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTCATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                              <  1:990328/148‑1 (MQ=255)
                  agatggtATAATAACAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                              <  1:222157/148‑1 (MQ=255)
                   gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGTCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                             >  1:864963/1‑149 (MQ=255)
                     tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a                                                                                                                                       <  1:153888/146‑1 (MQ=255)
                     tggtataATAACAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a                                                                                                                                       <  1:303724/146‑1 (MQ=255)
                      ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑a                                                                                                                                    <  1:198052/147‑1 (MQ=255)
                      ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑a                                                                                                                                    <  1:378791/147‑1 (MQ=255)
                      ggtataATAATAGAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑a                                                                                                                                    <  1:125056/147‑1 (MQ=255)
                       gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGTCTTAGGAAACGC‑‑‑‑a‑‑ac                                                                                                                                   >  1:653214/2‑148 (MQ=255)
                       gtataATAATAGAGTACCCTATCTCTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGGTAATCATTTTCGTTAACATCTGTGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑ac                                                                                                                                   <  1:541532/148‑2 (MQ=255)
                        tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑At                                                                                                                                   >  1:1030441/1‑149 (MQ=255)
                        tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGA‑‑‑‑A‑‑At                                                                                                                                   >  1:1030367/1‑149 (MQ=255)
                        tataATAATAGAGTCGCCTATCTCTCAGGGGTGGATTAGCGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑aca                                                                                                                                  >  1:825801/1‑147 (MQ=255)
                        tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑aca                                                                                                                                  <  1:350166/149‑3 (MQ=255)
                         atattaGTAGGGTGGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATa                                                                                                                                  <  1:143702/149‑1 (MQ=255)
                         ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaa                                                                                                                                 <  1:1035910/149‑4 (MQ=255)
                         ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATa                                                                                                                                  <  1:934251/149‑1 (MQ=255)
                            ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGATTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                              >  1:293289/5‑143 (MQ=255)
                            ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                              >  1:65428/5‑143 (MQ=255)
                            ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGGTTGCGTATTATCCTTATTGCCTTAGTAAACGC‑‑‑‑a‑‑acaagaa                                                                                                                              >  1:882551/5‑143 (MQ=255)
                            ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                              <  1:796739/149‑7 (MQ=255)
                            ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaata                                                                                                                              <  1:292042/149‑7 (MQ=255)
                             tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGAGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACaa‑‑a                                                                                                                            >  1:587387/1‑149 (MQ=255)
                             taatagCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataa                                                                                                                             <  1:856797/146‑8 (MQ=255)
                             taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTAATATACTTATTGGCTTAGTAAACGC‑‑‑‑accaaataa                                                                                                                               >  1:52079/4‑141 (MQ=255)
                               caattgTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatt                                                                                                                           >  1:1056579/6‑140 (MQ=255)
                                tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGCGCTGACTATGTGTGTTCACAAGATAGCACCAGTCATCAGTGGCTGTGCCGTTGCGTTTTTTTCTTATTGGCGAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATa                                                                                                                         >  1:1103212/1‑149 (MQ=255)
                                aagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataatta                                                                                                                          <  1:1050621/148‑11 (MQ=255)
                                 agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaa                                                                                                                         >  1:286889/1‑138 (MQ=255)
                                 agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataagtaa                                                                                                                         >  1:178157/1‑138 (MQ=255)
                                  gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTTGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAgg                                                                                                                       >  1:677342/1‑149 (MQ=255)
                                  gagTCGCCTATCGCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTGCTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGa‑‑‑‑a‑‑acaaataattaat                                                                                                                        >  1:801831/1‑137 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:167390/1‑136 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:122637/1‑136 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:112014/1‑136 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:243241/1‑136 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:895275/1‑136 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:701008/1‑136 (MQ=255)
                                   agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑a‑‑acaaataattaata                                                                                                                       >  1:853568/1‑136 (MQ=255)
                                    gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGA‑‑‑‑A‑‑ATACAA‑‑CATAGGTg                                                                                                                     >  1:589612/1‑149 (MQ=255)
                                      cGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGAc                                                                                                                   >  1:665176/1‑149 (MQ=255)
                                                 cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGGTAAGTATACGC‑‑‑‑A‑‑ATACAA‑‑AATCGGTAGCTTATAGCGGCa                                                                                                    >  1:244120/1‑149 (MQ=255)
                                                         aaTTTAGACGCAGAGGGGAGGTGTATAAGGTGATGCTAATTTCTGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATaaa                                                                                            <  1:866081/149‑1 (MQ=255)
                                                                  gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATGCAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc                                                                                   >  1:274713/1‑149 (MQ=255)
                                                                  gCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGGTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc                                                                                   >  1:65542/1‑149 (MQ=255)
                                                                      agagGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCATGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCTATAAAAATCCCCTCACTa                                                                               <  1:814817/149‑1 (MQ=255)
                                                                                       ggTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                                <  1:427533/149‑1 (MQ=255)
                                                                                         gTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAGTGgggg                                                               >  1:10834/1‑149 (MQ=255)
                                                                                           tGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                              <  1:354428/149‑1 (MQ=255)
                                                                                                cTTATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑C‑‑ATGCAT‑‑GGGAGGTGACATATAGCCGCACCCATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                                         <  1:214472/149‑1 (MQ=255)
                                                                                                         gTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGt                                                >  1:1074633/1‑149 (MQ=255)
                                                                                                          ttCACATCATAGCACCAGTCATCAGTGGCTGTTCCATTGCGTTTTTATCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                                               >  1:721197/1‑149 (MQ=255)
                                                                                                          ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                                               >  1:899438/1‑149 (MQ=255)
                                                                                                          ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                                               >  1:323030/1‑149 (MQ=255)
                                                                                                           tCACATCATAGCACCAGTCATCAGTGGCTGTGCCAATGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                              <  1:854300/149‑1 (MQ=255)
                                                                                                           ccacaTCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAA‑‑ATACTT‑‑AATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGgtgcatt                                                   <  1:84037/148‑5 (MQ=255)
                                                                                                            cacaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaaa                                             >  1:42761/1‑149 (MQ=255)
                                                                                                             acaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCCATAAAAATCCCCTCACTACCGCAAATAGTGTGGAGAATGGTGTATAAGTAAAt                                            >  1:432377/1‑149 (MQ=255)
                                                                                                                   tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGC‑‑‑‑AACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGAAAAAATCCCCACACTACTGCAATAGTGTGGGGATTGGTGcattaagtgtta                                           >  1:975997/1‑140 (MQ=255)
                                                                                                                   tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt                                      <  1:550275/149‑1 (MQ=255)
                                                                                                                        ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCg                                 >  1:936024/1‑149 (MQ=255)
                                                                                                                                  gTGGCTGTGCCATTGGTTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                       <  1:87447/149‑1 (MQ=255)
                                                                                                                                  gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                       <  1:648345/149‑1 (MQ=255)
                                                                                                                                  ctGGCTGTCCCATTGCGTTATTTTCTTAGTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                       <  1:481020/148‑1 (MQ=255)
                                                                                                                                          gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGc               >  1:373903/1‑149 (MQ=255)
                                                                                                                                             cgtGCGTTTTTTTCTTAGTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa            <  1:478425/147‑1 (MQ=255)
                                                                                                                                             aTTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa            <  1:741225/149‑1 (MQ=255)
                                                                                                                                                       tttCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAAGCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtttt  >  1:509531/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                                
CTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGC‑‑‑‑A‑‑ATACAA‑‑AATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTT  >  CP000730/1929781‑1930076

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: