Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,223,605 C→T 22.3% Y202Y (TAC→TAT priA → primosomal protein N prime

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,223,6050CT22.3% 30.8 / 5.6 18Y202Y (TAC→TATpriAprimosomal protein N prime
Reads supporting (aligned to +/- strand):  ref base C (7/7);  new base T (2/2);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

AAAGATGATATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGA  >  CP000730/1223471‑1223730
                                                                                                                                      |                                                                                                                             
aaaGATGATATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTg                                                                                                                         <  1:256161/141‑1 (MQ=255)
     tgatATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGc                                                                                                                    <  2:117489/141‑1 (MQ=255)
                  gaaAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTgaa                                                                                                       <  2:151200/141‑1 (MQ=255)
                         cGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGaataaa                                                                                                >  1:208993/1‑140 (MQ=255)
                                             aCTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACt                                                                            <  1:11940/141‑1 (MQ=255)
                                                       ccAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAgg                                                                  <  2:58787/141‑1 (MQ=255)
                                                       ccAAGCGTACTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAgg                                                                  <  2:59206/141‑1 (MQ=255)
                                                       ccAAGCGTACTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAgg                                                                  <  2:59208/141‑1 (MQ=255)
                                                                      gtgtCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCaa                                                   <  2:26201/141‑1 (MQ=255)
                                                                      gtgtCATTGAAGGGTATCATCCTGATGAAGTATTAGATAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCaa                                                   <  2:55443/141‑1 (MQ=255)
                                                                                     atcatcCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATg                                    >  1:184862/1‑141 (MQ=255)
                                                                                        atcCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGAc                                 >  2:235804/1‑141 (MQ=255)
                                                                                                             aGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTg            >  1:254800/1‑141 (MQ=255)
                                                                                                             aGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTg            >  2:53440/1‑141 (MQ=255)
                                                                                                                  gagaAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGaaaaa       >  2:69841/1‑141 (MQ=255)
                                                                                                                   agaAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAat      >  1:199055/1‑141 (MQ=255)
                                                                                                                       aGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGa  >  1:203390/1‑141 (MQ=255)
                                                                                                                       aGTTATTAAACAATATGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGa  >  1:203335/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                             
AAAGATGATATCGTTGAAGAAAAAACGATTCTCACACAAAATATAACTAAAAAAACCAAGCGTGCTGTTCGTGTCATTGAAGGGTATCATCCTGATGAAGTATTAGCTAAGTTGGAGAAAGTTATTAAACAATACGATTTGTATGCTTACTTGTCTGAAGAACAACATAAAACAATATTTTTAACTGATATTGAGGATATGGGCTTTTCAAAATCCAGTTTAGATGGACTTATCAAAAAAGGTTATGTTGAAAAATATGA  >  CP000730/1223471‑1223730

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