Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,641,025 Δ1 bp coding (84/687 nt) yjtD → putative rRNA methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,641,0230T.97.2% 139.4 / ‑1.8 36coding (82/687 nt)yjtDputative rRNA methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (16/19);  total (16/20)
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.25e-01

GTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGA  >  NC_000913/4640891‑4641150
                                                                                                                                    |                                                                                                                               
gTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGt                                                                                                                                 >  1:554249/1‑133 (MQ=255)
          aTCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTg                                                                                                                      <  2:25456/133‑1 (MQ=255)
           tCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGc                                                                                                                     >  2:515882/1‑133 (MQ=255)
              cccGGTTCTCGGTTATTTTATCAAGAAGAGCGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGa                                                                                                                  <  1:236936/133‑1 (MQ=255)
                  gTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGAt                                                                                                                 >  2:62605/1‑130 (MQ=255)
                  gTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGAt                                                                                                                 <  1:62605/130‑1 (MQ=255)
                  gTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGt                                                                                                              >  2:411456/1‑133 (MQ=255)
                  gTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAAGGAAGACGATGGGGTTAACGATCTGCGGATTGt                                                                                                              >  2:236936/1‑133 (MQ=255)
                         gTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGt                                                                                                       <  1:419501/133‑1 (MQ=255)
                                      gaagaGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGca                                                                                                 >  1:480084/1‑126 (MQ=255)
                                      gaagaGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGca                                                                                                 <  2:480084/126‑1 (MQ=255)
                                            tGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGagccag                                                                                    >  1:538701/1‑133 (MQ=255)
                                                     gCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCt                                                                           <  1:458367/133‑1 (MQ=255)
                                                        tataACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTggg                                                                        <  1:515882/133‑1 (MQ=255)
                                                         ataACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTggg                                                                        <  2:499724/133‑1 (MQ=255)
                                                          taACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTc                                                                      <  2:44073/133‑1 (MQ=255)
                                                          taACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTc                                                                      <  2:554249/133‑1 (MQ=255)
                                                           aaCGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCt                                                                           >  1:383129/1‑127 (MQ=255)
                                                           aaCGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGACAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCt                                                                           <  2:383129/127‑1 (MQ=255)
                                                                   ttCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGcaca                                                                 <  2:592614/129‑1 (MQ=255)
                                                                   ttCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGcaca                                                                 >  1:592614/1‑129 (MQ=255)
                                                                    tCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGAt                                                            <  1:115000/133‑1 (MQ=255)
                                                                    tCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGCATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGAt                                                            <  1:529760/133‑1 (MQ=255)
                                                                        gTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCTGCCACCCGCTGGGTCGCACATGGATCTgg                                                        >  2:516763/1‑133 (MQ=255)
                                                                        gTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTgg                                                        >  1:151015/1‑133 (MQ=255)
                                                                          cGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTgg                                                        >  2:414020/1‑131 (MQ=255)
                                                                          cGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTgg                                                        <  1:414020/131‑1 (MQ=255)
                                                                             aCCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGata                                                   >  1:553314/1‑133 (MQ=255)
                                                                             aCCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGata                                                   >  1:516969/1‑133 (MQ=255)
                                                                             aCCCGCCAGAGCAGAAAATATTGGAGCAACGGTGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGata                                                   <  2:183862/133‑1 (MQ=255)
                                                                                cGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGAtatta                                                <  2:553314/133‑1 (MQ=255)
                                                                                    agagCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGa                                            <  1:551294/133‑1 (MQ=255)
                                                                                     gagCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGat                                           <  1:411456/133‑1 (MQ=255)
                                                                                                         gcggcgCGGGCAATGAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCTGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGAc                       >  1:545282/1‑133 (MQ=255)
                                                                                                               cgGGCAATAAAAACGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGc                 >  2:121290/1‑133 (MQ=255)
                                                                                                                      taaaaaCGATGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCg          <  1:121290/131‑1 (MQ=255)
                                                                                                                              aTGGGGTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGa  >  2:322217/1‑133 (MQ=255)
                                                                                                                                    |                                                                                                                               
GTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGATATTATTGATAATATTAAAGTTTTCCCGACATTGGCTGAATCGTTACACGA  >  NC_000913/4640891‑4641150

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: