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

GGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001772‑2002067
                                                                                                                                                     |                                                                                                                                                      
ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAG‑cccc                                                                                                                                                   <  1:1144854/151‑1 (MQ=255)
   ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCATTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCccgcc                                                                                                                                                >  2:507737/1‑151 (MQ=255)
     agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCccgcc                                                                                                                                                <  1:377301/149‑1 (MQ=255)
     agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCgg                                                                                                                                              >  1:235165/1‑151 (MQ=255)
         tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTgg                                                                                                                                           >  1:255994/1‑150 (MQ=255)
             ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                                      >  1:767005/1‑151 (MQ=255)
               caagaaaaaCAGCCCTGCGTTATATGAGTTATGGAAATGTTTATCCTTTTCGGGGGGTTTTTT‑AATCGAACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGT‑TTAAAGGAAAATTCCAGAA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGaa                                                                                                                                    <  2:1250183/145‑1 (MQ=255)
                      aaCAGCCTTGCGTTATATAATTTATCGCCATAATTATCCGTTTCTG‑AGGGTTTTTAAACCGGACAATTGGTGGTTAAAAGGAGGTTTTATTGGGGGT‑TTACAGGTAATTCCCAGGC‑AAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                             <  2:1032983/151‑1 (MQ=255)
                        cAGCCGGGCTTTATAGGATTAATCGGAATGATTACCCGTTTCTCGAGGGTTTTA‑AACCGGACGATAAGGGGGAAAAAGGAGGGGTTTTTGGGGGG‑TTGCGGGTAATTCCCAGCC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                            <  1:145065/150‑1 (MQ=255)
                        cAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                                           >  1:504960/1‑151 (MQ=255)
                         aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                                          >  2:174503/1‑151 (MQ=255)
                                    atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                               >  1:1048524/1‑151 (MQ=255)
                                         aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                           >  2:797590/1‑150 (MQ=255)
                                           ttaccggagtGATTATCGTTTTCTGCAGGTTTTTT‑AATGGGATGTTTGGGGGGTAAAATGGGGGGTTTTTGGGGGT‑TTACGGGAAAATCCAAGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                         <  2:538792/141‑1 (MQ=255)
                                             accGGAATGATTATCCGTTTCTGCGGTTTTTTA‑AATCGAACGATTGGGGGGTAAAAGGGGGGTTTATTGGGGGT‑TAACAGGAAAATCCCAGGC‑AAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGATTTATg                                                                                                        <  2:332539/147‑1 (MQ=255)
                                                ggCATTTTTTTCCTTTTCTGCGGGTTTTTT‑AATGGGACGATTGGGGGGTAAAAGGGGGGTTTTTTTGGGGG‑TAACAGGAAAATCCAAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                     <  1:1085306/149‑1 (MQ=255)
                                                 gaattaTTATCCTTTTCTGCAGGTTTTTA‑AATCGGAGAATTAGTGGGTAAAAGGGGGGTTTTTTTGGGGG‑TTACAGAAAAATCCAAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                  <  1:1234714/146‑1 (MQ=255)
                                                   aTGATTATCCTTTTCGGCAGGTTTTTT‑AAGCGGACGATTAGTGGGGAAAATGAGGGTTTATTGGGGGG‑TTACAGGAAAATCCCAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                                <  2:1222839/151‑1 (MQ=255)
                                                    tGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGACATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                               >  2:1005476/1‑151 (MQ=255)
                                                          tCCGTTTCTGCAGGTTTTTT‑AATCGGACGATTGGGGGGTGAAATGGGGGGTTTTTTGGGGG‑TTACAGGTAAATCCCAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                           <  1:349589/149‑1 (MQ=255)
                                                                   gCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGATGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                >  1:1207339/1‑151 (MQ=255)
                                                                    cAGGTTTTTA‑AATCGGACGTTTAGTGGGTAAAATG‑GGGGTTTTTTGGGGGTTTACAGGTAAATTCCAGGA‑AGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                                   <  2:352990/147‑1 (MQ=255)
                                                                    cAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  1:702773/1‑151 (MQ=255)
                                                                     aGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  2:805627/1‑150 (MQ=255)
                                                                         ttttt‑AATGGGAGGTTTAGTGGGGAAATTGGGGGTTTTTTTGGGGT‑TTACAGTTAATCCCCTGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                           <  2:780876/150‑1 (MQ=255)
                                                                         tttta‑AATGGGATGATGGGTGGGTGAAAGGAGGGTTTTTTGGGGGT‑TTACAGGAAAATTCAGGGC‑AAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                           <  2:982087/150‑1 (MQ=255)
                                                                               aatggaaCGATTGGTGGTTGAAAGGGGGGTTTTTTTGGGGT‑TTACAGTTAAATCCCAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                       <  1:671564/143‑1 (MQ=25)
                                                                                   gaaCGTTTAGGGGGTGAAATGGGGGGTTATTTGGGGG‑TTAGAGAT‑AATTCCAGGCAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                   <  1:368750/147‑1 (MQ=37)
                                                                                     aCGATTGGTGGGTAAAAGGAGGGTTTTTTTGGGGG‑TGAAAGGAAAATCCCAGAA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                                <  1:980995/150‑1 (MQ=25)
                                                                                         tttgTGGGTAAAAGGGGGGTTTATTTGGGGG‑TTACAGGTAACTCCCAGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                            <  2:1086440/147‑1 (MQ=37)
                                                                                         ttAGGGGGTAAAAGGGGGGTTTATTGGGGGT‑TTACAGGTAATTCCCAGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                              <  2:621663/148‑1 (MQ=37)
                                                                                         tgaGTGGGGGAAGGGGGGGGTTTTTGGGGAG‑TTAAAGTTAAATTCCGGAA‑AGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                            <  2:1094014/148‑1 (MQ=25)
                                                                                            gTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                         >  1:152504/1‑150 (MQ=255)
                                                                                              gggTAAAATGGGGGTTTTTTTGGGGG‑TTAAGGAAAAATCCCAGGA‑AAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                       <  2:149143/150‑1 (MQ=25)
                                                                                                gTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                     >  1:1244537/1‑150 (MQ=255)
                                                                                                  gAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   >  1:767100/1‑150 (MQ=255)
                                                                                                   acaTGAGGGGTTATTTGGGGG‑TTACAGGTAAATCCCAGGC‑AGAAAAAAA‑CCGCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGGTGCGGCGTAAACGCCTTATCCGGCCTCCAAAAATGTGCAAATTCaa                                                  >  1:432899/3‑150 (MQ=255)
                                                                                                   aaagaaGGGTTTATTTGGGGG‑TGAGAGGAAAATTCAGGGA‑AAAAAAGAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                    <  2:802570/143‑1 (MQ=25)
                                                                                                          gggTTTTTTGGGGGTTTAGAGAAAACTCCCAGGC‑AAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                             <  2:775801/149‑1 (MQ=25)
                                                                                                             ttttttGGGGG‑TAACAGGAAATTCCCAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                        <  1:965123/147‑1 (MQ=25)
                                                                                                              tATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                      >  2:136910/1‑151 (MQ=255)
                                                                                                                 ttGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                    >  2:150575/1‑150 (MQ=255)
                                                                                                                  tGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                  >  1:496562/1‑151 (MQ=255)
                                                                                                                      gg‑TTACAGTAAAATCCCAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                              <  1:174500/151‑1 (MQ=37)
                                                                                                                      gg‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGGACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGCAAATTCAATAAATTGCAATTCAACTTg                               >  1:414134/1‑150 (MQ=255)
                                                                                                                       ggtttacagaaaaTTCCCGGG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                               <  2:634811/141‑1 (MQ=25)
                                                                                                                        ttaaaagaaaaaTCCCAGGA‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAGGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGAAAATTCAATAAAGGGCAATTCAACTGGt                              <  1:1071619/141‑1 (MQ=25)
                                                                                                                          tACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAACTTGCAATTCAACTTGTAgg                           >  2:897177/1‑151 (MQ=255)
                                                                                                                          tACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCAGCCGGTGGCGGGGAAGCACGTTGCTGCCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTgaa                             >  1:1160192/1‑147 (MQ=255)
                                                                                                                           acgggTAATTTCCAGGA‑AGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                             <  1:1110131/145‑1 (MQ=255)
                                                                                                                            cAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                         >  2:621152/1‑151 (MQ=255)
                                                                                                                                 aattcccAGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                      <  1:344643/144‑1 (MQ=255)
                                                                                                                                 aaaTCCCGGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                      <  1:417341/149‑1 (MQ=255)
                                                                                                                                 aaaTCCCAGGC‑AAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                     <  1:797555/150‑1 (MQ=255)
                                                                                                                                    ttCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc                 >  1:1161624/1‑151 (MQ=255)
                                                                                                                                       ccggcaaaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg                <  1:289699/142‑1 (MQ=255)
                                                                                                                                          ggaaaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc                 <  2:142971/141‑1 (MQ=255)
                                                                                                                                             aGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc          >  1:478676/1‑150 (MQ=255)
                                                                                                                                               aaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt        <  2:1048585/150‑1 (MQ=255)
                                                                                                                                                 aaaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc       >  2:929105/1‑149 (MQ=255)
                                                                                                                                                  aaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg    <  2:557452/151‑1 (MQ=255)
                                                                                                                                                  aaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg    >  1:192962/1‑151 (MQ=255)
                                                                                                                                                     a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa  >  2:217109/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  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: ‑