Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8480CA61.1% 32.7 / 91.8 72intergenic (‑931/‑237)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (12/16);  new base A (27/17);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCG  >  CP000730/1929710‑1929976
                                                                                                                                           |                                                                                                                                   
gAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCag                                                                                                                                   <  2:255884/141‑1 (MQ=255)
      gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAgagga                                                                                                                             <  2:203942/141‑1 (MQ=255)
                  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:207455/135‑4 (MQ=255)
                                                                    tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                           <  1:207285/135‑4 (MQ=255)
                                                                     gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                           <  2:93420/141‑1 (MQ=255)
                                                                     gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                           <  2:183901/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                                                 <  1:405295/124‑6 (MQ=255)
                                                                                            tggtataATAATAGAGTCGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                 >  2:405295/4‑122 (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                                                                   <  1:54759/92‑1 (MQ=255)
                                                                                                             cGCCTATCTCTCAGGCGTC‑AATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                   >  2:54759/1‑92 (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                                                                        >  2:152406/1‑60 (MQ=255)
                                                                                                                                          gCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                        <  1:152406/60‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
GAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTC‑AATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCG  >  CP000730/1929710‑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: ‑