Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,102,557 C→A 17.5% P98P (CCG→CCT trmI ← tRNA m(7)G46 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,102,5570CA17.5% 77.9 / 6.0 34P98P (CCG→CCTtrmItRNA m(7)G46 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (13/15);  new base A (3/3);  total (16/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 = 9.87e-01

TGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCT  >  NC_000913/3102411‑3102703
                                                                                                                                                  |                                                                                                                                                  
tcccACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGAAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCgg                                                                                                                                                  <  1:528862/147‑1 (MQ=255)
             ggAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAACCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCa                                                                                                                                     <  2:377806/149‑1 (MQ=255)
               aaaaaGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAAt                                                                                                                                   >  2:672267/1‑149 (MQ=255)
                aaaaGAGCTGCACCATGCGCAATGAATTGTCAGGTATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATg                                                                                                                                  >  1:378142/1‑149 (MQ=255)
                                  gcAACGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCa                                                                                                                <  2:704604/149‑1 (MQ=255)
                                    aaTGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAgg                                                                                                              >  2:355372/1‑149 (MQ=255)
                                         aTTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGATTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGAt                                                                                                         <  1:360538/149‑1 (MQ=255)
                                             tCAGGAATCATTTTATGCTGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCttt                                                                                                     <  2:203448/149‑1 (MQ=255)
                                                 gAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCCGAAGCCAGGCACGCACCAACGCCAGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTCGcc                                                                                                 >  1:475744/1‑149 (MQ=255)
                                                 gAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCAGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGcc                                                                                                 >  1:475745/1‑149 (MQ=255)
                                                 gAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCAGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGcc                                                                                                 >  1:475854/1‑149 (MQ=255)
                                                       ttttATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAACGCACTTCAATGCCGAGGAAGTCCTGCCCAGGCCGCTCTTCAGCCATTGcc                                                                                           >  2:357993/1‑149 (MQ=255)
                                                                gcaCTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGAc                                                                                  >  1:207822/1‑149 (MQ=255)
                                                                  acctCACCCGCATCGTGACACATCCCGCGCAGGTTGCTTGAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCTCTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGc                                                                                <  2:177922/146‑1 (MQ=255)
                                                                                gTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAAc                                                                  >  1:123562/1‑149 (MQ=255)
                                                                                 ggACACATCACGCGCAGGTTGCTTAACCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGACGACGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAAcc                                                                 <  2:36461/148‑1 (MQ=255)
                                                                                          aCGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAg                                                        >  1:63050/1‑149 (MQ=255)
                                                                                          aCGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCCCAAg                                                        >  1:63007/1‑149 (MQ=255)
                                                                                           cgcgCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGc                                                       <  1:348109/149‑1 (MQ=255)
                                                                                             cgcAGGTTGCTTAAACCTTCGTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGt                                                     <  1:672267/149‑1 (MQ=255)
                                                                                                aGGTTGCTTAAACCTGCTTCATGCGCAGTAGCCAGGCACGCCCCAACGCCCGGTGACTGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                                                  <  1:531308/149‑1 (MQ=255)
                                                                                                         aaaCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGGTCTTTGGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAACCTCAAGCGTCACCGGCCCTTc                                         >  2:268910/1‑149 (MQ=255)
                                                                                                          aaCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCa                                        <  1:644432/149‑1 (MQ=255)
                                                                                                           ccctgcgtcATGCGCAAAAGCCAGGCCCGGAGCAACGCCCGGCGAAGGGACGTCAATGGCGCGGAAGCCCTGCTCCGGGCGATCTTTAGCCACTGCCACCAGCGCCGCCCCCATGCCAACACCAATCTCAAGCGTCACCGGCGCTTCAc                                       <  2:641300/142‑1 (MQ=255)
                                                                                                           aCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCAc                                       <  2:627043/149‑1 (MQ=255)
                                                                                                            ccTTCTTCCTGCTCAGAAGCCAGGCACGCACCAACGCCAGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAACCCAATCTCAAGCGTCACCGGCGCTTCACg                                      <  2:475744/149‑1 (MQ=255)
                                                                                                            ccTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCAGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACg                                      <  2:475745/149‑1 (MQ=255)
                                                                                                            ccTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCAGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACg                                      <  2:475854/149‑1 (MQ=255)
                                                                                                                    aTGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCaaaaa                              <  1:355372/149‑1 (MQ=255)
                                                                                                                              gCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGACACCCTGACGGAACCAATCTCAAGCGTCCCCGGCGCTTCACGGCCAAAAAGCGCGGGGaa                    >  1:328865/1‑149 (MQ=255)
                                                                                                                                   gcacgcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCa               >  1:123144/1‑149 (MQ=255)
                                                                                                                                     acgcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCACGCCAAAACCAACCTCAAGCGTCACCGGCGCTTCACCGACCAAAAGCGCGGGGACATCCAGc             >  2:433773/1‑149 (MQ=255)
                                                                                                                                         acCAACGCCCGGGGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGGTCTTTAGCCAGTGCCACCAGCGACGCCCCCATGCCAAAACCACTCACAAGCGTCACGGGCGCTTCACGGCCAAAAAGCGCGGGGAGACCCAGCagat         >  2:359294/1‑149 (MQ=255)
                                                                                                                                               gCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGc   <  2:123562/149‑1 (MQ=255)
                                                                                                                                                cccGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCt  <  1:675113/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCCACGGGTCAGGGAAAAAGAGCTGCACCATGCGCAATGAATTGTCAGGAATCATTTTATGCAGCACTTCAACCGCATCGTGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAAAAAGCGCGGGGAAATCCAGCATATCTTCGCT  >  NC_000913/3102411‑3102703

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