Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 510,729 G→A S217S (AGC→AGT copA ← Cu(+) exporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913510,7290GA96.3% 85.3 / NA 27S217S (AGC→AGTcopACu(+) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (14/12);  total (15/12)
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

ACGCCAGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGGCTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAGTGCGACAATTGCCTGCCAGCGGAAGCGCTTCATCGTAGCGACGGCGGT  >  NC_000913/510595‑510851
                                                                                                                                      |                                                                                                                          
acgccaGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCg                                                                                                                       <  1:169779/140‑1 (MQ=255)
 cgccaGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACt                                                                                                                          <  2:597774/136‑1 (MQ=255)
 cgccaGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACt                                                                                                                          >  1:597774/1‑136 (MQ=255)
 cgccaGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCgg                                                                                                                      >  1:133995/1‑140 (MQ=255)
     aGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACTAACCACAGACTCCGGTTGt                                                                                                                  <  2:421826/140‑1 (MQ=255)
              cgcCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAAACAACCAAAGACTGGGGTTGTCAGCGGTGa                                                                                                         >  1:353876/1‑140 (MQ=255)
                 caccaGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGATTATCAGGTCGAAAACTAACCACAGGCTGCCGTTGTCAGCGGTGACca                                                                                                      >  2:221580/1‑140 (MQ=255)
                  accaGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACcat                                                                                                     <  2:463239/140‑1 (MQ=255)
                           tCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATc                                                                                            >  1:477442/1‑140 (MQ=255)
                              aTCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGcc                                                                                         <  2:543660/140‑1 (MQ=255)
                                   cgccgcACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGatca                                                                                    <  2:159368/140‑1 (MQ=255)
                                   cgccgcACCGTTCAGCAGGCTTGTCAATGAACTGCGGAAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGatca                                                                                    <  1:206893/140‑1 (MQ=255)
                                            gTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGa                                                                           >  2:582581/1‑140 (MQ=255)
                                             ttCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGAc                                                                          <  2:477442/140‑1 (MQ=255)
                                                                   tGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGcc                                                    >  2:702285/1‑140 (MQ=255)
                                                                    gCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCGTCCCCCAGACCATCACCGGGa                                                               <  2:410806/128‑1 (MQ=255)
                                                                    gCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGa                                                               >  1:410806/1‑128 (MQ=255)
                                                                     cGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGa                                                               <  2:145698/127‑1 (MQ=255)
                                                                     cGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGa                                                               <  2:682180/127‑1 (MQ=255)
                                                                     cGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGa                                                               >  1:682180/1‑127 (MQ=255)
                                                                     cGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGa                                                               >  1:145698/1‑127 (MQ=255)
                                                                     cGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAg                                                  >  2:609023/1‑140 (MQ=255)
                                                                     cGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAg                                                  >  2:610784/1‑140 (MQ=255)
                                                                                             aCCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAGTGCGACAATTGCCTGCCAGCGGaa                          <  1:702285/140‑1 (MQ=255)
                                                                                                       ccAGGGTTATCAGGCCGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAGTGCGACAATTGCCTGCCAGCGGAAGCGCTTCATc                >  1:28380/1‑140 (MQ=255)
                                                                                                                     ccGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAGTGCGACAATTGCCTGCCAGCGGAAGCGCTTCATCGTAGCGACGGCGGt  >  1:526378/1‑140 (MQ=255)
                                                                                                                     ccGATAACCAACCACAGACTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAGTGCGACAATTGCCTGCCACCGGCAGTGCTTCATCGTAGCGACGGCGGt  >  1:359182/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                          
ACGCCAGTACCCAGCGCCACCAGCGTATCCATCGTCGCCGCACCGTTCAGCAGGCTTTTCCATGCACTGCGGTAAAAATGGCCGCCGGCGAAAACCATCACTGCCAGGGTTATCAGGCCGATAACCAACCACAGGCTGCGGTTGTCAGCGGTGACCATCATGTTATCGCCGATCATCCCCCAGACCATCACCGGGATACCCACCGCCAGTGCGACAATTGCCTGCCAGCGGAAGCGCTTCATCGTAGCGACGGCGGT  >  NC_000913/510595‑510851

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