Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC40.9% 175.1 / 12.9 89intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (46/6);  new base C (0/36);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 8.49e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001771‑2002064
                                                                                                                                                          |                                                                                                                                                  
cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAAccc                                                                                                                                                 <  1:979035/151‑1 (MQ=255)
 ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAAcccc                                                                                                                                                >  1:102878/1‑151 (MQ=255)
  gTATAGAGTGAATTCACGATAAACAGCCCTGTGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCccg                                                                                                                                               <  2:1254729/151‑1 (MQ=255)
  gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCccg                                                                                                                                               <  2:44640/151‑1 (MQ=255)
   tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCccgc                                                                                                                                              >  2:328153/1‑151 (MQ=255)
     taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCccgccg                                                                                                                                            >  1:502404/1‑151 (MQ=255)
     taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCccgccg                                                                                                                                            >  2:507039/1‑151 (MQ=255)
           gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCg                                                                                                                                      >  2:1285144/1‑151 (MQ=255)
            aaTTCACGATAAACAGCCCTGCTTTATAGGAGTTATCGGCATGTTTACCCTTTTCGGCGGGGTTTTA‑AACCGGGAGATTAGGGGGTGAAAGGAGGGT‑TTATTGGGGGG‑TT‑ACAGG‑TAAATCCC‑G‑GGAAAAAAAAACCCCCGCCGGTGGCg                                                                                                                                      <  2:1188351/150‑1 (MQ=255)
                  cGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGGTTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGTAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGCAGc                                                                                                                               >  1:310561/1‑151 (MQ=255)
                     tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                            >  1:440275/1‑151 (MQ=255)
                     tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                            >  2:134519/1‑151 (MQ=255)
                        aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                          >  2:450800/1‑150 (MQ=255)
                            cccTGTTTTATATGATTTATCGGCATGATTTTCCGTTTCTGCAGGTTTTTT‑AATGGGACGATTGGTGGGTAAAATGGGGGG‑TTATTGGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGa                                                                                                                     <  1:870069/151‑1 (MQ=255)
                            cccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                      >  2:688959/1‑150 (MQ=255)
                               tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                   >  1:612753/1‑150 (MQ=255)
                                gtgTTATAGGAGTTATGGGAATGTTTATCGGTTTCTGCGGGGTTTTA‑AATCGGACAATGGGGGGGTAAAATGGGGGG‑TTATTGGGGGG‑TT‑ACAGG‑TAAATCCA‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                   <  2:1104516/147‑1 (MQ=255)
                                 cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                                 >  2:447896/1‑150 (MQ=255)
                                  gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                               >  1:215459/1‑151 (MQ=255)
                                   ttATATGGGTTATCGGCATTTTTACCCGTTTCTGCGGGTTTTTA‑AAGCGGACGATTAGGGGGTGAAATGAGGGG‑TTTTTTGGGGG‑TT‑ACAGG‑AAAATCCC‑G‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                <  1:169850/149‑1 (MQ=255)
                                   ttATATGATTTATCGGCATGATTATCCGTTTCTGCGGGTTTTTT‑AAGCGGACGATTGGTGGGTAAAAGGGGGGT‑TTTTTGGGGGG‑TT‑ACAGG‑AAAATTCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                  <  1:757105/147‑1 (MQ=255)
                                     atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                            >  1:228892/1‑151 (MQ=255)
                                            ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg                                                                                                     >  2:989194/1‑151 (MQ=255)
                                              aTCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    >  2:831802/1‑150 (MQ=255)
                                              aTCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGCGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACGGG‑TAAAATCC‑A‑GGCAGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                    >  2:491507/1‑150 (MQ=255)
                                                 ggCATGATTATCCGTTTCTGCGGGTTTTTT‑AATCGGACGTTTGGGGGGTGAAAGGGGGGG‑TTTTTGGGGGG‑TT‑AAAGG‑TAAATCCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                  <  2:141528/149‑1 (MQ=255)
                                                 ggCATGATTATCCGTTTCTGCGGGGTTTTA‑AATCGGACGTTTAGGGGGTGAAATGGGGGT‑TTATTGGGGGT‑TT‑ACAGG‑AAAATCCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                  <  1:127594/149‑1 (MQ=255)
                                                  gCATGATTATCCGTTTCTGCAGGGATTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                               >  2:876088/1‑151 (MQ=255)
                                                       aTTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                          >  1:470422/1‑151 (MQ=255)
                                                       aTTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                          >  2:837343/1‑151 (MQ=255)
                                                        ttATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                         >  2:411975/1‑151 (MQ=255)
                                                         tttcgtttgcTGCAGGTTTTTA‑AATGGAACGATTGGGGGGAAAAATGAGGGG‑TT‑TTTGGGGG‑TTAAAAGG‑TAATTCCA‑A‑GGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                         <  2:1270499/141‑1 (MQ=255)
                                                          aTCCTTTTCTGCGGGGTTTTT‑AATCGGACGATTAGTGGGTAAAATGAGGGT‑TTTTTTGGGGT‑TT‑ACAGG‑AAAATTCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                         <  1:735571/149‑1 (MQ=255)
                                                          aTCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                       >  1:805263/1‑151 (MQ=255)
                                                               tttCTGCAGGTTTTTT‑AATCGGACGATTAGGGGGTGAAATGAGGGG‑TTTTTTGGGGT‑TT‑ACAGT‑TAAATCCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                     <  2:740407/148‑1 (MQ=255)
                                                               tttCTGCAGGGTTTTT‑AACCGGACTATTAGTGGGTGAAAGGAGGGG‑TTATTGGGGGG‑TT‑CCAGG‑TAAATCCC‑A‑GGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                   <  1:556902/150‑1 (MQ=255)
                                                                ttCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                 >  1:1031895/1‑151 (MQ=255)
                                                                ttCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                 >  2:309519/1‑151 (MQ=255)
                                                                   tGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               >  1:843103/1‑150 (MQ=255)
                                                                     gcgggTTTTTAAACCGGACGATTAGGGGGT‑AAAAGGGGGG‑TTTTTTGGGGGTTT‑ACAGT‑TAAATTCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                               <  1:757525/147‑1 (MQ=255)
                                                                        gtttttt‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                         >  1:232555/3‑151 (MQ=255)
                                                                           tttt‑AGTCGGAAG‑TTGGGGGGGGAAAGGGGGGGTTTAGGTGGGGG‑TA‑ACGGG‑TAAATTCA‑A‑GACAAAAAAAACCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                        <  2:1093375/149‑1 (MQ=19)
                                                                           tttt‑AATCGGGCGTTAGGTGGGTGAAAGGGGGTT‑TTATTGGGGGT‑TC‑AAAGG‑TCAATACC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcgg                                                                       <  2:1343987/150‑3 (MQ=255)
                                                                             tt‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                    >  1:1191731/1‑151 (MQ=255)
                                                                              t‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                    >  2:29394/1‑150 (MQ=255)
                                                                                acggggACGATTAGTGGGTAAAATGGGGGG‑TTATTGGGGGT‑TA‑ACAGG‑TAAATTCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                   <  1:643534/146‑1 (MQ=38)
                                                                                aaTCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                   >  1:1160895/1‑150 (MQ=255)
                                                                                  tCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                 >  1:278537/1‑150 (MQ=255)
                                                                                    ggaggaTTGGTGGGTAAAATG‑GGGG‑TTTTTTGGGGGTTT‑ACAGG‑TAAATCCA‑A‑GCCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                              <  1:749505/151‑1 (MQ=25)
                                                                                    ggaggaTTGGGGGGT‑AAAAGGGGGG‑TTTTTTGGGGGTTT‑ACAGG‑TAAATCCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                              <  1:507047/151‑1 (MQ=25)
                                                                                    ggATGATTGGTGGGTGAAAGGAGGGG‑TTATTTGGGGT‑TT‑AAAGG‑TAAATCCC‑G‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                <  2:843021/149‑1 (MQ=37)
                                                                                          ttggggggaaaaaTGGGGGG‑TTATTGGGGGT‑TT‑ACAGG‑TAAATTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                            <  1:737904/137‑1 (MQ=37)
                                                                                          ttAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGA‑TAAATCCC‑A‑GAAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt                                                         <  1:497161/150‑1 (MQ=255)
                                                                                          ttAGGGGTTAAAATGGGGGT‑TTATTTGGGGG‑TT‑A‑AAG‑TAAATTCCAA‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          <  2:1087075/149‑1 (MQ=25)
                                                                                          ttAGGGGGTGAAATG‑GGGG‑TTTTTTGGGGGTTT‑ACAGG‑TAAATTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          <  1:899046/149‑1 (MQ=37)
                                                                                            agggggTAAAATG‑GGGG‑TTTATTGGGGG‑TT‑ACAGGATAAATCCC‑A‑GCAAAAAAAACTCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                       <  1:215052/147‑1 (MQ=25)
                                                                                             gggggTAAAAGGGGGGT‑TTATTGGGGGG‑TT‑ACAGG‑AAAATCCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                        <  1:1281329/146‑1 (MQ=37)
                                                                                             gTGGGTGAAGTGGGGGT‑TTTTTTGGGGT‑TT‑ATGGT‑TAAATTCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                          <  1:233579/146‑1 (MQ=37)
                                                                                               gggTGAAAGGAGGGG‑TTATTGGGGGG‑TT‑ACAGT‑AAAATCCC‑A‑GGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                   <  1:505280/151‑1 (MQ=38)
                                                                                               gggTGAAA‑GGGGGG‑TTTTTTGGGGG‑TTAAAAGA‑AAAATCCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                       <  2:1038008/147‑1 (MQ=25)
                                                                                                    aaaggggggt‑TTATTGGGGGG‑TT‑ACAGT‑AAAATTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                                  <  2:1145608/141‑1 (MQ=37)
                                                                                                    aaaTGAGGGG‑TTTTTGGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                              <  2:187806/151‑1 (MQ=255)
                                                                                                      aTGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa                                             >  2:1076181/1‑150 (MQ=255)
                                                                                                       tGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                            >  2:677325/1‑150 (MQ=255)
                                                                                                       ggggggt‑TTTTTTGGGGT‑TT‑CCAGG‑TAATTCCC‑A‑GAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                            <  1:282352/147‑1 (MQ=25)
                                                                                                        gggggt‑TTATTGGGGGG‑TT‑ACAGG‑AAAATTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                            <  1:271827/147‑1 (MQ=37)
                                                                                                          gggg‑TTATTGGGGGT‑TT‑ACAGG‑AAAATCCC‑G‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                        <  1:1100719/151‑1 (MQ=25)
                                                                                                           ggg‑TTATTGGGGGG‑TT‑ACAGG‑TAAATCCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                       <  1:1104051/151‑1 (MQ=255)
                                                                                                               ttttttGGGGT‑TT‑ACGGT‑TAATTTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                      <  2:660230/146‑1 (MQ=25)
                                                                                                                 aTTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                  >  1:1307862/1‑151 (MQ=255)
                                                                                                                   tgggggg‑TT‑AAAGT‑AAAATCCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                  <  1:55576/147‑1 (MQ=25)
                                                                                                                    tGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                >  1:1076015/1‑150 (MQ=255)
                                                                                                                    tGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                >  1:568289/1‑150 (MQ=255)
                                                                                                                    tGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                               >  2:511576/1‑151 (MQ=255)
                                                                                                                      gggg‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                             >  1:1302097/1‑151 (MQ=255)
                                                                                                                       ggg‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                            >  2:313964/1‑151 (MQ=255)
                                                                                                                        gg‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           >  1:903697/1‑151 (MQ=255)
                                                                                                                              aaaGG‑TAAATCCC‑G‑GTAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                       <  2:1031874/149‑1 (MQ=25)
                                                                                                                              aCAGT‑TAAATTCC‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                         <  1:825193/149‑1 (MQ=255)
                                                                                                                               cAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                      >  2:664518/1‑151 (MQ=255)
                                                                                                                                aGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                     >  2:738138/1‑151 (MQ=255)
                                                                                                                                 gg‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                    >  2:1137877/1‑151 (MQ=255)
                                                                                                                                  g‑TAAATTCC‑AGGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                             <  2:1007522/142‑1 (MQ=255)
                                                                                                                                     aaaTTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                  >  1:104399/1‑150 (MQ=255)
                                                                                                                                     aaaTCCG‑A‑GGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                           <  2:374425/141‑1 (MQ=255)
                                                                                                                                     aaaTCCC‑A‑GGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt                     <  1:677606/147‑1 (MQ=255)
                                                                                                                                       aTTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                 >  1:1329521/1‑149 (MQ=255)
                                                                                                                                             a‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc            >  1:1140032/1‑149 (MQ=255)
                                                                                                                                               ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag          <  1:1156586/145‑1 (MQ=255)
                                                                                                                                               ggaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca           <  1:1220139/144‑1 (MQ=255)
                                                                                                                                                  aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc       >  1:227536/1‑150 (MQ=255)
                                                                                                                                                    aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc    >  1:481166/1‑151 (MQ=255)
                                                                                                                                                     aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa   <  1:850215/151‑1 (MQ=255)
                                                                                                                                                      aaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAg  <  2:146682/151‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                                  
CGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGG‑TT‑ACAGG‑TAAATTCC‑A‑GGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAG  >  NC_000913/2001771‑2002064

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: ‑