Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,947,465 Δ1 bp 11.8% intergenic (+2/‑32) metZ → / → metW tRNA‑Met/tRNA‑Met

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4650T.11.8% 98.4 / 7.4 34intergenic (+2/‑32)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (7/23);  new base . (2/2);  total (9/25)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCA  >  NC_000913/2947325‑2947597
                                                                                                                                            |                                                                                                                                    
gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                              <  4:230501/149‑1 (MQ=255)
   aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                            <  2:183632/148‑1 (MQ=255)
          tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTa                                                                                                                     >  7:120562/1‑148 (MQ=255)
            aaaaacaGTTTGCATTAAAACGCGAGAGGATTATATTGAGCATCAAATTACGCGGGGTGGAGCAGCCTGGTAGCACGTAGGGCTCATAACAGGAAGGTCGTCGGTTCAAATCCGGCCAACGCAACCAATTAAAATTTGATGAATTaaa                                                                                                                   <  5:372237/148‑1 (MQ=17)
            aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                   <  5:393772/148‑1 (MQ=255)
            aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                   <  3:137476/148‑1 (MQ=255)
            aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                   <  4:215109/148‑1 (MQ=255)
            aaaaaaaGGTGGCAGGAAAACGCGTGCGGAGTTTAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCTGGCTCATAACCCGAAGGTCGTCGGTTCAAATGCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                                   <  6:359623/148‑1 (MQ=18)
              aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCTCGGACGCGGGGTGTAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAGCAATTAAAATTTGATGAAGTAAAGCa                                                                                                                <  6:287780/149‑1 (MQ=255)
              aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                                <  5:26443/149‑1 (MQ=255)
                 aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                                             <  7:9253/149‑1 (MQ=255)
                   ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTAGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                           <  7:15531/149‑1 (MQ=255)
                   ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                           <  6:381567/149‑1 (MQ=255)
                   ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                           <  6:222846/149‑1 (MQ=255)
                    gTTTCATGAAAAGGGGATCGTAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGATGGTCGTCGGTTCAAATCCGGGCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                           <  1:63590/148‑1 (MQ=18)
                    gTTGCATGATAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCTGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACgg                                                                                                          <  5:407347/149‑1 (MQ=255)
                    gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACgg                                                                                                          <  1:77221/149‑1 (MQ=255)
                     ttGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                                                                         <  6:159459/149‑1 (MQ=255)
                         aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                            >  5:193044/1‑125 (MQ=255)
                         aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                                                                                                                            <  6:193044/125‑1 (MQ=255)
                          tGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                             >  1:368124/1‑124 (MQ=255)
                          tGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                             <  2:368124/124‑1 (MQ=255)
                             aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCggg                                                                                                 <  3:188007/149‑1 (MQ=255)
                               aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                               <  8:120562/149‑1 (MQ=255)
                                             taGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                                                                                                        <  7:97438/109‑1 (MQ=255)
                                             taGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                                                                                                        >  8:97438/1‑109 (MQ=255)
                                                      tCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                       <  7:404283/118‑1 (MQ=255)
                                                      tCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                       >  8:404283/1‑118 (MQ=255)
                                                             aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCgggg                                                                                                >  4:294926/1‑118 (MQ=255)
                                                                                        gtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                           <  7:195674/80‑1 (MQ=255)
                                                                                        gtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                           >  8:195674/1‑80 (MQ=255)
                                                                                                     cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                            >  4:45990/1‑82 (MQ=255)
                                                                                                                                      aaCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCTTCGGGCTCATAAACAGAAGGTCGTCGGTTCAAATACTGCACCCGCAACCACTCAAATTTGATGAAGTAAAAGCa  >  2:342727/1‑139 (MQ=18)
                                                                                                                                      aaCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCa  <  1:342727/139‑1 (MQ=255)
                                                                                                                                            ttAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCaaa                                                          >  2:331073/1‑77 (MQ=255)
                                                                                                                                            |                                                                                                                                    
GTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCA  >  NC_000913/2947325‑2947597

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