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

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG  >  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)
        ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                     >  1:1260180/1‑151 (MQ=255)
         gTATAGAGTGAATTCACGATAAACATCCCTGCGTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGTGTTGCAGGGAAATTCCAGGCAGAAAAAATcccc                                                                                                                                     <  2:594644/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                  <  1:637045/151‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                   <  1:1268073/150‑1 (MQ=255)
              gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                >  1:1019385/1‑150 (MQ=255)
              gagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGt                                                                                                                               >  1:803592/1‑151 (MQ=255)
                 tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGc                                                                                                                            >  1:77749/1‑151 (MQ=255)
                      tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                        >  1:555986/1‑150 (MQ=255)
                       cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                       >  2:99880/1‑150 (MQ=255)
                          gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                    >  2:487104/1‑150 (MQ=255)
                             aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTTGGGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt                                                                                                                <  1:126696/151‑1 (MQ=255)
                              aaCAGCCTGGGGTTATATGAGTTGTGGGATTGATTACCTGTTTGTGGAGGTTTTTAAAGGGGCGGATTGGTGGGTGAAAGGGGGGGTTGTTGGGGGGTACCGGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                               <  2:388358/151‑1 (MQ=255)
                                   ccttGCTTTATATGATTTTGCGGCATGTTTATCCGTTTGTGCGGGTTTTTTAATCGGACGATGAGGGGGTGAAAGGGGGGGTTTTTTGGGGGTTACAGTAAAATTCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                           <  2:176451/147‑1 (MQ=255)
                                   cccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGc                                                                                                             >  2:848619/1‑148 (MQ=255)
                                    ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                         >  1:233695/1‑151 (MQ=255)
                                        cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACaaa                                                                                                      >  2:176444/1‑150 (MQ=255)
                                          ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                    >  1:270471/1‑150 (MQ=255)
                                           tataTGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTCTTTGGGGGTTACAGGTAAATTGCAGGCAGAAAAAAACCCCGCCGGGGGCGGGGCAGCACGTTGCTGACAAATTg                                                                                                   >  2:1254803/1‑150 (MQ=255)
                                                gatttATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAAGGAGGGGTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                              <  1:377274/147‑1 (MQ=255)
                                                 aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                 >  2:684018/1‑146 (MQ=255)
                                                   ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                          >  1:865160/1‑151 (MQ=255)
                                                        gccATGATTATCCTTTTCTGAAGGGTTTTTAAGCGTACTATTAGTGGTTAAAATGGGGGGTTTTTTGGGGGTAACGGTTAATTACCAGCAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                       <  2:1041795/147‑1 (MQ=255)
                                                          cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                   >  1:684469/1‑151 (MQ=255)
                                                            tGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGGTTACAGGAAATTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                   <  1:107742/149‑1 (MQ=255)
                                                            tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                 >  1:1087705/1‑151 (MQ=255)
                                                               tttcccGTTGTTGCGGGTTTTTTAATGAAAAGATTGGTGGGTAAAATGGGGGGTTATTTGGGGTGTAAAGGTAAATTCAAGGAAAAAAAAAGCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGag                                                                                 <  2:802722/144‑2 (MQ=255)
                                                                tatcgttttCTGCAGGTTTTTTAATCGAACTATTTGTGTGTGAAATGAGGGTTTATTTGGGGGTTACAGGAAAATTCCGGAAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTGTATGTTGCCGGAt                                                                                 <  2:1198356/141‑1 (MQ=255)
                                                                    cGTTTCTGCAGGGTTTTTAATCGGCCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCCGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACCAATTGCGCTTTATGTTGCCGGATgcggcg                                                                           >  2:189589/1‑149 (MQ=255)
                                                                     gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                         >  2:692772/1‑150 (MQ=255)
                                                                       ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                      >  1:1093874/1‑151 (MQ=255)
                                                                        tCTGCGGGTTTTTAAATCGGACGTTTAGTGGGTGAAATGGGGGGTTATTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                      <  1:32642/150‑1 (MQ=255)
                                                                           gcgggttttttAATCGGACGATTAGGGGGTGAAAGGAGGGGTTATTGGGGGTTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                    <  1:442994/143‑1 (MQ=255)
                                                                               ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                              >  1:830993/1‑151 (MQ=255)
                                                                                      aatggggagtttggTGGGTAAAATGAGGGTTTTTTTGGGGGTTAAGGGTAATTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                          <  1:1112954/135‑1 (MQ=37)
                                                                                          ggaggaTTGTTGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                        <  2:915566/146‑1 (MQ=38)
                                                                                             ggtttaaggggAGAAATGGGGGGTTTTTTGGGGTTTAAAGGTATATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                 <  2:270662/142‑1 (MQ=25)
                                                                                             cGATTAGGGGGAGAAATGAGGGTTTATTGGGGGTTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                <  1:354265/151‑1 (MQ=38)
                                                                                              gATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                               >  1:1208878/1‑151 (MQ=255)
                                                                                              gATTAGTGGGTAAAGTGGGGGGTTATTTGGGGGTTACGGGTAAATTCCGGGCAGAAAAAACCCCCGCCGGTGTCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACaaaa                                                  <  2:925571/148‑1 (MQ=38)
                                                                                                ttggTGGGTGAAATGAGGGTTTATTGGGGGTTTGCGGGTAATTCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGTGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAAAAAAAt                                                <  2:821709/145‑1 (MQ=37)
                                                                                                ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                             >  1:80652/1‑151 (MQ=255)
                                                                                                  aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                            >  2:949879/1‑150 (MQ=255)
                                                                                                   ggggggaaaaTGGGGGGTGATTGGGGGTTGACGGTAAAATTCCACGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTTCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                           <  1:1005836/143‑1 (MQ=25)
                                                                                                   gggggTAAAAGGAGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                          <  1:382206/149‑1 (MQ=37)
                                                                                                     gggTGAAATGAGGGGTTTTTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                           <  1:271253/148‑1 (MQ=255)
                                                                                                          aaaggggggtTTATTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                   <  1:792805/145‑1 (MQ=37)
                                                                                                          aaaTGAGGGTTTTTTTGGGGGTTACAGGTAATTTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                      <  2:110279/148‑1 (MQ=37)
                                                                                                            aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                  >  1:340868/1‑150 (MQ=255)
                                                                                                             tGAGGGGATATTTGGGGGTTACAGGTAAATTCCAGGCAGAACAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACCAATTGCGCTTTATGTTGCCGGATGAGGCGTACACGCCTTATCCGCCCTACAAAAATGTGCAAATTCaa                                    >  2:897877/1‑147 (MQ=255)
                                                                                                             ggAGGGTTTTTTTGGGGGTTAAAGGTAATTCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                <  2:1119576/150‑1 (MQ=25)
                                                                                                              gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                               >  2:958609/1‑151 (MQ=255)
                                                                                                                 ggtttATTGGGGGTTTACAGAAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                 <  1:136220/143‑1 (MQ=25)
                                                                                                                 gggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                            >  2:1218019/1‑151 (MQ=255)
                                                                                                                   gTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                           >  2:395859/1‑150 (MQ=255)
                                                                                                                    ttATTTGGGGTTTACGGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                            <  2:138556/148‑1 (MQ=37)
                                                                                                                    ttATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAACACCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGAGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACACCAATGTGCCAATTCAATAACTTGCAAt                         >  2:581945/1‑151 (MQ=38)
                                                                                                                       tttGGGGGTTCCAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                        <  2:488922/149‑1 (MQ=255)
                                                                                                                       tggggggtTTACAGGTAAATTCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                        <  1:684502/146‑1 (MQ=37)
                                                                                                                          gggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                   >  1:889579/1‑151 (MQ=255)
                                                                                                                             ggTTACAGTAAAATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                <  1:1265423/151‑1 (MQ=37)
                                                                                                                                 aCAGTTAAGTTCCAGGAGAAATACGCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCTGCGTAAACGCCTTATCCGGTCTCCAAAAATGTGCAAATTCAATAAAGTGCAATGCAACTTGTAgga            <  2:243183/151‑2 (MQ=255)
                                                                                                                                    ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg         >  1:616949/1‑151 (MQ=255)
                                                                                                                                    ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg         >  1:40379/1‑151 (MQ=255)
                                                                                                                                    ggTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt          <  2:1177659/150‑1 (MQ=255)
                                                                                                                                      tAAATTCCAGGCAGAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGACTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt       >  2:1267961/1‑151 (MQ=255)
                                                                                                                                       aaaTTCCAGGCAGAAACAAACCCCGCCGGTGGCGGGGAAGCACGATGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAGGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa        >  2:407945/1‑149 (MQ=255)
                                                                                                                                       aaaTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt          <  1:487523/147‑1 (MQ=255)
                                                                                                                                       aaaTTCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg              <  1:81474/143‑1 (MQ=255)
                                                                                                                                        aaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa     <  1:849164/151‑1 (MQ=255)
                                                                                                                                           tCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcg  >  2:1015093/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCG  >  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: