Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,451,483 G→A 17.4% P142P (CCG→CCA ytfR → putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,451,4830GA17.4% 55.3 / 5.4 23P142P (CCG→CCAytfRputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (7/12);  new base A (2/2);  total (9/14)
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

ATAATCTATTTATAGGCCGCGAACCCAAACGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTC  >  NC_000913/4451365‑4451625
                                                                                                                      |                                                                                                                                              
aTAATCTATTTATAGGCCGCGAACCCAAACGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCaaa                                                                                                                 <  1:97438/150‑1 (MQ=255)
aTAATCTATTTATAGGCCGCGAACCCAAACGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCaaa                                                                                                                 <  1:97432/150‑1 (MQ=255)
    tCTATTTATAGGCCGCGAACCCAAACGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAAtcgt                                                                                                             >  2:170160/1‑150 (MQ=255)
                   cgAACCCAAACGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGc                                                                                              <  1:168706/150‑1 (MQ=255)
                           aaCGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCACTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGAtc                                                                                      >  1:179210/1‑150 (MQ=255)
                           aaCGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCACTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGAtc                                                                                      >  1:179211/1‑150 (MQ=255)
                                     ccTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCaa                                                                            >  1:155526/1‑150 (MQ=255)
                                     ccTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCaa                                                                            >  1:88292/1‑150 (MQ=255)
                                      cTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCaaa                                                                           <  1:116583/150‑1 (MQ=255)
                                      cTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCaaa                                                                           <  1:240856/150‑1 (MQ=255)
                                                 cAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCt                                                                <  1:115616/150‑1 (MQ=255)
                                                    agagATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGa                                                             >  1:29769/1‑150 (MQ=255)
                                                        aTGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGaa                                                         <  2:232717/150‑1 (MQ=255)
                                                               aGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCg                                                  <  1:11955/150‑1 (MQ=255)
                                                                    gcCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGt                                             <  1:217844/150‑1 (MQ=255)
                                                                                     aTCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGACGt                            >  1:232013/1‑150 (MQ=255)
                                                                                     aTCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGt                            >  1:232019/1‑150 (MQ=255)
                                                                                        ttACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGa                         >  1:223866/1‑150 (MQ=255)
                                                                                          aCGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGt                       <  1:109002/150‑1 (MQ=255)
                                                                                            ggTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTa                     <  1:113927/150‑1 (MQ=255)
                                                                                                    ccTCGACGTGCGCGAACCACTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGa             <  2:179211/150‑1 (MQ=255)
                                                                                                    ccTCGACGTGCGCGAACCACTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGa             <  2:179210/150‑1 (MQ=255)
                                                                                                               gcgAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTc  <  2:33852/150‑1 (MQ=255)
                                                                                                                      |                                                                                                                                              
ATAATCTATTTATAGGCCGCGAACCCAAACGCTTCGGCCTTCTACGCCGCAAAGAGATGGAAAAGCGCGCCACCGAACTGATGGCATCTTACGGTTTCTCCCTCGACGTGCGCGAACCGCTCAACCGCTTTTCAGTCGCGATGCAGCAAATCGTCGCTATTTGTCGGGCTATCGATCTCTCTGCCAAAGTGCTGATCCTCGATGAACCCACCGCCAGTCTCGACACCCAGGAAGTGGAGTTACTGTTTGACCTGATGCGTC  >  NC_000913/4451365‑4451625

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