Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8730CA60.3% 24.1 / 82.8 64intergenic (‑956/‑212)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/9);  new base A (21/17);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAT  >  CP000730/1929734‑1930018
                                                                                                                                             |                                                                                                                                                 
gATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATGGAGACGCAGAGAGGAGGGGTA‑TAA‑GGTGAAGCTTAGTTTc                                                                                                                                           >  1:275481/1‑148 (MQ=255)
        cTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGACATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCAcat                                                                                                                                   >  2:53708/1‑148 (MQ=255)
           tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGTTTATTTTCGTTCAcatcat                                                                                                                                >  1:101695/1‑148 (MQ=255)
           tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCAcatcat                                                                                                                                >  2:201278/1‑148 (MQ=255)
                  gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                                         >  1:38192/1‑148 (MQ=255)
                  gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                                         <  1:579/148‑1 (MQ=255)
                      cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGtca                                                                                                                     <  1:316117/148‑1 (MQ=255)
                              ttATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                             >  1:177856/1‑148 (MQ=255)
                              ttATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                             >  1:302640/1‑148 (MQ=255)
                              ttATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                             >  2:328741/1‑148 (MQ=255)
                              ttATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                             >  2:239465/1‑148 (MQ=255)
                              ttATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                             >  2:173584/1‑148 (MQ=255)
                                aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                           <  1:169396/148‑1 (MQ=255)
                                     aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCa                                                                                                      <  2:251631/148‑1 (MQ=255)
                                       gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                                                    >  2:101087/1‑148 (MQ=255)
                                               atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                                           <  1:6538/145‑6 (MQ=255)
                                                tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                          <  2:129645/146‑7 (MQ=255)
                                                 atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                         >  1:308741/2‑144 (MQ=255)
                                                       gtatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                   >  1:59196/12‑148 (MQ=255)
                                                        tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                  <  2:371230/138‑1 (MQ=255)
                                                            ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                              <  2:279963/142‑3 (MQ=255)
                                                            ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                              >  2:11682/7‑146 (MQ=255)
                                                            ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtta                                                                              >  2:257843/7‑146 (MQ=255)
                                                             tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                             >  1:227408/6‑145 (MQ=255)
                                                              tGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGGTAAg                                                                             >  2:117982/1‑148 (MQ=255)
                                                              cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                            <  1:362772/144‑1 (MQ=255)
                                                                cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                          <  1:305347/146‑2 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTCCGTTCACATCATAGCAACAGTCCTCAGTGGCAGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         >  2:382852/2‑146 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         <  2:312688/147‑3 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         >  1:137948/2‑146 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         <  1:280287/147‑3 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         >  1:372082/2‑146 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         <  1:11970/147‑3 (MQ=255)
                                                                 agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                         >  2:208644/2‑146 (MQ=255)
                                                                    tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                      <  2:222396/145‑1 (MQ=255)
                                                                     ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                     >  1:218581/3‑148 (MQ=255)
                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                  <  2:47108/148‑3 (MQ=255)
                                                                        ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                  <  2:327037/148‑3 (MQ=255)
                                                                           ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                               >  1:381175/5‑143 (MQ=255)
                                                                           ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                               >  2:264232/5‑143 (MQ=255)
                                                                           ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                               >  1:142730/5‑143 (MQ=255)
                                                                           ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                               >  1:368747/5‑143 (MQ=255)
                                                                           ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                               >  2:126641/5‑143 (MQ=255)
                                                                           ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                               >  2:219711/1‑143 (MQ=255)
                                                                            tagtagAGTCGCCTATCTCTCAGGCGACAATTTCGACGCAGAGAGGAGGCGTA‑TAA‑GGCGATGATTATGTTCGTGCACATCATGGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGGCTAAGTAGACGCAATACaa                                                               <  2:328291/148‑1 (MQ=255)
                                                                            taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                              >  2:321650/4‑142 (MQ=255)
                                                                             aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                             >  2:379018/3‑141 (MQ=255)
                                                                             aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                             >  2:177323/3‑141 (MQ=255)
                                                                              gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                                             <  1:212971/148‑1 (MQ=255)
                                                                               tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                           <  1:92273/148‑10 (MQ=255)
                                                                               tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                           <  2:331011/148‑10 (MQ=255)
                                                                                  agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                        <  2:183425/148‑13 (MQ=255)
                                                                                  agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                        >  2:370969/1‑136 (MQ=255)
                                                                                   gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                        >  1:102863/1‑148 (MQ=255)
                                                                                   gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                       >  1:379556/1‑135 (MQ=255)
                                                                                   gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTACTAAACGCaacaaataattaata                                                       <  1:380798/148‑14 (MQ=255)
                                                                                                cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑AAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                           <  1:65201/148‑1 (MQ=255)
                                                                                                 aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                                          <  2:342756/148‑1 (MQ=255)
                                                                                                 aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTCATGCTTATTTTCGTACACATCATAGCACCAGTCACCAGTAGCTGTGCCATTGCGTTTTTTTCTTACTGGCTAAGTAGGCGCCATACATACTAGGTGACATATAGCCGCa                                          <  2:255193/148‑1 (MQ=255)
                                                                                                           ttAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                <  2:221979/148‑1 (MQ=255)
                                                                                                                 gCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                          >  1:344543/1‑148 (MQ=255)
                                                                                                                 gCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTAGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                          >  1:219464/1‑148 (MQ=255)
                                                                                                                             tgtA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTTCCATGGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAAAGGTGACATATAGCCGCACCAATAAAAATCTCCTCACTACCGTaaa              >  1:164400/1‑148 (MQ=255)
                                                                                                                                           tGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAt  >  2:127199/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAT  >  CP000730/1929734‑1930018

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