Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,219,155 A→G 16.0% D229G (GAC→GGC) ‡ PP_1064 → alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,219,1550AG16.0% 68.8 / 5.3 31D229G (GAC→GGC) ‡PP_1064alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base A (13/13);  new base G (2/3);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

ACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTAC  >  NC_002947/1219008‑1219303
                                                                                                                                                   |                                                                                                                                                    
aCGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGAcc                                                                                                                                                    >  1:395870/1‑150 (MQ=255)
         gCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTgg                                                                                                                                           >  1:510946/1‑150 (MQ=255)
          cGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCggg                                                                                                                                          <  1:327063/150‑2 (MQ=255)
          cGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGt                                                                                                                                          <  2:510946/150‑1 (MQ=255)
                        aTGTATTGGGTGTGAGCATGGGCGTCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGccc                                                                                                                            >  2:24587/1‑150 (MQ=255)
                                        cATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGgcgccggg                                                                                                                                          >  1:132091/1‑116 (MQ=255)
                                        cATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCggg                                                                                                                                          <  2:132091/120‑2 (MQ=255)
                                               ggcATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGaa                                                                                                     <  2:203240/150‑1 (MQ=255)
                                                  aTGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTg                                                                                                  <  1:429402/150‑1 (MQ=255)
                                                        gCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                                                                            <  1:498721/150‑1 (MQ=255)
                                                        gCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                                                                            <  1:253356/150‑1 (MQ=255)
                                                                  tGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                                  >  2:163955/1‑150 (MQ=255)
                                                                   ggcggcGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGc                                                                                 <  1:163955/150‑1 (MQ=255)
                                                                    gcggcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGcc                                                                                >  1:255301/1‑150 (MQ=255)
                                                                         gATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGAGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGACCGGCGCAGCGCACCGAACCGTCAGGTGGCAAGCGAACGACAGGCCTACct                                                                           >  2:261161/1‑150 (MQ=255)
                                                                               gcCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGTCGGCctgctgct                                                                     <  2:381423/150‑3 (MQ=255)
                                                                                cccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcg                                                                    >  2:126181/1‑150 (MQ=255)
                                                                                                    gACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCt                                                                <  1:432093/134‑1 (MQ=255)
                                                                                                    gACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCt                                                                >  2:432093/1‑134 (MQ=255)
                                                                                                                 tCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcgc                                   >  2:447252/1‑150 (MQ=255)
                                                                                                                 tCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcgc                                   >  2:318191/1‑150 (MQ=255)
                                                                                                                     gTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAgg                               >  2:400005/1‑150 (MQ=255)
                                                                                                                      ttCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGgt                              <  2:126699/150‑1 (MQ=255)
                                                                                                                           gggCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCGGCCCTGAGGTACGGGACGACCCCGAGGTGGtgc                         >  2:195380/1‑150 (MQ=255)
                                                                                                                                     gCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGTAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggc               <  1:64405/150‑1 (MQ=255)
                                                                                                                                         ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCt                                                     <  2:115267/108‑1 (MQ=255)
                                                                                                                                         ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCt                                                     >  1:115267/1‑108 (MQ=255)
                                                                                                                                          cagcgcCGGGTCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGag                                                                                         >  2:401199/3‑71 (MQ=255)
                                                                                                                                          cagcgcCGGGTCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGag                                                                                         <  1:401199/69‑1 (MQ=255)
                                                                                                                                                cGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCt    <  1:400005/150‑1 (MQ=255)
                                                                                                                                                  gACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGACCTTCCGCTGCACCAGGCGGCGCAGGCCTAc  <  1:323948/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTAC  >  NC_002947/1219008‑1219303

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