Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,4960CT20.0% 97.9 / 14.6 40intergenic (+236/+184)yjjV/yjjWputative DNase YjjV/putative glycyl‑radical enzyme activating enzyme YjjW
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base T (5/3);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACG  >  NC_000913/4614359‑4614620
                                                                                                                                             |                                                                                                                            
tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCgg                                                                                                                            >  1:1074027/1‑140 (MQ=255)
tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCgg                                                                                                                            <  2:1074027/140‑1 (MQ=255)
tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCgg                                                                                                                            <  1:799321/140‑1 (MQ=255)
        aTAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAAt                                                                                                                    <  1:565339/140‑1 (MQ=255)
                          aTCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGt                                                                                                  <  1:1042712/140‑1 (MQ=255)
                              ggCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcgcg                                                                                              <  2:1045299/140‑1 (MQ=255)
                                         cacaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCagg                                                                                    >  1:254812/1‑138 (MQ=11)
                                         cacaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCagg                                                                                    >  1:468428/1‑138 (MQ=11)
                                         cacaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGc                                                                                   >  2:972828/1‑140 (MQ=255)
                                          acaATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGcc                                                                                  >  2:706013/1‑140 (MQ=255)
                                            agtGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTa                                                                                <  2:439545/138‑1 (MQ=255)
                                                 cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGt                                                                           <  1:706013/140‑1 (MQ=255)
                                                  tGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTc                                                                                            <  2:1015856/122‑1 (MQ=255)
                                                  tGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTc                                                                                            >  1:1015856/1‑122 (MQ=255)
                                                    aTGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTtg                                                                        >  1:509700/1‑140 (MQ=255)
                                                      gCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgtg                                                                      >  2:173016/1‑140 (MQ=255)
                                                        gACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGcc                                                                    >  1:663530/1‑140 (MQ=255)
                                                         aCGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGcct                                                                   >  1:439545/1‑139 (MQ=255)
                                                              tGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCg                                                              <  2:509700/140‑1 (MQ=255)
                                                                    gTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGcc                                                        >  2:500288/1‑140 (MQ=255)
                                                                               ccTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt                                             <  2:125592/140‑1 (MQ=255)
                                                                                cTACAAGTCTGTGCCGACCCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                            <  1:486086/140‑1 (MQ=255)
                                                                                cTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                            >  2:231660/1‑140 (MQ=255)
                                                                                  aCAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcg                                          >  1:126270/1‑140 (MQ=255)
                                                                                   cAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                                                 <  2:366009/132‑1 (MQ=255)
                                                                                   cAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                                                 >  1:366009/1‑132 (MQ=255)
                                                                                    aaGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcacg                                        >  2:1005045/1‑140 (MQ=255)
                                                                                      gtttgtgCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc                                      <  1:1102142/137‑1 (MQ=2)
                                                                                       tttgtgCCGAACCGTAGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcgcg                                        <  1:697028/138‑1 (MQ=14)
                                                                                           tgCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATc                                 >  1:109180/1‑140 (MQ=255)
                                                                                             ccGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCg                               >  2:1155587/1‑140 (MQ=255)
                                                                                             ccGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCg                               >  2:134203/1‑140 (MQ=255)
                                                                                              cGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCgg                              >  2:767432/1‑140 (MQ=255)
                                                                                                aaCCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCa                            >  2:1102142/1‑140 (MQ=17)
                                                                                                aaCCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCa                            >  2:94811/1‑140 (MQ=255)
                                                                                                 aCCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAg                           <  1:1155587/140‑1 (MQ=255)
                                                                                                     t‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCAACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTgg                       <  2:730967/140‑1 (MQ=16)
                                                                                                     t‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTgg                       >  1:820369/1‑140 (MQ=255)
                                                                                                                          cACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCg     <  1:94811/140‑1 (MQ=255)
                                                                                                                             gccgcATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACg  <  1:500288/140‑1 (MQ=255)
                                                                                                                                             |                                                                                                                            
TAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGT‑‑‑AGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACG  >  NC_000913/4614359‑4614620

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: