Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,088,189 T→C intergenic (+329/‑95) metK → / → galP S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,1890TC100.0% 113.8 / NA 31intergenic (+329/‑95)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/9);  total (22/9)

ATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCG  >  NC_000913/3088047‑3088335
                                                                                                                                              |                                                                                                                                                  
atTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGAt                                                                                                                                             <  1:652860/150‑1 (MQ=255)
  tttCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTa                                                                                                                                           <  2:588631/150‑1 (MQ=255)
   ttCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTGTTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTac                                                                                                                                          >  1:619905/1‑150 (MQ=255)
     cATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACt                                                                                                                                       >  2:13456/1‑151 (MQ=255)
                         acaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATc                                                                                                                   >  2:514391/1‑151 (MQ=255)
                               gttACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCttttt                                                                                                              >  2:480139/1‑150 (MQ=255)
                                               gTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTAt                                                                                             >  2:156175/1‑151 (MQ=255)
                                                 aaaGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCt                                                                                           >  1:423195/1‑151 (MQ=255)
                                                  aaGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCtt                                                                                          <  1:132536/151‑1 (MQ=255)
                                                     tCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAAt                                                                                       >  2:534307/1‑151 (MQ=255)
                                                        tAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAAt                                                                                       <  1:374899/148‑1 (MQ=255)
                                                         aaTTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTcac                                                                                   >  2:61199/1‑151 (MQ=255)
                                                          aTTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTcac                                                                                   >  1:40631/1‑150 (MQ=255)
                                                            tGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACaa                                                                                 <  2:497096/150‑1 (MQ=255)
                                                            tGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAAt                                                                                >  2:146708/1‑151 (MQ=255)
                                                              cAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTTAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATa                                                                               >  1:357228/1‑150 (MQ=255)
                                                              cAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATaa                                                                              <  2:591947/151‑1 (MQ=255)
                                                                           ttAATCTATGATGATATAACTCAATTATTTTCATGCTCTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAACTCACAATAAAAAATAACCata                                                                  <  2:142457/150‑1 (MQ=255)
                                                                             aaTCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTg                                                               >  1:467780/1‑151 (MQ=255)
                                                                                  atgatgATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAgg                                                           <  2:661332/150‑1 (MQ=255)
                                                                                          ataACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGc                                                  >  2:138957/1‑151 (MQ=255)
                                                                                                      ttttCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTaaa                                       >  2:181506/1‑150 (MQ=255)
                                                                                                       tttCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTaaaa                                      >  1:262934/1‑150 (MQ=255)
                                                                                                         tCATGCACTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAAc                                   >  1:6630/1‑151 (MQ=255)
                                                                                                               aCTTAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGc                             >  1:620420/1‑151 (MQ=255)
                                                                                                                  tAAATCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCgg                           <  1:272186/150‑1 (MQ=255)
                                                                                                                     aTCATAACTAAGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCa                        >  1:571866/1‑150 (MQ=255)
                                                                                                                              aaGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCa               >  2:451658/1‑150 (MQ=255)
                                                                                                                               aGATAAATGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAAt             >  2:159370/1‑151 (MQ=255)
                                                                                                                                    aaTGTTAGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGt        >  1:437035/1‑151 (MQ=255)
                                                                                                                                          aGTGCAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGCCAAACAAGGCAATGACGTTTTTCg  >  1:351899/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCG  >  NC_000913/3088047‑3088335

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