Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,178,277 G→T 60.8% E27* (GAA→TAA)  yqiC → BMFP family putative fusogenic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,178,2770GT60.8% 1.7 / 16.3 18E27* (GAA→TAA) yqiCBMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base T (7/4);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

AATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAA  >  CP009273/3178128‑3178408
                                                                                                                                                     |                                                                                                                                   
aatggcgaagaaacggcttaaataggctgattgaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCgggg                                                                                                                                     <  1:119526‑M2/24‑1 (MQ=255)
                                 cAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGGAAAAGATCCGCCAAACCCTAc                                                                                                    >  2:24132/1‑150 (MQ=255)
                                    tGGTACTCAACACCTATCGCGACAGGGAAAGTACAATGATTGCCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAg                                                                                                 <  1:175857/150‑1 (MQ=255)
                                        tgacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCa                                                                                             <  2:97601‑M2/64‑1 (MQ=255)
                                          acgggcgaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGc                                                                                           <  2:94882‑M2/66‑1 (MQ=255)
                                                  tATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAATTCGGGTGAGATGTGGAGAAAAAAATCCGCCAAACGCTACAAGCGGAGCTGACgcgc                                                                                   >  1:179011/1‑150 (MQ=255)
                                                                         aatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAaga                                                            <  1:47811‑M2/97‑1 (MQ=255)
                                                                          atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAATGAATCAGGGAGTTCGGGGAAGCTGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTTGACCGTGTCAGCCGTGAagag                                                           >  2:274133‑M2/53‑150 (MQ=255)
                                                                                              gCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAGGCCGTGAAGAGTTCGACGTCCAAACGCAAGt                                       >  2:15271/1‑150 (MQ=255)
                                                                                              gCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGt                                       >  2:322975/1‑150 (MQ=255)
                                                                                               cAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGGGGTCGACGTCCAAACGCAAGTc                                      >  2:186209/1‑150 (MQ=255)
                                                                                               cAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTc                                      >  2:261644/1‑150 (MQ=255)
                                                                                               cAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTc                                      >  2:313802/1‑150 (MQ=255)
                                                                                                gaagaatttccccgcttattcgcaccgtccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTcc                                     >  2:214124‑M2/31‑150 (MQ=255)
                                                                                                       tcgcCCGGTTCACCAATCCATGCCTAAAGGAATCAGGGGGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCGACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACg                              <  1:34666/150‑1 (MQ=255)
                                                                                                        cgcCAGGTTCACGAATCCATGCCTAAAGGTATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTCACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACgg                             <  1:264085/150‑1 (MQ=255)
                                                                                                                               tAAAGGAAGCAGGGAGTTCGGGTAAGATGGGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGCCGCGCCTCAACCTTGTAAGCCGTGAAGACTTCGACGTCCAACCGCAAGGGCTGTTACGGACGTTGGAAACACTGGCgctgct      <  1:339703/150‑1 (MQ=255)
                                                                                                                                  aGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa   >  2:281890/1‑150 (MQ=255)
                                                                                                                                   ggAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTCCGGACTCGTGAAAAACTGGCGCTGCTTGaa  >  1:29803/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                   
AATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAA  >  CP009273/3178128‑3178408

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