Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,432,105 C→T L120L (CTG→CTA smg ← DUF494 domain‑containing protein Smg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,1050CT93.1% 77.8 / ‑0.5 29L120L (CTG→CTAsmgDUF494 domain‑containing protein Smg
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (14/13);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

AATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCCAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCATAGCAGGGCATTGTAGATATC  >  NC_000913/3431891‑3432350
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
aaTAACCAGTTCAGCCCCGCACTTTGAGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCCGCATACCTTCTTTCACTTCCAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCTGATCATCTCGTTCGAACTCTCCGTTATCCAGCTCGAGCACTCGc                                                                                                                                                                                                                   >  1:207415/1‑251 (MQ=255)
                ccGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCAc                                                                                                                                                                                                   >  2:334001/1‑251 (MQ=255)
                                   cATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCa                                                                                                                                                                                 <  1:115594/250‑1 (MQ=255)
                                     ttattaCGCACCGTGAACAGTGCTGATTTAGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCt                                                                                                                                                                              >  2:205292/1‑251 (MQ=255)
                                                    aaCAGTGCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAgg                                                                                                                                                                >  1:327595/1‑250 (MQ=255)
                                                          gCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCa                                                                                                                                                         <  1:214704/251‑1 (MQ=255)
                                                             gATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAg                                                                                                                                                        <  1:240621/249‑1 (MQ=255)
                                                             gATTTCGCCATAACAACTTATGCGGAATTAAAGAAATGATTAATGCAGCATAGCTTCATTGACTTCAAGGAGTAATTCTTCCATTTGCTGGGACGCATTTTCGCAGGCCGGAGGATTGAACACCACCGTCAGAATCACCCCGTTCAGATCATCTAATTCGAACTCTGCGTTATCGAGCGGGAGCACTCGCTCTATCACCATTTCACGAGTCTCGAGGTTGTGCACCGGAATcgggtcagggataatcata                                                                                                                                                       >  1:325244/1‑232 (MQ=255)
                                                              aTTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGATCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGaaa                                                                                                                                                     <  1:84767/251‑1 (MQ=255)
                                                                                            aGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTcac                                                                                                                        >  1:305256/1‑250 (MQ=255)
                                                                                              aaaTGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACAct                                                                                                                     >  1:40308/1‑251 (MQ=255)
                                                                                                  gATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGTAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTc                                                                                                                 <  1:122262/251‑1 (MQ=255)
                                                                                                     tAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCTCCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCg                                                                                                               <  1:271617/250‑1 (MQ=255)
                                                                                                       aTGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGgtg                                                                                                            >  1:183788/1‑251 (MQ=255)
                                                                                                          cagcaTACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGtat                                                                                                         <  1:205292/251‑1 (MQ=255)
                                                                                                                                aatagaaattttctcATTTGCTGGAACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCCCCCATTTCAGAACATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACGCGCTCTATCACAATTGCACGAGTTTCAAGGTTGTGCACCTGAATCTGCTCAAGGAAAAGCAGAATCCCACGGCAGCTGGCATCCAGTCGTTCACTCTCTTCCGGTGTTTAAATACGCATGGGGAGATGATc                                                                                   <  2:200153/237‑1 (MQ=255)
                                                                                                                                  gagTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGTCCGGAATATTGAACAACACCATAAGGATCACACATTTCAGATCATACAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCAACGGCAGCTGGCATCAAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCg                                                                                 >  2:153620/1‑251 (MQ=255)
                                                                                                                                        ttcttcCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGgag                                                                                           >  1:118649/1‑235 (MQ=255)
                                                                                                                                        gtcttcCAGTTGCTGGTGCGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCAGCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTAGCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGgag                                                                                           <  2:118649/234‑1 (MQ=255)
                                                                                                                                            tcCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATGTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAg                                                                          <  2:40308/248‑1 (MQ=255)
                                                                                                                                                             agtgtctcAGCCCGGAATATTGATCAACACCATCAGGAGCACCCATGTCAGAGCATCTATTTCGAACTCGGCGTTAACCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGGTGAGCAACTGAATCTGCTCAAGGAAAAGCAGAAACCCACTGCATCTGGCATACAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCTAGTGGCATCGGTTCTGCAAAcc                                                      <  2:225801/244‑1 (MQ=255)
                                                                                                                                                                                ttGAACAACACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGATTCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCGTCCTGCTAATCAGCa                                     >  1:150900/1‑249 (MQ=255)
                                                                                                                                                                                   acaatcttgaacAGGTTCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATACAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCACGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCTGTCGTTCACACTCTTCAGGAGTATAAATAAGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGGTCTGCCAACCATTCCTGAAAATCAGCAAGttt                                >  2:188519/12‑251 (MQ=255)
                                                                                                                                                                                   aacaacACCATCAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGttt                                >  1:239731/1‑251 (MQ=255)
                                                                                                                                                                                              cAGGATCACCCATTTCAGATCATCTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCAt                      >  1:240777/1‑250 (MQ=255)
                                                                                                                                                                                                   ccccccATTTCAGATCATCTAGTTCGACCCCTGCGTTATCCAGCGCGAGCCCTCGCTTTATCCCCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCATAGCAgg                <  2:327595/248‑1 (MQ=255)
                                                                                                                                                                                                      cccATTTCAGATCATCCAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGTGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCATAGCAGGGCa             <  1:58067/251‑1 (MQ=255)
                                                                                                                                                                                                      cccATTTCAGAACATCTAGTTCGAACGCTTCGTTATCCAGCGCGAGCACTCGCTCTAGCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTGCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCATAGCAGGGCa             <  1:331040/251‑1 (MQ=255)
                                                                                                                                                                                                                 tcatcTAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCGCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGACAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCATAGCAGGGCATTGTAGATAtc  >  1:181367/1‑251 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
AATAACCAGTTCAGCCCCGCACTTTGGGCAGGACTCATTATTACGCACCGTGAACAGTGCTGATTTCGCCATAACAACTTATGCTGAATTAAAGAAATGATTAATGCAGCATACCTTCATTCACTTCAAAGAGTAATTCTTCCATTTGCTGGTACGCATTTTCGCAGCCCGGAATATTGAACAACACCATCAGGATCACCCATTTCAGATCATCCAGTTCGAACTCTGCGTTATCCAGCGCGAGCACTCGCTCTATCACCATTTCACGAGTTTCAAGGTTGAGCACCTGAATCTGCTCAAGGAAAAGCAGAAACCCACGGCAGCTGGCATCCAGTCGTTCACACTCTTCCGGTGTATAAATACGCATGGAGAGAGGATCCGAGGCCAGTTGCATCGGTTCTGCCAACCCTTCCTGATAATCAGCAAGTTTTTCCAGCCATAGCAGGGCATTGTAGATATC  >  NC_000913/3431891‑3432350

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