Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,397,592 A→G 9.5% R10R (CGA→CGG PP_3002 → shikimate dehydrogenase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,397,5920AG 9.5% 111.2 / 3.3 42R10R (CGA→CGGPP_3002shikimate dehydrogenase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (18/20);  new base G (2/2);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CCCGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAT  >  NC_002947/3397452‑3397738
                                                                                                                                            |                                                                                                                                                  
cccGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAAc                                                                                                                                           >  1:582922/1‑150 (MQ=255)
cccGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAAc                                                                                                                                           <  2:313603/150‑1 (MQ=255)
cccGGACCTTGGACCGCACCACGCAGTAGTGGCCGCCTTGCTTTGCTTCAGCCCGGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGTTCGGCCGACCGATCAAc                                                                                                                                           <  2:615210/150‑1 (MQ=255)
 ccGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAAcc                                                                                                                                          >  1:1031066/1‑150 (MQ=255)
 ccGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAAcc                                                                                                                                          <  2:582922/150‑1 (MQ=255)
       cTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGTCGGCCGATCactgc                                                                                                                                         <  1:1088115/145‑5 (MQ=255)
       cTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGTCGGCCGATCactgc                                                                                                                                         >  2:1088115/1‑141 (MQ=255)
          ggCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGt                                                                                                                                 <  2:172670/150‑1 (MQ=255)
            ccAGCACAACGCCGAAGTGGCCGCCGGGCTTAGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGCTCAACCACACCAAGTcc                                                                                                                               >  1:241813/1‑150 (MQ=255)
                      gCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTc                                                                                                                     >  2:1024324/1‑150 (MQ=255)
                                         ttggcttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGc                                                                                                  <  1:1024324/150‑1 (MQ=255)
                                             cttCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGc                                                                                              >  1:934933/1‑150 (MQ=255)
                                                     cTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACcagca                                                                                      >  2:2208/1‑150 (MQ=255)
                                                       gaagaaATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCagcagc                                                                                    >  1:666733/1‑150 (MQ=255)
                                                        aagaaATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCt                                                                                   <  2:666733/150‑1 (MQ=255)
                                                         agaaATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTg                                                                                  <  1:245740/150‑1 (MQ=255)
                                                           aaaTCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGc                                                                                                                           <  2:481226/107‑1 (MQ=255)
                                                           aaaTCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGc                                                                                                                           >  1:481226/1‑107 (MQ=255)
                                                                  tGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACggg                                                                         >  1:970087/1‑150 (MQ=255)
                                                                        gCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAAc                                                                                           <  1:1077905/126‑1 (MQ=255)
                                                                        gCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAAc                                                                                           >  2:1077905/1‑126 (MQ=255)
                                                                        gCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATc                                                                   <  1:2208/150‑1 (MQ=255)
                                                                        gCAGGCGGATGGCGGGGTGATTACCAAAGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGGCCCCGCTGATTCATGGCCTCTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATc                                                                   <  2:241813/150‑1 (MQ=255)
                                                                                   gCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCt                                                        <  1:376837/150‑1 (MQ=255)
                                                                                      tgttgtTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGaa                                                                               >  1:1150537/1‑124 (MQ=255)
                                                                                      tgttgtTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGaa                                                                               <  2:1150537/124‑1 (MQ=255)
                                                                                      tgttgtTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGga                                                     <  2:934933/150‑1 (MQ=255)
                                                                                            ttACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCACCCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACtt                                               <  1:530647/150‑1 (MQ=255)
                                                                                             tACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGt                                                                                                            <  1:929848/88‑1 (MQ=255)
                                                                                             tACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGt                                                                                                            >  2:929848/1‑88 (MQ=255)
                                                                                                ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACg                                                                           <  2:825502/118‑1 (MQ=255)
                                                                                                ccACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACg                                                                           >  1:825502/1‑118 (MQ=255)
                                                                                                 cACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAg                                          >  1:1111271/1‑150 (MQ=255)
                                                                                                  aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAAt                                                                              <  2:49634/113‑1 (MQ=255)
                                                                                                  aCGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAAt                                                                              >  1:49634/1‑113 (MQ=255)
                                                                                                    ggAGATGCGCGATGAGCGACCGCTACGCAGTGATCGGCCGACCGCTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGCCGACTTCGAAGccc                                       >  2:146402/1‑150 (MQ=255)
                                                                                                       gaTGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAgg                                    <  2:970087/150‑1 (MQ=255)
                                                                                                          gCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTAc                                 >  1:283678/1‑150 (MQ=255)
                                                                                                                         gCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCg                  >  2:787150/1‑150 (MQ=255)
                                                                                                                               ggttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                           <  2:137722/116‑1 (MQ=255)
                                                                                                                               ggttGATCGGTCGGCCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGaa                                           >  1:137722/4‑119 (MQ=255)
                                                                                                                                   gATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGaa        <  2:1031066/150‑1 (MQ=255)
                                                                                                                                         ccgaccgaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAt  >  1:15553/1‑150 (MQ=255)
                                                                                                                                           gaccgaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                >  2:825995/1‑54 (MQ=255)
                                                                                                                                           gaccgaTCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCa                                                                                                <  1:825995/54‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CCCGGACCTTGGCCAGCACAACGCCGAAGTGGCCGCCGGGCTTGGCTTCAGCCCTGAAGAAATCACTGCCTTGCAGGCGGATGGCGTGTTGTTTACCCACGGAGATGCACGATGAGCGACCGCTACGCAGTGATCGGCCGACCGATCAACCACACCAAGTCCCCGCTGATTCATGGCCTGTTCGCCCAGGCCAGCAACCAGCAGCTGGAATACGGGGCTATCGAAGGCTCGCTGGACGACTTCGAAGCCCAGGTACTGCAGTTTCGTAGCGAAGGCGGGAAGGGCAT  >  NC_002947/3397452‑3397738

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: