Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA44.8% 50.8 / 40.8 58intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (27/5);  new base A (0/26);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 8.16e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001764‑2002059
                                                                                                                                                      |                                                                                                                                                     
ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGaa                                                                                                                                                     >  2:1071681/1‑150 (MQ=255)
 cTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGaaa                                                                                                                                                    >  2:991570/1‑150 (MQ=255)
        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)
                                                                                                                                                     aGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  1:478676/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGC‑AGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001764‑2002059

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