Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,702,802 G→A P23S (CCG→TCG)  dppF ← dipeptide ABC transporter ATP binding subunit DppF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,8020GA87.5% 28.4 / ‑0.7 17P23S (CCG→TCG) dppFdipeptide ABC transporter ATP binding subunit DppF
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (5/9);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

CAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACG  >  NC_000913/3702611‑3703030
                                                                                                                                                                                               |                                                                                                                                                                                                                                    
cAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCAGTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTGACCGCAGCCAGAGTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGTCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGc                                                                                                                                                                                           <  1:172691/235‑1 (MQ=255)
                       gCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGcc                                                                                                                                                                                  <  1:13305/221‑1 (MQ=255)
                                      gCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATTCCTTTCTTCACCGAATAATGTTTTTTCAGGGCGATAGCCTGCAACAGCGGTGGTTGCAGGGTGGCCTCTTGTGTACTCATAGTGTCGGCGTCCCGGCATCATCAAGTGGGt                                                                                                                                     <  1:176823/251‑1 (MQ=255)
                                                         cAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGACCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAAACAGACGTTCCGGCGCGAACATGCCTTTCTTCCCCGCATAATGTTTTTTCAGGGCGATAGCCTGCCACAGCGGTTGTTGCAGGGTGGCCTCTTGCGAACTCATAGTGTCg                                                                                                                                                             >  1:246073/1‑208 (MQ=255)
                                                                    nggTCGATTTACCGCAGCAAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTggg                                                                                                                                      >  1:261763/2‑220 (MQ=255)
                                                                                  cAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAgcgcn                                                                                         <  1:80879/251‑2 (MQ=255)
                                                                                                        ccAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTAc                                                                                                                                                                        >  1:273045/1‑150 (MQ=255)
                                                                                                                       gTTCAAGGTTAAACGAAACGCCATCCAGCTCTTTAACCAGACGTTCCGGCGCGAACAGGGCTTTCTTCACCGAATAATGTTTTTTCAGGTCTATAGCCTTCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATan                                                                                                                                                                  <  1:65264/141‑2 (MQ=255)
                                                                                                                             ggTTACACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTc                                                                                                                                                                                                                    >  1:76190/1‑85 (MQ=255)
                                                                                                                                tAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATTGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCg                                                                                                                                                             <  1:122891/137‑1 (MQ=255)
                                                                                                                                                      tttAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTc                                                                                                                                                                      <  1:70496/106‑1 (MQ=255)
                                                                                                                                                         aaCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTg                                                                    <  1:91498/201‑1 (MQ=255)
                                                                                                                                                          aCCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATACTGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCCAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAg                                                        >  1:123830/1‑212 (MQ=255)
                                                                                                                                                           ccAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGAATAATGTTTTGTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAg                                                                                                                                          <  1:264097/129‑1 (MQ=255)
                                                                                                                                                                          cgAACATGCCTTTCTTCACCGAATAATGTTTTTTCAGGTCGATAGCCTGAAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCACGCGGGTAATGGCCTTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACg  >  1:218345/1‑250 (MQ=255)
                                                                                                                                                                                tGCCTTTCTTCACCGAAGAATGTTTTTTCAGGTCGATAGCCTGCaan                                                                                                                                                                                                       <  1:66378/47‑2 (MQ=255)
                                                                                                                                                                                tGCCTNNNTTNACNGNANAATGTTTTTTCAGGTCGATAGCCTGCAAc                                                                                                                                                                                                       >  2:66378/1‑47 (MQ=38)
                                                                                                                                                                                            ccGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAg                                                                                                                                                                  >  1:54643/1‑72 (MQ=255)
                                                                                                                                                                                               |                                                                                                                                                                                                                                    
CAGATCCTGCCCCTGGTAATACAGCTCGCCACCGGTGGGCATTTCAATCATCGTCAGCAACCGACCGAGGGTCGATTTACCGCAGCCAGATTCGCCCACTACTGCCAGCGTTTTGCCACGTTCAAGGTTAAACGAAACGCCATCCAGCGCTTTAACCAGACGTTCCGGCGCGAACATGCCTTTCTTCACCGGATAATGTTTTTTCAGGTCGATAGCCTGCAACAGCGGTTGTTGCAGGGTGGCCTCTTGCGTACTCATAGTGTCGGCCTCCCGGCATCATCAAGTGGGTAATGGCATTTGGACTGACGCCCGTCAGCGAGCATATTCAGCGCCGGTTCTTCAGCGCGACATCTGTCAGTGGCATAGGGGCAGCGCGGGTTAAGCAGGCAGCCGTTCGGGCGGTCGTACTTGCCGGGAACG  >  NC_000913/3702611‑3703030

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