Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8610TC59.5% 48.2 / 69.7 74intergenic (‑944/‑224)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/16);  new base C (27/17);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCG  >  CP000730/1929728‑1929976
                                                                                                                                      |                                                                                                                      
gtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑Gt                                                                                                              >  2:329965/1‑141 (MQ=255)
        ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTa                                                                                                      <  2:387984/141‑1 (MQ=255)
            aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTAtttt                                                                                                  >  1:327432/1‑141 (MQ=255)
                    gattgAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTcaca                                                                                          >  1:248743/5‑141 (MQ=255)
                        gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCAcatcat                                                                                      >  2:68206/1‑141 (MQ=255)
                            cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCa                                                                                  >  2:127105/1‑141 (MQ=255)
                                     tATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatc                                                                         <  1:204100/141‑1 (MQ=255)
                                           tcggaagtagatatgCTAGTATATTCACAGATGGTATAATAATAGAGTCGCCTACCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGAAGGTGCA‑TTACT‑TTGCTAATCATTTACGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                  >  2:330566/14‑141 (MQ=255)
                                           aaGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                   <  1:16825/141‑1 (MQ=255)
                                              gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                               <  1:44075/131‑1 (MQ=255)
                                               aaGTGACTGTGCTTGTATATTTGCCAATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                               >  1:290106/1‑141 (MQ=255)
                                                 gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                            >  2:168459/8‑139 (MQ=255)
                                                 gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                            >  1:71290/8‑139 (MQ=255)
                                                  tagatatgCTTGTATATTCACAGATTGTATAATAATAGAGTCTCCTATCTCTCAGGCGGC‑AAGT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                           <  1:350334/135‑4 (MQ=255)
                                                  tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                           <  1:207285/135‑4 (MQ=255)
                                                  tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                           <  1:207455/135‑4 (MQ=255)
                                                   gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                           <  2:183901/141‑1 (MQ=255)
                                                   gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                           <  2:93420/141‑1 (MQ=255)
                                                   agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                          >  2:62631/6‑137 (MQ=255)
                                                   agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGGC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                          >  2:392242/6‑137 (MQ=255)
                                                    gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                         <  1:125257/137‑6 (MQ=255)
                                                    gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                         >  1:157456/5‑136 (MQ=255)
                                                    aCTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                          <  1:10158/141‑1 (MQ=255)
                                                      tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                       >  1:410236/3‑141 (MQ=255)
                                                      tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                       >  2:33219/3‑141 (MQ=255)
                                                       atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                      >  2:124294/2‑141 (MQ=255)
                                                       atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                      <  2:285202/140‑1 (MQ=255)
                                                       atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCACAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                      >  1:363226/2‑141 (MQ=255)
                                                           ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                  <  1:318392/141‑5 (MQ=255)
                                                             gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                                                 <  1:3211/141‑1 (MQ=255)
                                                               atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc                                               <  1:417919/141‑1 (MQ=255)
                                                                 attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactta                                            >  2:181286/8‑134 (MQ=255)
                                                                 atTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                             <  1:124244/141‑1 (MQ=255)
                                                                    cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                         >  2:78306/5‑141 (MQ=255)
                                                                     acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                 >  1:130734/4‑116 (MQ=255)
                                                                     acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                 <  2:130734/113‑1 (MQ=255)
                                                                      cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                       <  1:300150/139‑1 (MQ=255)
                                                                      cagatggtATAATAATAGAGTCGCCTATCTCTCAGGAGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                       >  1:59774/3‑141 (MQ=255)
                                                                         atggtATAATAAGAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTCATTGGCtta                                    >  2:85646/1‑139 (MQ=255)
                                                                         atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTCTTTCTCTTATTGTCTa                                     >  2:389478/1‑141 (MQ=255)
                                                                          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                 >  2:405295/4‑122 (MQ=255)
                                                                          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                 <  1:405295/124‑6 (MQ=255)
                                                                          tggtataATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                   <  2:198921/138‑1 (MQ=255)
                                                                           ggtataATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                  <  1:120233/139‑1 (MQ=255)
                                                                           gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                   >  2:59458/1‑141 (MQ=255)
                                                                           gatATAATAGTAGAGTCGCCTAGCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGGGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGCTTTTTTCTTATTGGCTGAg                                   >  2:217532/1‑141 (MQ=255)
                                                                            gtataATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                 >  2:31455/2‑141 (MQ=255)
                                                                                 ataataGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                            <  2:365298/141‑1 (MQ=255)
                                                                                 ataGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                             >  2:385191/1‑141 (MQ=255)
                                                                                   aatagAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                          <  1:70452/139‑1 (MQ=255)
                                                                                    gtagAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                          <  2:273971/141‑1 (MQ=255)
                                                                                      agagTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                       >  1:289441/1‑138 (MQ=255)
                                                                                       gagTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                       <  1:329965/141‑1 (MQ=255)
                                                                                       gagTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                   >  1:191104/1‑109 (MQ=255)
                                                                                       gagTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                   <  2:191104/112‑4 (MQ=255)
                                                                                       gagTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                      <  2:131607/141‑5 (MQ=255)
                                                                                       gagTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                      >  2:63469/1‑137 (MQ=255)
                                                                                        agTCGCCTATCTCTCAGGCGTCAAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                       <  2:345464/141‑1 (MQ=255)
                                                                                           cGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                   >  2:54759/1‑92 (MQ=255)
                                                                                           cGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                   <  1:54759/92‑1 (MQ=255)
                                                                                           cGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                  >  2:171395/1‑133 (MQ=255)
                                                                                            gCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                 >  1:173528/1‑132 (MQ=255)
                                                                                             ccTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTGGTCGGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                >  1:227300/1‑131 (MQ=255)
                                                                                              cTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa               >  2:137690/1‑130 (MQ=255)
                                                                                               tATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat              >  1:317704/1‑129 (MQ=255)
                                                                                                aTCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata             >  2:143909/1‑128 (MQ=255)
                                                                                                   tctcAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGAc           >  1:368849/1‑141 (MQ=255)
                                                                                                            tC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg  <  1:127105/141‑1 (MQ=255)
                                                                                                                        gCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                        <  1:152406/60‑1 (MQ=255)
                                                                                                                        gCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                        >  2:152406/1‑60 (MQ=255)
                                                                                                                             gaggagGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttt                                                    <  2:163449/75‑1 (MQ=255)
                                                                                                                             gaggagGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttt                                                    >  1:163449/1‑75 (MQ=255)
                                                                                                                                 aggTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                          <  2:414565/97‑1 (MQ=255)
                                                                                                                                 aggTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCaa                          >  1:414565/1‑97 (MQ=255)
                                                                                                                                     gtA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                >  1:387433/1‑71 (MQ=255)
                                                                                                                                     gtA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                <  2:387433/71‑1 (MQ=255)
                                                                                                                                      |                                                                                                                      
GTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCG  >  CP000730/1929728‑1929976

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