Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,425,164 C→A A24S (GCG→TCG)  hisM ← lysine/arginine/ornithine ABC transporter/histidine ABC transporter, membrane subunit HisM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,425,1640CA97.3% 122.5 / NA 37A24S (GCG→TCG) hisMlysine/arginine/ornithine ABC transporter/histidine ABC transporter, membrane subunit HisM
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (16/20);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACCGCGAAAAATATAGGTAAATAACCAGATTGGAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACA  >  NC_000913/2425027‑2425299
                                                                                                                                         |                                                                                                                                       
cgTACCGCGAAAAATATAGGTAAATAACCAGATTGGAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGAca                                                                                                                                       >  1:564779/1‑140 (MQ=255)
                      aaTAACCAGATTGGAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAGACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTc                                                                                                                 >  2:274806/1‑140 (MQ=255)
                             aGATTGGAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCcac                                                                                                          >  1:466469/1‑140 (MQ=255)
                              gATTGGAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCcaca                                                                                                         >  1:254905/1‑140 (MQ=255)
                                  ggAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTGCCGACAGAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCAcagcag                                                                                                     <  1:276302/140‑1 (MQ=255)
                                    aaaCTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCg                                                                                                   <  2:570906/140‑1 (MQ=255)
                                      accggATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGt                                                                                                 >  1:213952/4‑140 (MQ=255)
                                      aCTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGt                                                                                                 <  2:353450/140‑1 (MQ=255)
                                       cTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGtt                                                                                                <  2:248116/140‑1 (MQ=255)
                                        tGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGttt                                                                                               >  2:431560/1‑140 (MQ=255)
                                        tGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGttt                                                                                               >  1:198006/1‑140 (MQ=255)
                                          gATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTcc                                                                                             >  2:452141/1‑140 (MQ=255)
                                               atttattACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTAt                                                                                        <  1:233611/140‑1 (MQ=255)
                                               atttattACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTAt                                                                                        <  1:210261/140‑1 (MQ=255)
                                                    ttaCTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATg                                                                                   >  2:136966/1‑140 (MQ=255)
                                                    ttaCTGGAGACACGACCAATCGCCAGAAACAACGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATg                                                                                   <  2:250351/140‑1 (MQ=255)
                                                      aCTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGCTTCCAGTATTCATGTa                                                                                 >  2:154981/1‑140 (MQ=255)
                                                                  aCCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATca                                                                     <  1:194192/140‑1 (MQ=255)
                                                                            aGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGcc                                                           <  1:406928/140‑1 (MQ=255)
                                                                             gAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGccc                                                          >  2:314478/1‑140 (MQ=255)
                                                                                   gcgcCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGATACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTca                                                    <  1:475484/140‑1 (MQ=255)
                                                                                    cgcCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTcac                                                   <  1:408474/140‑1 (MQ=255)
                                                                                         ggACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCa                                              <  1:431560/140‑1 (MQ=255)
                                                                                               cGCCTATCACTACCGACAAAATAAGCAGACACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGt                                        >  2:259400/1‑140 (MQ=255)
                                                                                                   tATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAgcg                                    <  2:483116/140‑1 (MQ=255)
                                                                                                    aTCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAgcgg                                   >  2:574423/1‑140 (MQ=255)
                                                                                                     tCACTACCGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAgcggc                                  <  2:466737/140‑1 (MQ=255)
                                                                                                           ccGACAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCaa                            >  1:96846/1‑140 (MQ=255)
                                                                                                               cAAAATAAGCAGCCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGaa                        <  1:136966/140‑1 (MQ=255)
                                                                                                                       gcagcCACAGCGTGATCGTCCAACTAGTAAAGCGATAACCGTCGGTCCACAGCCTCGGTTCCCACTATTCATGTAAGATCTCGATCACAAGTCAGCCCACTTCTCACCCGCGGAGTAGCGGCGCTCAAGGAACAgcagca                <  2:198006/140‑1 (MQ=255)
                                                                                                                       gcagcCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAgcagca                <  2:269394/140‑1 (MQ=255)
                                                                                                                       gcagcCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAgcagca                <  2:96846/140‑1 (MQ=255)
                                                                                                                         agcCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGcaca              <  2:45637/140‑1 (MQ=255)
                                                                                                                         agcCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGcaca              <  2:25631/140‑1 (MQ=255)
                                                                                                                         agcCACAGAGTGATCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGcaca              <  1:314478/140‑1 (MQ=255)
                                                                                                                                     aTCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACa  >  2:254804/1‑140 (MQ=255)
                                                                                                                                     aTCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACa  <  2:564779/140‑1 (MQ=255)
                                                                                                                                     aTCGACACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACa  >  1:256/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CGTACCGCGAAAAATATAGGTAAATAACCAGATTGGAAACTGGATGTATTTATTACTGGAGACACGACCAATCGCCAGAAACAGCGCCAGGACTCCGCCTATCACTACCGACAAAATAAGCAGCCACAGAGTGATCGCCACACCAGTAAAGCGATAACCGTCGGTCCACAGCAGCGGTTTCCAGTATTCATGTAAGATTTCGATCACAGGTCAGCCCTCTTCACACCCACGGAGTAGCGGCGCTCAAGGAACAGCAGCACACCATTGGAAACA  >  NC_000913/2425027‑2425299

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: