Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,947,713 (T)7→8 intergenic (+30/+44) metV → / ← amiC tRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,7061.T92.9% 84.4 / ‑0.0 28intergenic (+23/+51)metV/amiCtRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase C
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (15/11);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTG‑‑TTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCGCAAGAATAGACTCCGCAACTTCCTGCTG  >  NC_000913/2947599‑2947837
                                                                                                             |                                                                                                                                     
tACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc                                                                                                      >  1:249564/1‑140 (MQ=255)
       aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGTACACCCAAAAGGGGGATTTTTTTGATTCTGGGCTCCCATAAAAAAGAGCCCTTc                                                                                               >  2:599302/1‑140 (MQ=17)
       aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTc                                                                                               >  2:517517/1‑140 (MQ=255)
         gcgGGGTGGAGCAGTCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAg                                                                                             <  2:905115/140‑1 (MQ=255)
          cgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAgc                                                                                            >  1:654956/1‑140 (MQ=255)
          cgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTG‑‑TTTTTTTGTTTTTTGTCTTCCATTAAAAAAGCCCCAtcaacc                                                                                            >  1:316395/1‑135 (MQ=18)
           ggggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAgcg                                                                                           <  1:599302/140‑1 (MQ=255)
                 gAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTG‑‑TTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCttttt                                                                                    >  2:728251/1‑140 (MQ=255)
                   gcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTa                                                                                   <  1:791803/140‑1 (MQ=255)
                       ccTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAgc                                                                                            >  2:817734/1‑127 (MQ=255)
                       ccTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAgc                                                                                            <  1:817734/127‑1 (MQ=255)
                             aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCtt                                                                         <  1:379569/140‑1 (MQ=255)
                                 cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGATCACCCTAACGGGGGATTATTTTGATTCGGGTCGCCCCTAAAAAAGAGCCATTCAGCCCCTTATTATCATCCCCTtctgt                                                                     >  2:341871/1‑138 (MQ=21)
                                 cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTtctc                                                                      >  2:605053/1‑140 (MQ=255)
                                 cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCg                                                                     >  2:204414/1‑140 (MQ=255)
                                   tcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTTACACGCTGAAGGGGGA‑‑TTTTTTTGATGCTGGGCGCCCATAAAAGAGCGCCATTCAGCCCCATTTTATCATCCCCTTCTCGcc                                                                   >  2:379569/1‑140 (MQ=18)
                                               cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCa                                                       <  2:654956/140‑1 (MQ=255)
                                                   aaGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCgg                                                   >  1:909146/1‑140 (MQ=255)
                                                   aaGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCgg                                                   >  1:261597/1‑140 (MQ=255)
                                                   aaGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCgg                                                   >  1:358888/1‑140 (MQ=255)
                                                     ggTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCa                                                 <  2:249564/140‑1 (MQ=255)
                                                     ggTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCa                                                 >  2:595032/1‑140 (MQ=255)
                                                         gtcgGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAAt                                             <  1:341871/140‑1 (MQ=255)
                                                          tcgGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATa                                            >  1:566909/1‑140 (MQ=255)
                                                          tcgGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATa                                            >  1:519227/1‑140 (MQ=255)
                                                                        gCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCg                              <  2:358888/140‑1 (MQ=255)
                                                                        gCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCg                              >  1:618053/1‑140 (MQ=255)
                                                                                  ccAATTATTGAACACCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCGCAAGAATAGa                    <  1:595032/140‑1 (MQ=255)
                                                                                              cacCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCGCAAGAATAGACTCCGCAACTTc        >  2:575216/1‑140 (MQ=255)
                                                                                               acCCTAACGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCGCAAGAATAGACTCCGCAACTTcc       <  2:566909/140‑1 (MQ=255)
                                                                                                     aaCGGGTGTTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCACCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCGCAAGAATAGACTCCGCAACTTCctgctg  >  1:121146/1‑140 (MQ=255)
                                                                                                             |                                                                                                                                     
TACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTG‑‑TTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAATCCCCGCAAGAATAGACTCCGCAACTTCCTGCTG  >  NC_000913/2947599‑2947837

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: