Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC43.6% 76.5 / 12.3 56intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (25/6);  new base C (0/24);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 2.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001772‑2002067
                                                                                                                                                     |                                                                                                                                                      
ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAGcccc                                                                                                                                                   <  1:1144854/151‑1 (MQ=255)
   ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCATTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                                >  2:507737/1‑151 (MQ=255)
     agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                                <  1:377301/149‑1 (MQ=255)
     agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                              >  1:235165/1‑151 (MQ=255)
         tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTgg                                                                                                                                           >  1:255994/1‑150 (MQ=255)
             ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                      >  1:767005/1‑151 (MQ=255)
               caagaaaaaCAGCCCTGCGTTATATGAGTTATGGAAATGTTTATCCTTTTCGGGGGGTTTTTTAATCGAACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGTTTAAAGGAAAATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGaa                                                                                                                                    <  2:1250183/145‑1 (MQ=255)
                      aaCAGCCTTGCGTTATATAATTTATCGCCATAATTATCCGTTTCTGAGGGTTTTTAAACCGGACAATTGGTGGTTAAAAGGAGGTTTTATTGGGGGTTTACAGGTAATTCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                             <  2:1032983/151‑1 (MQ=255)
                        cAGCCGGGCTTTATAGGATTAATCGGAATGATTACCCGTTTCTCGAGGGTTTTAAACCGGACGATAAGGGGGAAAAAGGAGGGGTTTTTGGGGGGTTGCGGGTAATTCCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                            <  1:145065/150‑1 (MQ=255)
                        cAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                           >  1:504960/1‑151 (MQ=255)
                         aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                          >  2:174503/1‑151 (MQ=255)
                                    atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                               >  1:1048524/1‑151 (MQ=255)
                                         aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                           >  2:797590/1‑150 (MQ=255)
                                           ttaccggagtGATTATCGTTTTCTGCAGGTTTTTTAATGGGATGTTTGGGGGGTAAAATGGGGGGTTTTTGGGGGTTTACGGGAAAATCCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                         <  2:538792/141‑1 (MQ=255)
                                             accGGAATGATTATCCGTTTCTGCGGTTTTTTAAATCGAACGATTGGGGGGTAAAAGGGGGGTTTATTGGGGGTTAACAGGAAAATCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGATTTATg                                                                                                        <  2:332539/147‑1 (MQ=255)
                                                ggCATTTTTTTCCTTTTCTGCGGGTTTTTTAATGGGACGATTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTAACAGGAAAATCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                     <  1:1085306/149‑1 (MQ=255)
                                                 gaattaTTATCCTTTTCTGCAGGTTTTTAAATCGGAGAATTAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGAAAAATCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                  <  1:1234714/146‑1 (MQ=255)
                                                   aTGATTATCCTTTTCGGCAGGTTTTTTAAGCGGACGATTAGTGGGGAAAATGAGGGTTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                                <  2:1222839/151‑1 (MQ=255)
                                                    tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGACATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                               >  2:1005476/1‑151 (MQ=255)
                                                          tCCGTTTCTGCAGGTTTTTTAATCGGACGATTGGGGGGTGAAATGGGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                           <  1:349589/149‑1 (MQ=255)
                                                                   gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGATGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                >  1:1207339/1‑151 (MQ=255)
                                                                    cAGGTTTTTAAATCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTGGGGGTTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                   <  2:352990/147‑1 (MQ=255)
                                                                    cAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  1:702773/1‑151 (MQ=255)
                                                                     aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  2:805627/1‑150 (MQ=255)
                                                                         tttttAATGGGAGGTTTAGTGGGGAAATTGGGGGTTTTTTTGGGGTTTACAGTTAATCCCCTGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                           <  2:780876/150‑1 (MQ=255)
                                                                         ttttaAATGGGATGATGGGTGGGTGAAAGGAGGGTTTTTTGGGGGTTTACAGGAAAATTCAGGGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                           <  2:982087/150‑1 (MQ=255)
                                                                               aatggaaCGATTGGTGGTTGAAAGGGGGGTTTTTTTGGGGTTTACAGTTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                       <  1:671564/143‑1 (MQ=25)
                                                                                   gaaCGTTTAGGGGGTGAAATGGGGGGTTATTTGGGGGTTAGAGATAATTCCAGGCAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                   <  1:368750/147‑1 (MQ=37)
                                                                                     aCGATTGGTGGGTAAAAGGAGGGTTTTTTTGGGGGTGAAAGGAAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                                <  1:980995/150‑1 (MQ=25)
                                                                                         tttgTGGGTAAAAGGGGGGTTTATTTGGGGGTTACAGGTAACTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                            <  2:1086440/147‑1 (MQ=37)
                                                                                         ttAGGGGGTAAAAGGGGGGTTTATTGGGGGTTTACAGGTAATTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                              <  2:621663/148‑1 (MQ=37)
                                                                                         tgaGTGGGGGAAGGGGGGGGTTTTTGGGGAGTTAAAGTTAAATTCCGGAAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                            <  2:1094014/148‑1 (MQ=25)
                                                                                            gTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                         >  1:152504/1‑150 (MQ=255)
                                                                                              gggTAAAATGGGGGTTTTTTTGGGGGTTAAGGAAAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                       <  2:149143/150‑1 (MQ=25)
                                                                                                gTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                     >  1:1244537/1‑150 (MQ=255)
                                                                                                  gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   >  1:767100/1‑150 (MQ=255)
                                                                                                   acaTGAGGGGTTATTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAAACCGCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGGTGCGGCGTAAACGCCTTATCCGGCCTCCAAAAATGTGCAAATTCaa                                                  >  1:432899/3‑150 (MQ=255)
                                                                                                   aaagaaGGGTTTATTTGGGGGTGAGAGGAAAATTCAGGGAAAAAAAGACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                    <  2:802570/143‑1 (MQ=25)
                                                                                                          gggTTTTTTGGGGGTTTAGAGAAAACTCCCAGGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                             <  2:775801/149‑1 (MQ=25)
                                                                                                             ttttttGGGGGTAACAGGAAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                        <  1:965123/147‑1 (MQ=25)
                                                                                                              tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                      >  2:136910/1‑151 (MQ=255)
                                                                                                                 ttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                    >  2:150575/1‑150 (MQ=255)
                                                                                                                  tGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                  >  1:496562/1‑151 (MQ=255)
                                                                                                                      ggTTACAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                              <  1:174500/151‑1 (MQ=37)
                                                                                                                      ggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGGACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGCAAATTCAATAAATTGCAATTCAACTTg                               >  1:414134/1‑150 (MQ=255)
                                                                                                                       ggtttacagaaaaTTCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                               <  2:634811/141‑1 (MQ=25)
                                                                                                                        ttaaaagaaaaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGAAAATTCAATAAAGGGCAATTCAACTGGt                              <  1:1071619/141‑1 (MQ=25)
                                                                                                                          tACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAACTTGCAATTCAACTTGTAgg                           >  2:897177/1‑151 (MQ=255)
                                                                                                                          tACAGGTAAATTCCAGGCAGAAAAAAACCCAGCCGGTGGCGGGGAAGCACGTTGCTGCCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTgaa                             >  1:1160192/1‑147 (MQ=255)
                                                                                                                           acgggTAATTTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                             <  1:1110131/145‑1 (MQ=255)
                                                                                                                            cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                         >  2:621152/1‑151 (MQ=255)
                                                                                                                                 aattcccAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                      <  1:344643/144‑1 (MQ=255)
                                                                                                                                 aaaTCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                      <  1:417341/149‑1 (MQ=255)
                                                                                                                                 aaaTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                     <  1:797555/150‑1 (MQ=255)
                                                                                                                                    ttCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc                 >  1:1161624/1‑151 (MQ=255)
                                                                                                                                       ccggcaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg                <  1:289699/142‑1 (MQ=255)
                                                                                                                                          ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc                 <  2:142971/141‑1 (MQ=255)
                                                                                                                                             aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc          >  1:478676/1‑150 (MQ=255)
                                                                                                                                               aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt        <  2:1048585/150‑1 (MQ=255)
                                                                                                                                                 aaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc       >  2:929105/1‑149 (MQ=255)
                                                                                                                                                  aaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg    <  2:557452/151‑1 (MQ=255)
                                                                                                                                                  aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg    >  1:192962/1‑151 (MQ=255)
                                                                                                                                                     aCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa  >  2:217109/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001772‑2002067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: