Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,940,308 A→G L102P (CTG→CCG)  hsrA ← putative multidrug or homocysteine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,940,3080AG100.0% 115.5 / NA 32L102P (CTG→CCG) hsrAputative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

CACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTG  >  NC_000913/3940175‑3940454
                                                                                                                                     |                                                                                                                                                  
caccaGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAAtt                                                                                                                                    >  2:132747/1‑150 (MQ=255)
                gCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGcc                                                                                                                    >  1:241588/1‑150 (MQ=255)
                  gACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCaga                                                                                                                 <  2:431972/151‑1 (MQ=255)
                     aaaaTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGaa                                                                                                              <  1:247096/151‑1 (MQ=255)
                        aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                           <  1:220346/151‑1 (MQ=255)
                        aTTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAAccc                                                                                                           <  1:624810/151‑1 (MQ=255)
                            aaTACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtg                                                                                                       >  2:309104/1‑151 (MQ=255)
                               aCTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAAtgtg                                                                                                     >  2:552140/1‑150 (MQ=255)
                                     agaagTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGAc                                                                                              <  2:408011/151‑1 (MQ=255)
                                               tACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCaa                                                                                     <  2:491878/150‑1 (MQ=255)
                                                aCGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAgg                                                                                   <  1:613170/151‑1 (MQ=255)
                                                  gAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAggg                                                                                  >  2:314416/1‑150 (MQ=255)
                                                   aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGt                                                                                 >  2:327686/1‑150 (MQ=255)
                                                   aGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAAGTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTa                                                                                <  1:412416/151‑1 (MQ=255)
                                                     gATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaa                                                                               >  2:626346/1‑150 (MQ=255)
                                                      aTAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTaaa                                                                              >  1:62001/1‑150 (MQ=255)
                                                          gcgcgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGc                                                                         >  1:171298/1‑151 (MQ=255)
                                                             cgcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAc                                                                      >  1:34636/1‑151 (MQ=255)
                                                              gcAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGAcg                                                                     >  2:442292/1‑151 (MQ=255)
                                                                       ccAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGaa                                                            <  1:275970/151‑1 (MQ=255)
                                                                          gCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCg                                                         <  2:174149/151‑1 (MQ=255)
                                                                              aGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCg                                                     <  2:39103/151‑1 (MQ=255)
                                                                                 aaCAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCg                                                     >  1:357124/1‑148 (MQ=255)
                                                                                        atcatcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATc                                           <  2:212944/151‑1 (MQ=255)
                                                                                          catcatTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAGCCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATcca                                         <  1:5398/151‑2 (MQ=255)
                                                                                                 gcgcCGCCTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCg                                  <  2:53011/151‑1 (MQ=255)
                                                                                                      gccTATCCCCTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAAt                              >  2:316488/1‑150 (MQ=255)
                                                                                                             cccTGAATAACCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCg                      >  2:5158/1‑151 (MQ=255)
                                                                                                                     aaCCCGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTc              >  1:548715/1‑151 (MQ=255)
                                                                                                                        ccGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGc           <  1:608963/151‑1 (MQ=255)
                                                                                                                         cGGAAGACAACCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCg          <  1:132743/151‑1 (MQ=255)
                                                                                                                                 aaCCGGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACtg  >  2:640314/1‑151 (MQ=255)
                                                                                                                                     |                                                                                                                                                  
CACCAGACCCGGCATGGCGACAAAATTCAATACTGGAAGAAGTTCATTACGAGGATAAGCGCGCAGTAAGGCCAGCCGAGCAACAGGCATCATCATTGCGCCGCCTATCCCCTGAATAACCCGGAAGACAACCAGCTGTGGTAGCGAATTAGAAAGTGCGCAGGCCAGAGAACCCAATGTGAACAGACTCACGGCAAGGGTAAAAATGCGACGCGTACCGAAGCGATCGGCTAGCCATCCGCTTACCGGAATAAGCATCGCCACCGTCAGCGTATAACTG  >  NC_000913/3940175‑3940454

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