Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,947,466 T→C 19.3% intergenic (+3/‑31) 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,4660TC19.3% 38.1 / 4.3 27intergenic (+3/‑31)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (6/13);  new base C (2/2);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACG  >  NC_000913/2947319‑2947602
                                                                                                                                                   |                                                                                                                                        
aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCATCCAAtt                                                                                                                                          >  1:168229/1‑148 (MQ=255)
    ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                     <  7:190241/149‑1 (MQ=255)
     cGTGATGGGGGTAAAAACTAGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCTCGGACGCTGGGTGGAGCAGCCTGTTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                                    <  2:236745/149‑1 (MQ=21)
      ttGATCGGGGTAATAAATAGGTTGCATGATAACGCGAGGTGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                   <  4:71142/148‑1 (MQ=21)
         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                <  8:54046/149‑1 (MQ=255)
         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                <  8:189911/149‑1 (MQ=255)
            ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                             <  2:236910/149‑1 (MQ=255)
             gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTATGAAGt                                                                                                                           >  6:289443/1‑149 (MQ=255)
                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                                     <  7:200882/149‑1 (MQ=255)
                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGCAGCAGCCTGGTAGCTCGTCGGGCTCATAACCAGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAATTTGAATAATTaaaag                                                                                                                      >  8:351034/1‑145 (MQ=18)
                                 gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                                                       >  8:215897/1‑118 (MQ=255)
                                 gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                                                       <  7:215897/118‑1 (MQ=255)
                                  aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgcg                                                                                                        <  1:175856/148‑1 (MQ=255)
                                        cgAGCGGAGTATAGTGCGCATCCACGGACGCGTGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                                                                 <  1:27914/149‑1 (MQ=255)
                                            cGGAGTATCCTGCGCATCCACGGACTCGGTGAGGAGCATCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTAGGTTCTTATAATTACCGCGGAACCAATAAAATTTTGATGAATATACGCACTACCCTGACGCGGGGTGGagca                                                                                              <  2:168229/148‑1 (MQ=9)
                                                 tataGTGCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                                                                                <  2:330856/108‑1 (MQ=11)
                                                 tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                    <  5:177523/105‑1 (MQ=255)
                                                 tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                    >  6:177523/1‑105 (MQ=255)
                                                 tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                    <  2:213352/105‑1 (MQ=255)
                                                 tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                                    >  1:213352/1‑105 (MQ=255)
                                                 tataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTAGGTTCAAATCCGGCCCCAGCAACCAATTAAAATTtgatgta                                                                                                                             >  2:264623/1‑110 (MQ=255)
                                                                                        tAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGTAACCAATAAAAATTTGATCAATTAAATCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                                          >  6:188728/1‑108 (MQ=14)
                                                                                         aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAATAAATTAAAGCAGTACGGTGACGCGGGGTGGagcagc                                                                                            >  2:240544/1‑105 (MQ=14)
                                                                                                                      tcgGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagc                                                                                               >  6:284859/1‑73 (MQ=255)
                                                                                                                                             aCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACTCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGATCATAACCCGAAGGTCGTCTGTTCAAATCCGGCCACCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg  <  7:64756/143‑1 (MQ=18)
                                                                                                                                             aCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCAGGACAAAGCAACAAATCAAATTTGATGAAGTAAAAGCAGTACg  >  8:64756/1‑143 (MQ=18)
                                                                                                                                                aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAttatt                         <  4:29397/117‑5 (MQ=25)
                                                                                                                                                aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAttatt                         >  3:29397/1‑113 (MQ=25)
                                                                                                                                                   |                                                                                                                                        
AAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACG  >  NC_000913/2947319‑2947602

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