Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,517,587 2 bp→TT intergenic (+25/‑60) ortT → / → ydcY orphan toxin OrtT/DUF2526 domain‑containing protein YdcY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,517,5870CT84.2% 48.9 / ‑0.3 19intergenic (+25/‑61)ortT/ydcYorphan toxin OrtT/DUF2526 domain‑containing protein YdcY
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (9/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009131,517,5880CT84.2% 46.5 / 0.4 19intergenic (+26/‑60)ortT/ydcYorphan toxin OrtT/DUF2526 domain‑containing protein YdcY
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (9/7);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAG  >  NC_000913/1517366‑1517828
                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                
ttCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATc                                                                                                                                                                                                                      >  1:438822/1‑251 (MQ=255)
             aCGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTcaca                                                                                                                                                                                                         <  1:95367/251‑1 (MQ=255)
                   gTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGc                                                                                                                                                                                                   >  1:413958/1‑251 (MQ=255)
                                tATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATAGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAgg                                                                                                                                                                                       >  1:213491/1‑249 (MQ=255)
                                 aTCAACACATGCTTGTTTTTTATGCGGTTATAGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTa                                                                                                                                                                                      <  1:20340/250‑1 (MQ=255)
                                                           gTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATcgcgc                                                                                                                                                           <  2:271808/251‑1 (MQ=255)
                                                                     caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAgg                                                                                                                                                                                              <  1:419888/206‑1 (MQ=255)
                                                                     caATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAgg                                                                                                                                                                                              >  2:419888/1‑206 (MQ=255)
                                                                           cATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGAt                                                                                                                                           <  1:97362/251‑1 (MQ=255)
                                                                                 ttATCACCTGGATTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGa                                                                                                                                        >  2:47375/1‑248 (MQ=255)
                                                                                           gTTTCTTTCTCACGCTAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCAtt                                                                                                                           >  1:46503/1‑251 (MQ=255)
                                                                                                                                           ggtggGGGCAACATGGCCGATGAGTTTCCCGAAGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGtt                                                                            >  1:434646/1‑250 (MQ=255)
                                                                                                                                                        tGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGa                                                                   >  2:468342/1‑246 (MQ=255)
                                                                                                                                                                         ggTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAgcagca                                              >  1:48727/1‑250 (MQ=255)
                                                                                                                                                                          gTGNCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTTGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACTCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAAc                                            <  2:438822/251‑1 (MQ=255)
                                                                                                                                                                            ggCGCTGTTGTTTTCACTGTTTTAACTCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAACGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCaaagt                                          >  2:405685/1‑248 (MQ=255)
                                                                                                                                                                                                    aaTCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGTCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACTCCCAGCAGCAACGTCTGCGCACAGCAATCGCCCATCAc                  <  2:213491/251‑1 (MQ=255)
                                                                                                                                                                                                     aTCTTACTCACTGGCTTGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACg                 >  1:154316/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGACGGCATTACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAaag  <  1:192421/250‑1 (MQ=255)
                                                                                                                                                                                                                             ||                                                                                                                                                                                                                                                
TTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCACGATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCTGAGCGATGACGCGGAGTTAAGTCGGGAAGATCGCTACACCCAGCAGCAACGTCTGCGCACAGCGATCGCCCATCACGGTCGCAAGCATAAAG  >  NC_000913/1517366‑1517828

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