Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,399,552 A→G I101V (ATC→GTC)  miaA → tRNA dimethylallyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,399,5520AG97.7% 151.3 / ‑4.0 43I101V (ATC→GTC) miaAtRNA dimethylallyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (9/33);  total (9/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

CGAAGCCGAACGCTGAAGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCA  >  NC_000913/4399415‑4399688
                                                                                                                                         |                                                                                                                                        
cGAAGCCGAACGCTGAAGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTc                                                                                                                                        >  2:213541/1‑140 (MQ=255)
 gAAGCCGAACGCTGAAGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTcc                                                                                                                                       <  1:280885/140‑1 (MQ=255)
         aCGCTGAAGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTa                                                                                                                               <  2:499856/140‑1 (MQ=255)
           gCTGAAGAGTCACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGAAATCACCGCGGCGGGGCGGATCCCACTGTTAGt                                                                                                                             <  1:1058769/140‑1 (MQ=255)
             tGAAGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTgg                                                                                                                           <  1:371232/140‑1 (MQ=255)
               aaGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCCCCGTTTGTgggg                                                                                                                         >  2:592532/1‑139 (MQ=255)
                agagTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCg                                                                                                                        <  2:992573/140‑1 (MQ=255)
                 gagTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCgg                                                                                                                       <  2:1103418/140‑1 (MQ=255)
                   gTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTa                                                                                                                     <  2:507355/140‑1 (MQ=255)
                   gTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTa                                                                                                                     <  1:435268/140‑1 (MQ=255)
                    ttACTCGCCGCGCCGCACCGATTGTTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGGAc                                                                                                                    >  1:905893/1‑140 (MQ=255)
                    ttACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCACAGTTAATGGGGGGGAc                                                                                                                    >  1:1095313/1‑140 (MQ=255)
                     tACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACg                                                                                                                   <  1:110146/140‑1 (MQ=255)
                      aCTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGa                                                                                                                  <  1:413735/140‑1 (MQ=255)
                      aCTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGa                                                                                                                  <  1:629497/140‑1 (MQ=255)
                      aCTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGGCCCCTGTTTGTggggggggcga                                                                                                                  >  1:586237/1‑135 (MQ=255)
                        tCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGAtg                                                                                                                <  1:225352/140‑1 (MQ=255)
                         cgccgcgccgCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCCCCGTTTGTGGGGGGGGCGAtgt                                                                                                               >  2:1071933/1‑140 (MQ=255)
                           ccgcgccgcACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGAtgttg                                                                                                             <  2:577417/140‑1 (MQ=255)
                            cgcgccgcACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCCCTGTTTGTGGGGGGGACGAtgttgt                                                                                                            >  1:740480/1‑140 (MQ=255)
                            cgcgccgcACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGAtgttgt                                                                                                            <  1:43819/140‑1 (MQ=255)
                            cgcgccgcACCGATTGCTGGATATTCGCGATCCGTCACAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTa                                                                                                                     >  1:523500/1‑131 (MQ=255)
                            cgcgccgcACCGATTGCTGGATATTCGCGATCCGTCACAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTa                                                                                                                     <  2:523500/131‑1 (MQ=255)
                             gcgccgcACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCCCCGTTAGGTGGGGGGGCGAtgttgtt                                                                                                           >  2:435268/1‑139 (MQ=255)
                               gccgcACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCCCCGGTAGTGGGGGGGACCAAGGTGTAtt                                                                                                         >  2:43819/1‑140 (MQ=255)
                                    aCCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAg                                                                                                    <  2:325443/140‑1 (MQ=255)
                                     ccGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAgg                                                                                                   <  1:213541/140‑1 (MQ=255)
                                     ccGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAgg                                                                                                   <  1:1117367/140‑1 (MQ=255)
                                     ccGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAgg                                                                                                   <  2:1123399/140‑1 (MQ=255)
                                       gATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGTGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCa                                                                                                 <  1:27107/140‑1 (MQ=255)
                                        aTTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCAt                                                                                                <  1:1130397/140‑1 (MQ=255)
                                                                        aCTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCa                                                                <  2:848182/140‑1 (MQ=255)
                                                                          tCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGa                                                              <  1:684748/140‑1 (MQ=255)
                                                                              ctgctgATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCCCTGTTTGTGGGGGGGACGATGTTGTATTTTAAAGCATTTCTGGGAGGGGTGTCGCCCCCACCGTCGGCAGAcccc                                                          >  2:8430/1‑139 (MQ=255)
                                                                              ctgctgATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGGCCCCCTGTTTGTGGGGGGGACGAAGTTGTATTTTAAAGGATTTCTGGGAGGGGTGTTGCCGCCACCCTCGGGAGAaccc                                                          >  2:463432/1‑139 (MQ=255)
                                                                                    aTTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTa                                                    <  1:592532/140‑1 (MQ=255)
                                                                                    aTTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTa                                                    <  1:463432/140‑1 (MQ=255)
                                                                                           ccgcGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCa                                             <  1:8430/140‑1 (MQ=255)
                                                                                                             aaTGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAggcgg                           <  2:307482/140‑1 (MQ=255)
                                                                                                             aaTGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAggcgg                           <  2:586237/140‑1 (MQ=255)
                                                                                                             aaTGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAggcgg                           <  1:768103/140‑1 (MQ=255)
                                                                                                             aaTGGCCGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAggcgg                           <  2:944983/140‑1 (MQ=255)
                                                                                                                  ccGATATCACCGCGGCGGGGCGGGGCCCCCTGGTTGTGGGGGGGACGGTGTTGTATTTTAAGGGATTGCTGGGAGGGGTGTCGCCGCCACCCTCGGGCGACCCGGGAGTACGGGCCCGAATTGAGCAAAAGGGGGGagag                      >  2:69824/1‑140 (MQ=255)
                                                                                                                   cGATATCACCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGc                     <  2:963594/140‑1 (MQ=255)
                                                                                                                          aCCGCGGCGGGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTg              <  1:220640/140‑1 (MQ=255)
                                                                                                                            cgcgGCGGGGCGGGGCCCCCCGTTTGTGGGGGGGACGATGTTGTATTTTAAGGCCTTGCTGGAAGGGGTGTCGCCCCCACCGGCGGCAGACCCGGAAGTACGGGGCAGAAATTAGCAACAAGCGGCAGAGCAAAGGTggg            >  2:425602/1‑140 (MQ=255)
                                                                                                                              cggcggGGCGGGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAg          <  1:214451/140‑1 (MQ=255)
                                                                                                                                      cggGTCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCa  <  1:1071933/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CGAAGCCGAACGCTGAAGAGTTACTCGCCGCGCCGCACCGATTGCTGGATATTCGCGATCCGTCGCAGGCTTACTCGGCTGCTGATTTTCGCCGCGATGCGCTGGCGGAAATGGCCGATATCACCGCGGCGGGGCGGATCCCACTGTTAGTGGGCGGTACGATGTTGTATTTCAAGGCATTGCTGGAAGGGTTGTCGCCGCTACCGTCGGCAGACCCGGAAGTACGGGCCAGAATTGAGCAACAGGCGGCAGAGCAAGGTTGGGAGTCATTGCA  >  NC_000913/4399415‑4399688

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: