Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA46.2% 94.8 / 60.7 91intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (46/3);  new base A (0/42);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 8.97e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001764‑2002059
                                                                                                                                                        |                                                                                                                                                     
ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGaaa                                                                                                                                                    <  1:985183/151‑1 (MQ=255)
  ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGaaaa                                                                                                                                                   >  1:1025848/1‑150 (MQ=255)
       cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑cc                                                                                                                                             <  1:281105/150‑1 (MQ=255)
          tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCAGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCccgc                                                                                                                                         >  2:840460/1‑151 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCccgccg                                                                                                                                       >  1:260950/1‑151 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCccgccg                                                                                                                                       >  2:1024409/1‑151 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCccgccg                                                                                                                                       >  1:924979/1‑151 (MQ=255)
              gagTGAATTCACGATAAACAGCCCTGCGTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAAAAAAAAC‑CCccgccg                                                                                                                                       <  2:433974/149‑1 (MQ=255)
               agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATACGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTg                                                                                                                                    >  1:1257283/1‑151 (MQ=255)
                     ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                              >  2:477262/1‑151 (MQ=255)
                       cacaaaaaccAGCCTGGCGTTATATGAGTTATGGGCATGATTTTCCGTTTCTGCGGGGTTTTTAATCGGACGTTTGGTGGGTGAAATGGGGGGTTTTTTGGGGG‑TT‑A‑C‑AGGT‑AAATCCCAGGAAAAAAAAAA‑CCCCGCCGGTGGCGGGGaa                                                                                                                            <  2:568281/142‑1 (MQ=255)
                        aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                             >  1:513374/1‑149 (MQ=255)
                           aaaaaCGGCCCTGCGTTATATGATTTATCGGCATGATTATCCTTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTGAAATGGGGGGTTATTGGGGGG‑TT‑A‑C‑AGGT‑AAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                        <  1:602848/149‑1 (MQ=255)
                           aTAACCACCCCTGGTTTATATGAGTTATCGGTATAATTATCCGTTTCCGATGGGTTTTAAATCGGATATTTGGTGGGTGAAAAGAGGGGTT‑TTTGGGGG‑TT‑T‑CAAGGT‑AAATCCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                           <  1:291236/148‑1 (MQ=255)
                             aaaCAGCCCGGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGAGGATTAGGGGGTAAAAGGGGGGGTTATTGGGGGG‑TT‑A‑A‑AGGA‑AAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACg                                                                                                                       <  2:378106/150‑1 (MQ=255)
                             aaaCACCCCTGCGTTATATGATTTATCGCATTTATTATCTGTTTCTGCAGGGTTTTAATTCGGACGATAAGTGGGTAAAAGGGGGGGTTATTGGGGGT‑TT‑A‑A‑AGGT‑AAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCa                                                                                                                         <  2:850674/148‑1 (MQ=255)
                                 aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCTGGTGGCGGGGAAGCACGTTGCt                                                                                                                  >  2:1364336/1‑151 (MQ=255)
                                       gctttATATGATTTATCGAAATAATTATCCTTTTCGGCAGGTTTTTTAATCGGACGTTTAGTGGGTAAAAGGGGGGTTTATTGGGGGG‑TT‑A‑C‑AGGA‑AAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                              <  1:1142034/146‑1 (MQ=255)
                                         gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                          >  1:596634/1‑151 (MQ=255)
                                          ttATATGAGTTATCGGCATGTTTTTCCTTTTCTGCAGGTTTTTTAATCGGAGGATTAGTGGGGAAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGA‑AAATTCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                            <  1:1113609/148‑1 (MQ=255)
                                           tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                        >  2:341337/1‑151 (MQ=255)
                                             tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                       >  1:1024081/1‑150 (MQ=255)
                                                gAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                      <  1:884256/148‑1 (MQ=255)
                                                gAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   >  2:134138/1‑151 (MQ=255)
                                                  gTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTGAAAGGGGGGTTTATTGGGGGG‑TT‑A‑C‑AGGA‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   <  2:375418/149‑1 (MQ=255)
                                                   ttttCGGAATGATTATCCTTTTCTGCGGGGTTTTTAAGCGAACAATTAGGGGGGAAAATGGGGGTTTATTGGGGGT‑TT‑A‑C‑AGGT‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   <  1:1354048/145‑1 (MQ=255)
                                                   ttATCGGAATGATTATCCTTTTCGGCGGGGTTTTAAATCGGACGATTGGTGGGTAAAATGGGGGTTTATTGGGGGG‑TT‑A‑C‑AGGA‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                  <  1:792538/149‑1 (MQ=255)
                                                    tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                               >  1:181863/1‑151 (MQ=255)
                                                      tCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                             >  1:1190337/1‑151 (MQ=255)
                                                         gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                          >  1:855354/1‑151 (MQ=255)
                                                           aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  1:1150227/1‑150 (MQ=255)
                                                              aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                      >  2:991335/1‑150 (MQ=255)
                                                               ttttCCGTTTCTGCGGGTTTTTTAATCGGACGTTTGGTGGGTAAAAGGGGGGTTTTTTTGGGGG‑TT‑A‑C‑AGTT‑AATTTCCAGCAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                     <  2:931114/147‑1 (MQ=255)
                                                               ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                    >  2:582703/1‑151 (MQ=255)
                                                                 aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                   >  2:1282828/1‑150 (MQ=255)
                                                                 aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                   >  2:814957/1‑150 (MQ=255)
                                                                   ccGTTTCTGCGGGTTTTTTAATCGGACGATTGGTGGGTGAAATGAGGGTTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                 <  2:569233/150‑1 (MQ=255)
                                                                   ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                >  1:1278609/1‑151 (MQ=255)
                                                                      tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                              >  1:420484/1‑150 (MQ=255)
                                                                          tGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGcc                                                                         >  1:688273/1‑151 (MQ=255)
                                                                             aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       >  2:295058/1‑150 (MQ=255)
                                                                               ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                    >  2:700017/1‑151 (MQ=255)
                                                                                 tttttAAGGGGAAGATAAGGGGTTAAAAGGAGGTGTTATTGGGGGT‑TA‑A‑C‑TGAT‑AATTTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       <  2:1063872/146‑1 (MQ=255)
                                                                                  ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                 >  2:780373/1‑151 (MQ=255)
                                                                                   tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                >  2:1379426/1‑151 (MQ=255)
                                                                                          ggaggataggTGGGTAAATTGAGGGTTTATGGGGGGG‑TT‑A‑A‑GGGT‑AAA‑TCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                           <  1:633409/140‑1 (MQ=25)
                                                                                          ggaggaTTGGGGGGTAAAATGGGGGGTT‑TTTGGGGG‑TT‑A‑A‑AAGTAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  2:1403056/150‑1 (MQ=25)
                                                                                          ggACGATTGGGGGGTGAAATGGGGGTTTTTTGGGGGG‑TT‑A‑C‑GGGA‑AAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  2:1142738/150‑1 (MQ=37)
                                                                                          ggACGATTAGGGGGTGAAATGAGGGGTTTTTTGGGGG‑TT‑A‑C‑AGGT‑AAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  2:395826/150‑1 (MQ=255)
                                                                                           ggacgatggggggTGAAAATGGGGGTTTTTTGGGGG‑TT‑T‑A‑AAGT‑AAAATCCAGGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                           <  2:1180403/139‑1 (MQ=25)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  2:1039451/1‑150 (MQ=255)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  2:815371/1‑150 (MQ=255)
                                                                                              gATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     >  2:166958/1‑151 (MQ=255)
                                                                                               aTTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     >  1:1207380/1‑150 (MQ=255)
                                                                                                ttAGTGGGGGAAAGGGGGGGTTTTTGGGGGT‑TT‑A‑C‑AGGA‑AAATCCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                    <  1:638426/150‑1 (MQ=25)
                                                                                                ttAGGGGGTAAAAGGGGGGTTTTTTGGGGGG‑TT‑A‑C‑AGAT‑AATTCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                      <  2:39351/148‑1 (MQ=25)
                                                                                                  aGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  1:783240/1‑150 (MQ=255)
                                                                                                   gggggTAAAATGGGGGTTTTTTGGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                   <  1:75998/146‑1 (MQ=37)
                                                                                                    gggggaaaaaggggggTTTTTTGGGGG‑TT‑AGC‑AGGA‑AAATCCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     <  1:1282856/134‑1 (MQ=25)
                                                                                                    ggggaaaaaTGGGGGGTTATTGGGGGT‑TT‑A‑C‑AGGA‑AAATTCGAGCGAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACAAAAAtgtg                                                   <  1:850950/141‑1 (MQ=37)
                                                                                                     gggaaaaaTGAGGGTTTATTGGGGGG‑TT‑A‑C‑AGGT‑AAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                               <  1:436352/145‑1 (MQ=38)
                                                                                                      ggTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                             >  2:1207073/1‑151 (MQ=255)
                                                                                                      ggTGAAAAGGGGGGTTTTTTGGGGGTTT‑A‑A‑AGTA‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  <  2:662362/147‑1 (MQ=25)
                                                                                                        tGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                            >  1:483701/1‑150 (MQ=255)
                                                                                                           aaaggggggGTATTTTGGGG‑TTAA‑C‑AGTA‑AAATCCAAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                         <  2:138356/146‑1 (MQ=25)
                                                                                                              gAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                     >  2:362132/1‑151 (MQ=255)
                                                                                                                gggtTTTTTTGGGGT‑TA‑A‑C‑AGGT‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                   <  2:924958/151‑1 (MQ=25)
                                                                                                                ggggTTATTTGGGGG‑TT‑A‑C‑AGGT‑AATTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                    <  2:762434/150‑1 (MQ=255)
                                                                                                                 gggTTTTTTGGGGG‑TT‑A‑C‑AGGA‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  <  1:352210/151‑1 (MQ=37)
                                                                                                                  ggTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  >  2:924465/1‑150 (MQ=255)
                                                                                                                    ttATTTGGGGT‑TT‑A‑C‑AGGT‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                 <  1:796370/149‑1 (MQ=255)
                                                                                                                     tttttGGGGGTTT‑A‑C‑AGGT‑AAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                <  1:1238714/148‑1 (MQ=37)
                                                                                                                     tATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                              >  2:954341/1‑151 (MQ=255)
                                                                                                                       tttGGGGT‑TT‑A‑C‑ATGT‑AAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  1:1367581/150‑1 (MQ=255)
                                                                                                                       ttgggggg‑TT‑A‑C‑AGGT‑AAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  1:367280/147‑1 (MQ=255)
                                                                                                                       ttgggggg‑TT‑A‑A‑AGGA‑AAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                             <  2:244217/147‑1 (MQ=37)
                                                                                                                              g‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                      >  1:370268/1‑150 (MQ=255)
                                                                                                                                 t‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   >  1:1117485/1‑151 (MQ=255)
                                                                                                                                 t‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   >  2:143227/1‑151 (MQ=255)
                                                                                                                                 t‑A‑C‑AGGT‑AAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                    <  1:134144/150‑1 (MQ=255)
                                                                                                                                            aattcccAGCCAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  2:1025827/145‑1 (MQ=25)
                                                                                                                                            aaaTTCCAGGAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  1:1176619/150‑1 (MQ=255)
                                                                                                                                            aaaTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             <  1:1132009/150‑1 (MQ=255)
                                                                                                                                             aaatCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                    <  2:260943/140‑1 (MQ=255)
                                                                                                                                             aaTTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa            >  2:1395926/1‑150 (MQ=255)
                                                                                                                                             aaTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa           <  1:991936/151‑1 (MQ=255)
                                                                                                                                               ttCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc         >  1:978365/1‑151 (MQ=255)
                                                                                                                                                 caaGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg        <  1:1395011/148‑1 (MQ=255)
                                                                                                                                                  cAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc       >  1:517350/1‑150 (MQ=255)
                                                                                                                                                      cAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  2:393621/1‑151 (MQ=255)
                                                                                                                                                      cAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  1:66690/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TT‑A‑C‑AGGT‑AAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001764‑2002059

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