Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,535,189 C→T intergenic (‑1136/‑1408) topAI ← / → ythB toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,535,1890CT93.1% 75.8 / ‑4.3 30intergenic (‑1136/‑1408)topAI/ythBtoxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/17);  minor base G (1/1);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

ACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTG  >  NC_000913/4534948‑4535414
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
aCTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAACAAGGACGATTTTCCCCGCAATGGGCTATTTACATATTGATGGTAGAGGCGTGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGACCGCCATTCCGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCAGTTTTCGGCGTAGGTTAATg                                                                                                                                                                                                                          >  2:379544/1‑251 (MQ=255)
aCTCACTTTGCTATCACGGTAAGATGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAAAAATAGCAGCATGGCTGGGGGCACATACGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTCGGTTAATg                                                                                                                                                                                                                          >  2:323652/1‑251 (MQ=255)
aCTCACGTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATg                                                                                                                                                                                                                          <  1:219358/251‑1 (MQ=255)
                      gTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTaa                                                                                                                                                                                                    <  1:284278/251‑1 (MQ=255)
                          cAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTg                                                                                                                                                                                                >  1:367343/1‑251 (MQ=255)
                                 aaaGAGGGATGGAACAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCCCATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTCCAACGTTGTTGGTGTTGTGACGCCATTTTCGGCGGAGGTTAATGCTCTCAAAATGTCATTGCGTACGCTGGAAAttt                                                                                                                                                                                         >  2:245063/1‑250 (MQ=255)
                                               aaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAg                                                                                                                                                                           <  1:323652/251‑1 (MQ=255)
                                               aaaGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGACAAGACg                                                                                                                                                                             >  1:293343/1‑249 (MQ=255)
                                                  ggACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGc                                                                                                                                                                          >  1:132595/1‑249 (MQ=255)
                                                       aTTTTCCCCGCAATGGGATATTTACATATTTATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCCTAGCCCTCTTGAAGCTGAAACAATAACCGCATGGCTGGTGGCACATAATGCTGATATAGAACGCCATTACGGGGAACCACTGACCAAAGTAGTTGGAGTTGTGAAGCCATTTTCGGCGAAGGTTACTGCTTACAAAATGACATTGCGTGAGCTGGCAATTAATgggacagcagagaagggtgt                                                                                                                                                                   >  2:260642/1‑233 (MQ=255)
                                                         tttCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCTCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGTAACCACTGTACAAAGTTGTTGGGGTTGTGACGCCTTTGGCGGCGTAGGTTAATGCTATCAAAATGGCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTa                                                                                                                                                                  <  1:218769/250‑1 (MQ=255)
                                                         tttCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTa                                                                                                                                                                  <  1:245063/250‑1 (MQ=255)
                                                                         tatTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACtgt                                                                                                                                                 >  1:163646/1‑251 (MQ=255)
                                                                                atatTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTgggg                                                                                                                                                      <  1:170189/239‑1 (MQ=255)
                                                                                atatTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAGCCCTCTTGAAGCTGAAACAATAGCAGAATGGCTGGTGGCACATCAGGATGCTATAGCGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTAATTGCGTAAGCTGGAAATTAATGGCAAAGGCGAGCAGGGTTTACTGCCAGTGGGGagctcgggag                                                                                                                                            >  2:170189/1‑240 (MQ=255)
                                                                                                                         cATAGCCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGt                                                                                                 <  1:46574/251‑1 (MQ=255)
                                                                                                                              cccTATTGAAGCTGAAACAATATACGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAACGGGCGATTGTTCTCTTTTCTCCTGTCTACt                                                                                            <  2:293343/251‑1 (MQ=255)
                                                                                                                               ccTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTTGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGATGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCATTGCAGGGGGCGGAAAGGGCGATGGTTCTCTTTTCTTCTGTCTACTc                                                                                           >  2:203444/1‑251 (MQ=255)
                                                                                                                                  cgtGAAGTTGAAACAATAGCAGCATGGCTGGTGTCACATAAGGTTGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGGAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGTCAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGTGCTACTGTTCTCTTTTCTCCTGTTTACTCa                                                                                          <  2:343899/247‑1 (MQ=255)
                                                                                                                                     gAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACa                                                                                     <  1:379544/251‑1 (MQ=255)
                                                                                                                                      aaGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGAGCCATTACGGGGTACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAATTGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTTGGGACTGTGAACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTCTCTACTCAAAACAt                                                                                    <  1:49016/251‑1 (MQ=255)
                                                                                                                                                   agcagcATGGCTGGTGGCACATAAGGATGATATAGAGAGCCCTTACGCGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGCAGATGGCAGGt                                                                       >  1:323802/1‑251 (MQ=255)
                                                                                                                                                         aTGGCTTGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTGTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTGCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGTTTTCTTg                                                                 <  2:163646/251‑1 (MQ=255)
                                                                                                                                                           ggcggggggtacctAAGGATGATAGTGAGCGCCTTTCTGTGGAGCCACTGTAAAAAGTTGTTGGGGTTGGGTCGCTAGTTTCGGTTTATGTTTATGCGATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGATTAGGGTTTTCTGACAGTGGGGACTGTGCACTCCTTTCAGGGTGCGGTAAGTGCGTTTGTTCTCTTTTTTCCTGTCTACTCAAAATATGAAGAGGGCAGGtttagtgat                                                               <  2:323802/238‑7 (MQ=255)
                                                                                                                                                                  ggCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAAc                                                         <  1:112373/250‑1 (MQ=255)
                                                                                                                                                                     acaTAAGGTTGATATAGAGCGCCATTACGGGGAACCACTGGACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGGTTACTGACAGTGGGGACTGTTCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCa                                                     <  1:246352/251‑1 (MQ=255)
                                                                                                                                                                        tAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGa                                                                                                                    >  1:53609/1‑185 (MQ=255)
                                                                                                                                                                        tAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGa                                                                                                                    <  2:53609/185‑1 (MQ=255)
                                                                                                                                                                                                     aaCCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTAGCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGTGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTAGCCCCAATGTTGCTGTCTCACGCGCTaaa                     <  1:142241/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggTGTTGTGACGCCATTTTCGGCGTAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTGTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACGATCCTCAATGGTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTg  <  2:175961/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
ACTCACTTTGCTATCACGGTAAGTTGCAGCCTAAAAGAGGGATGGAAAAAGGACGATTTTCCCCGCAATGGGATATTTACATATTGATGGTAGAGGCATGAAGCCAAATGGAGGGAGCCGTCATAACCCTCTTGAAGCTGAAACAATAGCAGCATGGCTGGTGGCACATAAGGATGATATAGAGCGCCATTACGGGGAACCACTGTACAAAGTTGTTGGTGTTGTGACGCCATTTTCGGCGCAGGTTAATGCTATCAAAATGTCATTGCGTAAGCTGGAAATTAATGGCAAAGACGAGCAGGGTTTACTGACAGTGGGGACTGTGCACTCCTTGCAGGGGGCGGAAAGGGCGATTGTTCTCTTTTCTCCTGTCTACTCAAAACATGAAGATGGCAGGTTTCTTGATAGCAACAGCACTATCCTCAATGTTGCTGTCTCACGCGCTAAAGATAGTTTTCTGGTATTTG  >  NC_000913/4534948‑4535414

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