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

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG  >  NC_000913/2001764‑2002053
                                                                                                                                                   |                                                                                                                                               
ccTTCAGCGGTATAGAGTGAATTCTCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaa                                                                                                                                              <  2:301300/151‑1 (MQ=255)
ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaa                                                                                                                                              <  1:431869/151‑1 (MQ=255)
  ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                             >  1:459049/1‑150 (MQ=255)
   tCAGCGGTATAGAGTTAATTCACGATAAACAGACCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGTGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                             <  2:891640/149‑1 (MQ=255)
        ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑cccc                                                                                                                                     >  1:1260180/1‑151 (MQ=255)
         gTATAGAGTGAATTCACGATAAACATCCCTGCGTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGTGTTGCAGGGAAATTCCAGGCAGAAAAAAT‑cccc                                                                                                                                     <  2:594644/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCccgcc                                                                                                                                  <  1:637045/151‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCccgc                                                                                                                                   <  1:1268073/150‑1 (MQ=255)
              gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCgg                                                                                                                                >  1:1019385/1‑150 (MQ=255)
              gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGt                                                                                                                               >  1:803592/1‑151 (MQ=255)
                 tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGc                                                                                                                            >  1:77749/1‑151 (MQ=255)
                      tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                        >  1:555986/1‑150 (MQ=255)
                       cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                       >  2:99880/1‑150 (MQ=255)
                          gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGc                                                                                                                    >  2:487104/1‑150 (MQ=255)
                             aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTTGGGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                <  1:126696/151‑1 (MQ=255)
                              aaCAGCCTGGGGTTATATGAGTTGTGGGATTGATTACCTGTTTGTGGAGGTTTTTAAAGGGGCGGATTGGTGGGTGAAAGGGGGGGTTGTTGGGGGGTACCGGGAAAATCCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                               <  2:388358/151‑1 (MQ=255)
                                   ccttGCTTTATATGATTTTGCGGCATGTTTATCCGTTTGTGCGGGTTTTTTAATCGGACGATGAGGGGGTGAAAGGGGGGGTTTTTTGGGGGTTACAGTAAAATTCCGGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                           <  2:176451/147‑1 (MQ=255)
                                   cccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                             >  2:848619/1‑148 (MQ=255)
                                    ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                         >  1:233695/1‑151 (MQ=255)
                                        cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGGGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  2:176444/1‑150 (MQ=255)
                                          ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                    >  1:270471/1‑150 (MQ=255)
                                           tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTCTTTGGGGGTTACAGGTAAATTGCAGGCAGAAAAAAA‑CCCCGCCGGGGGCGGGGCAGCACGTTGCTGACAAATTg                                                                                                   >  2:1254803/1‑150 (MQ=255)
                                                gatttATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAAGGAGGGGTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                              <  1:377274/147‑1 (MQ=255)
                                                 aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                 >  2:684018/1‑146 (MQ=255)
                                                   ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                          >  1:865160/1‑151 (MQ=255)
                                                        gccATGATTATCCTTTTCTGAAGGGTTTTTAAGCGTACTATTAGTGGTTAAAATGGGGGGTTTTTTGGGGGTAACGGTTAATTACCAGCAAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                       <  2:1041795/147‑1 (MQ=255)
                                                          cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                   >  1:684469/1‑151 (MQ=255)
                                                            tGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAATTCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                   <  1:107742/149‑1 (MQ=255)
                                                            tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                 >  1:1087705/1‑151 (MQ=255)
                                                               tttcccGTTGTTGCGGGTTTTTTAATGAAAAGATTGGTGGGTAAAATGGGGGGTTATTTGGGGTGTAAAGGTAAATTCAAGGAAAAAAAAAG‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGag                                                                                 <  2:802722/144‑2 (MQ=255)
                                                                tatcgttttCTGCAGGTTTTTTAATCGAACTATTTGTGTGTGAAATGAGGGTTTATTTGGGGGTTACAGGAAAATTCCGGAAAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTGTATGTTGCCGGAt                                                                                 <  2:1198356/141‑1 (MQ=255)
                                                                    cGTTTCTGCAGGGTTTTTAATCGGCCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCCGAAAAAAA‑CCCCGCCGGGGGCGGGGAAGCACGTTGCTGACCAATTGCGCTTTATGTTGCCGGATgcggcg                                                                           >  2:189589/1‑149 (MQ=255)
                                                                     gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                         >  2:692772/1‑150 (MQ=255)
                                                                       ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                      >  1:1093874/1‑151 (MQ=255)
                                                                        tCTGCGGGTTTTTAAATCGGACGTTTAGTGGGTGAAATGGGGGGTTATTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                      <  1:32642/150‑1 (MQ=255)
                                                                           gcgggttttttAATCGGACGATTAGGGGGTGAAAGGAGGGGTTATTGGGGGTTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                    <  1:442994/143‑1 (MQ=255)
                                                                               ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                              >  1:830993/1‑151 (MQ=255)
                                                                                      aatggggagtttggTGGGTAAAATGAGGGTTTTTTTGGGGGTTAAGGGTAATTTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                          <  1:1112954/135‑1 (MQ=37)
                                                                                          ggaggaTTGTTGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                        <  2:915566/146‑1 (MQ=38)
                                                                                             ggtttaaggggAGAAATGGGGGGTTTTTTGGGGTTTAAAGGTATATTCCAGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                 <  2:270662/142‑1 (MQ=25)
                                                                                             cGATTAGGGGGAGAAATGAGGGTTTATTGGGGGTTTACAGGTAAATTCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                <  1:354265/151‑1 (MQ=38)
                                                                                              gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                               >  1:1208878/1‑151 (MQ=255)
                                                                                              gATTAGTGGGTAAAGTGGGGGGTTATTTGGGGGTTACGGGTAAATTCCGGGCAGAAAAAAC‑CCCCGCCGGTGTCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACaaaa                                                  <  2:925571/148‑1 (MQ=38)
                                                                                                ttggTGGGTGAAATGAGGGTTTATTGGGGGTTTGCGGGTAATTCCCAGGCAAAAAAAAA‑CCCCGCCGGTGGCGGTGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAAAAAAAt                                                <  2:821709/145‑1 (MQ=37)
                                                                                                ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                             >  1:80652/1‑151 (MQ=255)
                                                                                                  aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                            >  2:949879/1‑150 (MQ=255)
                                                                                                   ggggggaaaaTGGGGGGTGATTGGGGGTTGACGGTAAAATTCCACGCAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTTCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                           <  1:1005836/143‑1 (MQ=25)
                                                                                                   gggggTAAAAGGAGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                          <  1:382206/149‑1 (MQ=37)
                                                                                                     gggTGAAATGAGGGGTTTTTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                           <  1:271253/148‑1 (MQ=255)
                                                                                                          aaaggggggtTTATTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                   <  1:792805/145‑1 (MQ=37)
                                                                                                          aaaTGAGGGTTTTTTTGGGGGTTACAGGTAATTTCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                      <  2:110279/148‑1 (MQ=37)
                                                                                                            aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                  >  1:340868/1‑150 (MQ=255)
                                                                                                             tGAGGGGATATTTGGGGGTTACAGGTAAATTCCAGGCAGAACAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACCAATTGCGCTTTATGTTGCCGGATGAGGCGTACACGCCTTATCCGCCCTACAAAAATGTGCAAATTCaa                                    >  2:897877/1‑147 (MQ=255)
                                                                                                             ggAGGGTTTTTTTGGGGGTTAAAGGTAATTCCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                <  2:1119576/150‑1 (MQ=25)
                                                                                                              gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                               >  2:958609/1‑151 (MQ=255)
                                                                                                                 ggtttATTGGGGGTTTACAGAAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                 <  1:136220/143‑1 (MQ=25)
                                                                                                                 gggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                            >  2:1218019/1‑151 (MQ=255)
                                                                                                                   gTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                           >  2:395859/1‑150 (MQ=255)
                                                                                                                    ttATTTGGGGTTTACGGGAAAATCCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                            <  2:138556/148‑1 (MQ=37)
                                                                                                                    ttATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAACA‑CCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGAGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACACCAATGTGCCAATTCAATAACTTGCAAt                         >  2:581945/1‑151 (MQ=38)
                                                                                                                       tttGGGGGTTCCAGGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                        <  2:488922/149‑1 (MQ=255)
                                                                                                                       tggggggtTTACAGGTAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                        <  1:684502/146‑1 (MQ=37)
                                                                                                                          gggggTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                   >  1:889579/1‑151 (MQ=255)
                                                                                                                             ggTTACAGTAAAATTCCAGAAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                <  1:1265423/151‑1 (MQ=37)
                                                                                                                                 aCAGTTAAGTTCCAGG‑AGAAATACGCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCTGCGTAAACGCCTTATCCGGTCTCCAAAAATGTGCAAATTCAATAAAGTGCAATGCAACTTGTAgga            <  2:243183/151‑2 (MQ=255)
                                                                                                                                    ggTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg         >  1:616949/1‑151 (MQ=255)
                                                                                                                                    ggTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg         >  1:40379/1‑151 (MQ=255)
                                                                                                                                    ggTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt          <  2:1177659/150‑1 (MQ=255)
                                                                                                                                      tAAATTCCAGGCAGAAAAACA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt       >  2:1267961/1‑151 (MQ=255)
                                                                                                                                       aaaTTCCAGGCAGAAACAAA‑CCCCGCCGGTGGCGGGGAAGCACGATGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAGGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa        >  2:407945/1‑149 (MQ=255)
                                                                                                                                       aaaTTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt          <  1:487523/147‑1 (MQ=255)
                                                                                                                                       aaaTTCCAGCCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg              <  1:81474/143‑1 (MQ=255)
                                                                                                                                        aaTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa     <  1:849164/151‑1 (MQ=255)
                                                                                                                                           tCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg  >  2:1015093/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG  >  NC_000913/2001764‑2002053

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