Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,331 C→T 29.4% intergenic (+42/‑80) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3310CT29.4% 22.1 / 10.0 17intergenic (+42/‑80)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base T (2/3);  total (8/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 = 1.00e+00

TTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGT  >  NC_000913/4612263‑4612453
                                                                    |                                                                                                                          
ttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                                                       >  4:998669/1‑74 (MQ=255)
ttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACt                                                                                                                       <  3:998669/74‑1 (MQ=255)
                    aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                              >  1:374439/1‑79 (MQ=255)
                    aCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGaataat                                                                                              <  2:374439/79‑1 (MQ=255)
                                                         aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                              >  1:545504/1‑90 (MQ=255)
                                                         aGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCaga                                              <  2:545504/90‑1 (MQ=255)
                                                          gCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                    >  3:30047/1‑67 (MQ=38)
                                                          gCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTc                                                                    <  4:30047/67‑1 (MQ=38)
                                                           ccAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTa                                                                                   <  3:611400/51‑1 (MQ=25)
                                                               tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                   >  3:301529/1‑111 (MQ=255)
                                                               tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                   <  4:301529/111‑1 (MQ=255)
                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                          >  3:632092/1‑55 (MQ=255)
                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                          <  4:632092/55‑1 (MQ=255)
                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt  <  2:970037/127‑1 (MQ=255)
                                                                ccAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGt  >  1:970037/1‑127 (MQ=255)
                                                                  aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                               >  3:385319/1‑80 (MQ=255)
                                                                  aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                               <  4:385319/80‑1 (MQ=255)
                                                                    cGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                        >  2:808140/1‑69 (MQ=255)
                                                                    cGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                        <  1:808140/69‑1 (MQ=255)
                                                                    |                                                                                                                          
TTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGT  >  NC_000913/4612263‑4612453

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