Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,088,095 T→A 100% intergenic (+235/‑189) 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,0950TA100.0% 74.8 / NA 21intergenic (+235/‑189)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (8/13);  total (8/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAA  >  NC_000913/3087820‑3088266
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                           
aCAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAAtt                                                                                                                                                                 >  2:330733/1‑288 (MQ=255)
              ctgcGCGATGCTGCCGTTCTGAAGTAATCTTTCTTCGCCAGCGTTCAAAGGCCAGCCTCGCGCTGGCCTGTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGt                                                                                                                                                   <  2:204518/288‑1 (MQ=255)
                         tgcCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTatga                                                                                                                                        <  1:481097/288‑1 (MQ=255)
                                cTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCGCGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATaa                                                                                                                                <  2:156719/289‑1 (MQ=255)
                                                       cggtCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACt                                                                                                            <  1:303737/283‑1 (MQ=255)
                                                       cGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGata                                                                                                                                   >  1:468442/1‑263 (MQ=255)
                                                       cGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGata                                                                                                                                   <  2:468442/263‑1 (MQ=255)
                                                       cGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTaa                                                                                                         >  1:396134/1‑289 (MQ=255)
                                                         ttCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAAt                                                                                                       >  2:534519/1‑289 (MQ=255)
                                                                 ccAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACt                                                                                                <  2:580960/288‑1 (MQ=255)
                                                                 ccAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAAAAGCGTTTAATCTATGATGATATAAATCAATTATTTCCATGCACTTAAATCATAAcct                                                                                               >  2:303737/1‑287 (MQ=255)
                                                                           gcTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGt                                                                                     <  1:330733/289‑1 (MQ=255)
                                                                                                               aTCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTcac                                                 <  1:658384/289‑1 (MQ=255)
                                                                                                                                     acaATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCAc                            <  1:175904/288‑1 (MQ=255)
                                                                                                                                               catTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCATATCTTTTTACGTCGTACTCACCTATCTTAAtt                 >  1:530252/1‑289 (MQ=255)
                                                                                                                                                 ttcttcCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATAATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTc                >  1:658084/1‑288 (MQ=255)
                                                                                                                                                            ccagTTGCACCGCAGATTATCAGCACAACTGACACCGATTACACCAACCACCACAGACAAAGATTGGTAATATTTACATATTATTATTCGGTTTTCACGGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCGATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAAATAAATGTTAGTGTAAGAGATTACACTTTTGTAATTTGTTTCACATCTTTTTACGGCGTACTCACCTGTCTTAATTCACAATAAAAAAt    <  2:36995/286‑1 (MQ=255)
                                                                                                                                                            cACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAAt    <  1:534519/289‑1 (MQ=255)
                                                                                                                                                              cGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa  <  2:396134/289‑1 (MQ=255)
                                                                                                                                                                                         tGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCAcc                           >  2:246763/1‑237 (MQ=255)
                                                                                                                                                                                         tGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGAAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCAcc                           <  1:246763/237‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                           
ACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAGATTTGTAATATTTTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAA  >  NC_000913/3087820‑3088266

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