Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,178,277 G→T 42.1% 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,2770GT42.1% 13.2 / 17.4 22E27* (GAA→TAA) yqiCBMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base G (6/6);  new base T (4/5);  total (11/11)
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

ATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTT  >  CP009273/3178129‑3178405
                                                                                                                                                    |                                                                                                                                
atggcgaagaaacggtctaaataggctgattcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGa                                                                                                                                 >  1:15850‑M2/126‑150 (MQ=255)
      gTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAg                                                                                                                               >  2:232811/1‑146 (MQ=255)
          aacggtctaaataggctgattcaaagcatgtacggaagaaaaaatcggctcaaagatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGgag                                                                                                                       >  1:145736‑M2/116‑150 (MQ=255)
            cggtctaaataggcagattcaaggcatttacgggagcatcaatagggtcaaccatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGaa                                                                                                                     <  1:115206‑M2/37‑1 (MQ=255)
                   aataggctgattcaaggcatttacgggagaaaaaatcgactcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGaaaaa                                                                                                                  <  2:80597‑M2/40‑1 (MQ=255)
                        gctgattcaaggcattgacgggagaaaaaatcggcccaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATcc                                                                                                             <  2:89178‑M2/45‑1 (MQ=255)
                        gccgaatcaaggcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCa                                                                                                         <  1:78954‑M2/49‑1 (MQ=255)
                               caaggcatttacgggagaaaaaatcagctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGATAAAAATCCGCCAaacccgt                                                                                                  >  1:195446‑M2/95‑148 (MQ=255)
                                        aacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCCGGGAGTTCGGGGGAGATGTGGAGACAAAAATCCGCCAAACCCTACAAGCGCAg                                                                                         >  1:184651‑M2/86‑150 (MQ=255)
                                         tCAACACCTATTAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGc                                                                                        <  1:223154/150‑1 (MQ=255)
                                                     atcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCg                                                                            <  1:84019‑M2/78‑1 (MQ=255)
                                                        aCAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGAcc                                                                         <  1:82285/150‑1 (MQ=255)
                                                         cAGGGAAACTACAATGATTGACCCGTAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCAAAACCCTACAAGCGCAGCTGACGCGCCTCGACCt                                                                        <  1:167719/150‑1 (MQ=255)
                                                             gAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCt                                                                        <  2:67231/146‑1 (MQ=255)
                                                                         aTTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAACAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAGGCCGTGa                                                            >  2:69820/1‑146 (MQ=255)
                                                                                   aaaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGAAATCCGGAAGTTCGGGAAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTcc                                              >  1:117313/1‑150 (MQ=255)
                                                                                                         ccAGGTTCTCGAAGCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTAc                            <  2:174557/146‑1 (MQ=255)
                                                                                                           cgcttattcgccccttccCTAAAGGAATGAGGGAGTTCGGGGAAGATGAGCAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGTCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACgg                          <  2:127249‑M2/128‑1 (MQ=255)
                                                                                                                 ttcgcaccttccCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaa                >  1:208442‑M2/13‑150 (MQ=255)
                                                                                                                             cTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCACCTCGACCTTGTAACCCGTGAAGAGTTCGACGTCCAAACGCCAGTCCTGTTACGGACCCGGGAACACCTGgcc        >  2:98758/1‑145 (MQ=255)
                                                                                                                               aaaGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAATCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt  >  1:38238/1‑150 (MQ=255)
                                                                                                                               aaaGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAAGGCAAGTCCTGTTACGGAGTCGTGAAAAAATGGCgct      >  2:54412/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                
ATTCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTT  >  CP009273/3178129‑3178405

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: