Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,532,626 T→C 12.9% M105T (ATG→ACG) ‡ PP_4020 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,532,6260TC12.9% 73.8 / 4.0 31M105T (ATG→ACG) ‡PP_4020oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (14/13);  new base C (2/2);  total (16/15)
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.61e-01

CGTAACGCTGCTGGAAGCCAGCCGGGTTGCCTGGGCGGCTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAAGAGGCCCTGCATGATGC  >  NC_002947/4532483‑4532758
                                                                                                                                               |                                                                                                                                    
cgTAACGCTGCTGGAAGCCAGCCGGGTTGCCTGGGCGGCTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTg                                                                                                                                <  2:200019/150‑1 (MQ=255)
        tgctgGAAGCCAGCCGGGTTGCCTGGGCGGCTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGccgccg                                                                                                                        >  2:173949/1‑150 (MQ=255)
                             ccTGGGCGGCTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGCCCAGCACATGCTggcc                                                                                                                             <  2:94423/124‑3 (MQ=255)
                             ccTGGGCGGCTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGCCCAGCACATGCTggcc                                                                                                                             >  1:94423/1‑122 (MQ=255)
                                 ggcggcTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACg                                                                                               >  2:385192/1‑150 (MQ=255)
                                      cTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACg                                                                                          <  1:173949/150‑1 (MQ=255)
                                       ttCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACgg                                                                                         >  2:451522/1‑150 (MQ=255)
                                       ttCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACgg                                                                                         <  1:265747/150‑1 (MQ=255)
                                                 aaCGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATc                                                                               >  2:243626/1‑150 (MQ=255)
                                                 aaCGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGCCCAGCACATGCTggccca                                                                                                                           >  1:446648/1‑102 (MQ=255)
                                                 aaCGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGCCCAGCACATGCTGGccca                                                                                                                           <  2:446648/106‑2 (MQ=255)
                                                   cGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGa                                                                             <  2:20111/150‑1 (MQ=255)
                                                                              gTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGc                                                  >  2:397813/1‑150 (MQ=255)
                                                                                     aTGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTg                                           >  2:348202/1‑150 (MQ=255)
                                                                                       gccgccGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGcc                                         >  1:88720/1‑150 (MQ=255)
                                                                                                   gCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGaa                             >  1:456622/1‑150 (MQ=255)
                                                                                                       gCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCtgt                                                         <  1:606997/118‑1 (MQ=255)
                                                                                                       gCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCtgt                                                         >  2:606997/1‑118 (MQ=255)
                                                                                                         ccTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATc                                                                               >  1:544196/1‑94 (MQ=255)
                                                                                                         ccTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATc                                                                               <  2:544196/94‑1 (MQ=255)
                                                                                                          cTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTc                      <  1:348202/150‑1 (MQ=255)
                                                                                                             ggCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGaa                   <  1:451522/150‑1 (MQ=255)
                                                                                                               cATCGAACGCACCGAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAaga                 <  1:29381/150‑1 (MQ=255)
                                                                                                                 tCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAAGAgg               <  1:385192/150‑1 (MQ=255)
                                                                                                                 tCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCAACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTCTTGCCGCGCCGCGCGAAAATGCTCGAagagc               >  2:276217/1‑149 (MQ=255)
                                                                                                                                        gCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAAGAGGCCCTGCATGATGc  >  1:60145/1‑140 (MQ=255)
                                                                                                                                        gCCAGCATGGGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAAGAGGCCCTGCATGATGc  <  2:60145/140‑1 (MQ=255)
                                                                                                                                          cAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAag                  <  2:521908/122‑1 (MQ=255)
                                                                                                                                          cAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAag                  >  1:521908/1‑122 (MQ=255)
                                                                                                                                            gCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTAcc                                                                        >  2:471489/1‑66 (MQ=255)
                                                                                                                                            gCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTAcc                                                                        <  1:471489/66‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                    
CGTAACGCTGCTGGAAGCCAGCCGGGTTGCCTGGGCGGCTTCCGGACGCAACGGTGGCCAGGCGCTGTTGGGCTGGTCGTGCGACATGCCGCCGCTGGAGCAGGCCCTTGGCATCGAACGCACCCAACGGCTATGGGCCAGCATGTGCTGGGCCGCCGAAGAGATGCGCGAACTGCCCCAACGCCACGGCTTCGACATCGACTACCGGCAGGGCAGCCTGTGGACTGCGGTATTGCCGCGCCGGGTGAAAATGCTCGAAGAGGCCCTGCATGATGC  >  NC_002947/4532483‑4532758

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