Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA40.2% 135.7 / 40.4 97intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (49/9);  new base A (0/39);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 8.36e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCA  >  NC_000913/2001762‑2002055
                                                                                                                                                       |                                                                                                                                                  
cacCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGa                                                                                                                                                   <  1:1008032/151‑1 (MQ=255)
cacCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGa                                                                                                                                                   <  1:345286/151‑1 (MQ=255)
    ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                                >  2:1240830/1‑150 (MQ=255)
    ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                                >  1:789166/1‑150 (MQ=255)
     tCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaa                                                                                                                                               >  2:1108698/1‑150 (MQ=255)
        gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                         <  1:554534/151‑1 (MQ=255)
          ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                       <  2:772795/151‑1 (MQ=255)
          ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                       >  1:835369/1‑151 (MQ=255)
          ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcccc                                                                                                                                       >  2:992582/1‑151 (MQ=255)
              taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                   >  2:884192/1‑151 (MQ=255)
               agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                  <  1:1457565/151‑1 (MQ=255)
                 agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                >  1:1386867/1‑151 (MQ=255)
                     aaTTCACGATAAACAGCCCTGCTTTATATGATTTATGGGAATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTTGGGGGTTAAAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGc                                                                                                                              <  2:549481/149‑1 (MQ=255)
                         cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                         >  1:172310/1‑150 (MQ=255)
                           cGAAAAACACCCTTGTTTTATATGATTTACCGGCTTGTTTATCCTTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTTGGGGTTTACAGGAAAATTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGc                                                                                                                      <  1:443359/151‑1 (MQ=255)
                               aaaCAGCCCTGCGTTATAGGAGTTATCGGCTTTTTTTTCCGTTTATGCAGTTTTTTTACTGGTTCTTTTAGTGGTAAAATTGTGGGGTTTTTTGGGTGTGACGGGAAATTTCCAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCa                                                                                                                     <  2:952907/148‑1 (MQ=255)
                               aaaCACCCCTGCGTTATATGATTTATCGGCTTGATTATCCGTTTCTGCGGGGTTTTTAATCGAACGATTAGGGGGTAAAAGGAGGGGTTATTTGGGGGTTAAAGTAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                    <  2:87245/149‑1 (MQ=255)
                                   aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt                                                                                                              >  1:1201443/1‑151 (MQ=255)
                                    gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                             >  1:1011225/1‑151 (MQ=255)
                                      ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                           >  1:821729/1‑151 (MQ=255)
                                        tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                          >  1:300498/1‑150 (MQ=255)
                                           gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                      >  1:666014/1‑151 (MQ=255)
                                            ttATATGAGTTATCGGCATGTTTACCCGTTTCTGCGGGGTTTTTAATGGGACGATTGGTGGGTGAAATGAGGGTTTATTGGGGGGTTAAAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                      <  1:582971/150‑1 (MQ=255)
                                            ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                     >  1:411609/1‑151 (MQ=255)
                                              atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                   >  1:1555619/1‑151 (MQ=255)
                                                     ttATGGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGGTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                               <  1:1564750/148‑1 (MQ=255)
                                                     ttATCGGCTTGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGTTTAGTGGGTGAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                   <  2:156213/144‑1 (MQ=255)
                                                      tATCGGATTGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGTTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                            <  1:134192/150‑1 (MQ=255)
                                                       atgggaaTGATTATCCGTTTCGGCAGGTTTTTTAATCGGAGGATTAGGGGGGAAAAGGGGGGTTTATTGGGGGTTTACGGATAAATTCCAGCAAAAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                               <  2:78930/140‑1 (MQ=255)
                                                            cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                     >  2:57812/1‑151 (MQ=255)
                                                                aTTATCTGTTTCTGCGGGGTTTTAAATCGGACGATTAGGGGGTAAAATGGGGGGTTATTGGGGGGTTAAAGGTAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                  <  1:1239211/150‑1 (MQ=255)
                                                                 ttccccGTTTTTGAGGGTTTTTTAGTCGGTGATTTAGTGGGTCAAATGAGGGTTCAGTGGGGGCTAACGGCTAATTTCCAGCAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                   <  1:509902/144‑1 (MQ=255)
                                                                 ttATCCTTTTCGGCAGGTTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGTTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                    <  1:1130764/147‑1 (MQ=255)
                                                                  tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                               >  1:755183/1‑151 (MQ=255)
                                                                  tATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                               >  2:715812/1‑151 (MQ=255)
                                                                   aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                              >  1:1393449/1‑151 (MQ=255)
                                                                    tCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                             >  1:721190/1‑151 (MQ=255)
                                                                      cGTTTCTGGAGGGTTTTTAAGCGGACGATTTGTGGGTGAATGGGGGGGTTTTTTGGGGTTTAAAGGTAAATACCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                  <  2:635481/144‑1 (MQ=255)
                                                                        tttCTGCAGGTTTTTTAATCGGACGATTAGTGGGTGAAAGGGGGGTTTTTTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                          <  2:647643/150‑1 (MQ=255)
                                                                             gCAGGTTTTTTAATCGGACGATTGGGGGGTGAAAGGAGGGGTTATTGGGGGGTTACAGTTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                    <  2:128214/151‑1 (MQ=255)
                                                                             gCAGGGTTTTAAATCGGACGATTGGGGGGTAAAATGAGGGTTTATTGGGGGGTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                       <  1:534550/148‑1 (MQ=255)
                                                                             gCAGGCATTTTAATCGACCGATTCGAGTGCGATATGTCGAGTTGTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGCGACAGGACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTTAACg                                                                       >  2:319163/1‑148 (MQ=255)
                                                                               aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                   >  1:1240483/1‑150 (MQ=255)
                                                                                ggttttttATTGGGACGATTGGGGGGTAAAATGAGGGTTTATTTGGGTTTTACAGAAAATTTCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                      <  1:1257184/143‑1 (MQ=255)
                                                                                ggttttttAATCGGACGATTAGTGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                       <  1:683654/142‑1 (MQ=255)
                                                                                gggTTTTTATCCGGATAATGAGTGGGGGAAAGGAGGGGTTATTTGGGGTTTACGGGAAAATCCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                        <  1:294418/144‑1 (MQ=255)
                                                                                gggTTTTTAATCGGACGATTAGTGGGTGAAATGCGGGGTTCTTTGGGGGTTACAGGTAAATTACAGGCTGAAAAAAACCCCGCCTTTGGCGGGGAAGCACGTTTCTGACAAATTACGCTTCATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                 >  2:873117/1‑151 (MQ=255)
                                                                                 ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                >  1:1104828/1‑151 (MQ=255)
                                                                                 ggTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                >  1:75801/1‑151 (MQ=255)
                                                                                   tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                              >  2:1416853/1‑151 (MQ=255)
                                                                                   tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                              >  1:1099963/1‑151 (MQ=255)
                                                                                    ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                             >  1:904394/1‑151 (MQ=255)
                                                                                       tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                          >  2:509679/1‑151 (MQ=255)
                                                                                       tAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                          >  1:1538029/1‑151 (MQ=255)
                                                                                        aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                          >  1:312585/1‑150 (MQ=255)
                                                                                           cGGACGTTTATGGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                        <  2:673153/149‑1 (MQ=38)
                                                                                            ggACGATTGGGGGGTAAAATGGGGGTTTTTTGGGGGTTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                     <  1:309434/151‑1 (MQ=37)
                                                                                            ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                      >  1:646219/1‑150 (MQ=255)
                                                                                             gACGATTAGTGGGAGAAACGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAACAACCCCGCCGGTGGCGGGGAAGCACGAAGCTGACAAATTGCGCTTTAAGTTGCCGGCGGCGGCGTACACACCCTATCTGGCCGACaaa                                                     >  2:109716/1‑150 (MQ=255)
                                                                                                  ttggggggaaaaaTGGGGGGTTTTTTGGGGTTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                   <  1:660530/137‑1 (MQ=25)
                                                                                                  ttggggggTAAAAGGAGGGGTTATTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                 <  2:985929/144‑1 (MQ=38)
                                                                                                  ttAGTGGGTGAAAGGGGGGTTTATTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                <  1:1238797/150‑1 (MQ=38)
                                                                                                     gTGGGAGAAATGGGGGGTTTTTTGGGGTTAACGGGAAAATGCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                <  1:178973/147‑1 (MQ=25)
                                                                                                      tGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa                                            >  2:945406/1‑150 (MQ=255)
                                                                                                       gggTGAAAGGGGGGTTTTTTTGGGGTTTACAGTTAAATTCCGGGAAAAAAAAATCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                           <  2:1234583/150‑1 (MQ=37)
                                                                                                       ggaaaaaaTGAGGGGTGAATGGGGGGTTAAATAAAACTCCCAGTAAAAAAAAAACACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGGTGCGGCATAAACGCCTTATCCCGCCTACAAAAATGTGCaaa                                           <  2:1204482/145‑1 (MQ=255)
                                                                                                              aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                    >  2:677910/1‑150 (MQ=255)
                                                                                                              aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                    >  1:801162/1‑150 (MQ=255)
                                                                                                                gggggtTTATTGGGGGTTTAAAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                   <  1:434839/147‑1 (MQ=37)
                                                                                                                gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                 >  1:646871/1‑151 (MQ=255)
                                                                                                                gAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                 >  2:965496/1‑151 (MQ=255)
                                                                                                                gAGGGGTTATTGGGGGTTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaaa                                  <  1:1513624/150‑1 (MQ=37)
                                                                                                                    gggtttttGGGGGTTAACAGGTAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                              <  2:1201427/146‑1 (MQ=25)
                                                                                                                     gTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                             >  1:1372126/1‑150 (MQ=255)
                                                                                                                      ttATTTGGGGGTTAGAGGTAAATTCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                   <  2:96459/143‑1 (MQ=255)
                                                                                                                      ttATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                           >  1:1539435/1‑151 (MQ=255)
                                                                                                                       ttttttGGGGGTAAAGGTAAATTCCGGGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                              <  2:858704/146‑1 (MQ=25)
                                                                                                                       tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                          >  2:612940/1‑151 (MQ=255)
                                                                                                                        aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                         >  2:106080/1‑151 (MQ=255)
                                                                                                                        aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                         >  2:217827/1‑151 (MQ=255)
                                                                                                                         ttggggggTTAAAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                         <  1:472559/147‑1 (MQ=37)
                                                                                                                            gggggTTAAAGATAAATTCAAGCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                          <  1:558641/146‑1 (MQ=37)
                                                                                                                             ggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                    >  1:373187/1‑151 (MQ=255)
                                                                                                                              ggggataacggttAAATTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                       <  1:1108713/137‑1 (MQ=25)
                                                                                                                              gggTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                   <  1:1214048/151‑1 (MQ=255)
                                                                                                                               ggTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                  <  1:367868/151‑1 (MQ=255)
                                                                                                                                     aCAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc              <  1:371481/151‑1 (MQ=255)
                                                                                                                                      cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc             >  2:4494/1‑151 (MQ=255)
                                                                                                                                      acggataaTTCCCAGCAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc              <  1:1304215/145‑1 (MQ=25)
                                                                                                                                           aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt         >  1:198176/1‑150 (MQ=255)
                                                                                                                                           aaaTCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa        <  2:789152/151‑1 (MQ=255)
                                                                                                                                           aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCCACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc             <  1:1088386/146‑1 (MQ=255)
                                                                                                                                            aaTTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa       <  1:205458/151‑1 (MQ=255)
                                                                                                                                            aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa        >  1:941893/1‑150 (MQ=255)
                                                                                                                                                ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc   >  1:839549/1‑151 (MQ=255)
                                                                                                                                                  aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc   >  2:692611/1‑149 (MQ=255)
                                                                                                                                                  aGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca  <  2:1172899/150‑1 (MQ=255)
                                                                                                                                                       gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg               >  2:5929/1‑132 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCA  >  NC_000913/2001762‑2002055

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