Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,111,7470GT30.0% 100.7 / 20.8 40N156K (AAC→AAAyiiQDUF1454 domain‑containing protein YiiQ
Reads supporting (aligned to +/- strand):  ref base G (27/1);  new base T (0/12);  total (27/13)
Fisher's exact test for biased strand distribution p-value = 2.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGG  >  NC_000913/4111608‑4111886
                                                                                                                                            |                                                                                                                                           
ttttttGCCATTTATTCAACCCTTCAAGGGATTCAAATAGCGCAAGCTTAACCGGTCCAACAGCGAAGGTCAGCCCTTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCCCGGTGTAATAACgtttg                                                                                                                                             <  2:234511/139‑1 (MQ=255)
ctttTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACgtttg                                                                                                                                             >  2:83386/1‑140 (MQ=255)
    ttGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGtttt                                                                                                                                         >  1:955048/1‑140 (MQ=255)
        tttttATTAAACCCTACAAGCGATTCAGATAGCGCCAGTTTAATCGGTCCAACAGCGAAGGTCACCCCTTTTTCCCCGTTGCCCGCGAAAACATACCGCAGTCCCCCTTCTTTCTGGGTGTAATAACTTTTTTTTTcccc                                                                                                                                    <  2:619559/138‑1 (MQ=255)
             tttaaacccTTCAAGCGATTCAGATAGCGCCAGTTTAATCGTTTCAACAGCGAAGGCCACCCCTTTTTCGCCGTTGCCCGCAACAACATAACGCAGTGCCCCTTCTGTCTGGGTGTAATAACTTTTTTTTTTCCccgccg                                                                                                                                <  2:545102/133‑1 (MQ=255)
                       tCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGa                                                                                                                      >  1:1171057/1‑140 (MQ=255)
                        cAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGAc                                                                                                                     >  1:696939/1‑140 (MQ=255)
                          agctttacAGATAGCGCCAGTTTAATCGTCACAACAGCGAAGGTCCCCCCTTTTTCGCATTTGTCCGCGACAACATAACGCAGTGCCCCTTCTTCCTGGGTTTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTg                                                                                                                   <  2:904955/133‑1 (MQ=255)
                             aaTTCAGATGCCGCCATTTAAATCGTTAAAACAGCGAGGGTCGCCCCTTTTCCCCCTTTGCCCGCGACAACAAAACGCAGTCCACCTTCTGTCGGGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAg                                                                                                                <  2:655060/139‑1 (MQ=255)
                                    aTAGCGCAAGCTTAATCGGTTAAACAGCGAGGGTCACCCCTTTTTCGCCGTTGTCCGCGAAAACATAACGCAGTCCCCCTTCTGTCTGGGTGTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTc                                                                                                         <  2:1057032/140‑1 (MQ=255)
                                         gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGcc                                                                                                    >  1:655060/1‑139 (MQ=255)
                                         gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                    >  1:152148/1‑140 (MQ=255)
                                         gcCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCt                                                                                                    >  1:142143/1‑140 (MQ=255)
                                              cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                               >  2:160766/1‑140 (MQ=255)
                                              cTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTtgtg                                                                                               >  1:1251495/1‑140 (MQ=255)
                                                     ggTTCAACAGCGAGGGCCACCCCTTTTTCGCGGTTGTCCGCGAAAACATAACGCAGTGCTCCTTCTGTCTGGGTTTAATACCTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGt                                                                                        <  2:152148/140‑1 (MQ=255)
                                                       ttCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCa                                                                                      >  2:1139992/1‑140 (MQ=255)
                                                         cAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCAtt                                                                                    >  1:530264/1‑140 (MQ=255)
                                                               cGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATggg                                                                              >  2:470042/1‑140 (MQ=255)
                                                                 aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            >  1:236617/1‑140 (MQ=255)
                                                                 aaGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            >  2:354033/1‑140 (MQ=255)
                                                                  agggccgCCCCTTTTCCGCTTTTGTCCGCGAAAACATAACGCAGTGCCCCTTCTGTCTGGGTTTAATACTGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            <  2:142143/134‑1 (MQ=255)
                                                                  agggccaccccTTTTCCGCGGTTGCCCGCAACAACATACCGCAGTCCCCCTTCTGTCTGGGTGTAATACCGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGt                                                                            <  2:218305/133‑1 (MQ=255)
                                                                              tttCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                               >  2:165815/1‑140 (MQ=255)
                                                                              tttCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATc                                                               >  1:1363054/1‑140 (MQ=255)
                                                                                   cccGTTGCCCGCGACAACATAACGCAGTGCCCCTTCTGTCTCGGTGAAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCg                                                           <  1:944292/139‑1 (MQ=255)
                                                                                      tttGCCCGCCACAACAAAACGCAGTGCCCCTTCTGTCTGGGTGTAATAACTTTTGTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                        <  1:1282956/139‑1 (MQ=255)
                                                                                      tttGCCCGCAAAAACATACCGCGGTGCCCTTTCTGTCTGGGTTTAATAACGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                        <  2:57377/139‑1 (MQ=255)
                                                                                      gTTGTCCGCAAAAACATAACGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                        <  1:605955/140‑1 (MQ=255)
                                                                                      gTTGGCCGCAAAAAAAAACCGTAGTGCCCTTTCTGTCGCGGTTTATAATGGTTTGTTTTCCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg                                                        <  2:1088376/140‑1 (MQ=255)
                                                                                        tGTCCGCGAAAACATAACGCAGTCCCCCTTTTGTCTGGGTGTAATAATTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc                                                      <  2:1090051/140‑1 (MQ=255)
                                                                                                 caacaTAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTcc                                             >  1:598369/1‑140 (MQ=255)
                                                                                                  aacaTAACGCAGCCCCCTTTCTTTCTGGGTGTAATATGGTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCt                                            <  2:151792/140‑1 (MQ=255)
                                                                                                         cGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgcttt                                     >  1:567695/1‑140 (MQ=255)
                                                                                                              tGCACCTTCTGTCTCGGTGCAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgctt                                >  2:556647/1‑140 (MQ=255)
                                                                                                               gCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAgctttcgcttt                               >  1:646714/1‑140 (MQ=255)
                                                                                                                   tttCTGTCTGCGTTTAATACTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTcgcg                           <  2:956155/139‑1 (MQ=255)
                                                                                                                        gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                      >  2:1126615/1‑140 (MQ=255)
                                                                                                                        gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                      >  2:1075111/1‑140 (MQ=255)
                                                                                                                        gTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCttt                      >  2:766885/1‑140 (MQ=255)
                                                                                                                             ggTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                 >  1:229230/1‑140 (MQ=255)
                                                                                                                             ggTGAAATACCTTTTTTTTTCCTCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGct                 <  2:1171057/140‑1 (MQ=255)
                                                                                                                              gtgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctc                >  1:730191/1‑140 (MQ=255)
                                                                                                                               ttaaataatTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct               <  2:955048/137‑1 (MQ=255)
                                                                                                                               tgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct               >  2:118537/1‑140 (MQ=255)
                                                                                                                                 taataaCGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTgg             >  1:515641/1‑140 (MQ=255)
                                                                                                                                     atttttttttttCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGcccc         <  2:236617/136‑1 (MQ=255)
                                                                                                                                       cGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTg       >  1:1334041/1‑140 (MQ=255)
                                                                                                                                         tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  1:120148/1‑140 (MQ=255)
                                                                                                                                         tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:1374345/1‑140 (MQ=255)
                                                                                                                                         tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:947167/1‑140 (MQ=255)
                                                                                                                                         tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:583443/1‑140 (MQ=255)
                                                                                                                                         tttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGa     >  2:461896/1‑140 (MQ=255)
                                                                                                                                           tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGctct               >  1:273467/1‑128 (MQ=255)
                                                                                                                                           tgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATg   >  2:633634/1‑140 (MQ=255)
                                                                                                                                            gtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATgg  <  1:633634/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CTTTTTGTCATTTATTCAAACCTTCAAGCGATTCAGATAGCGCCAGCTTAATCGGTTCAACAGCGAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCCATGTATTCCTGAGCTTTCGCTTTCGCGGCTTTTTGCTCTGGCCCCTGGATGG  >  NC_000913/4111608‑4111886

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: