Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,268,815 A→G 10.2% A451A (GCA→GCG)  PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,268,8150AG10.2% 104.5 / 3.4 39A451A (GCA→GCG) PP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base A (16/19);  new base G (2/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

CTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGC  >  NC_002947/5268674‑5268964
                                                                                                                                             |                                                                                                                                                     
cTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTACAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGa                                                                                                                                               >  2:312357/1‑150 (MQ=255)
       caccacTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCAcc                                                                                                                                        <  2:331220/150‑1 (MQ=255)
           acTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATcc                                                                                                                                    >  2:324032/1‑150 (MQ=255)
            cTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCa                                                                                                                                   <  2:72977/150‑1 (MQ=255)
                     cttgTTCGAGACCCTCCTCCACTCCGGGATTAGCTCGGCCGCGATCATGGCCATACTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGa                                                                                                                          <  2:94317/147‑1 (MQ=255)
                                        cACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGAt                                                                                                       <  2:235345/150‑1 (MQ=255)
                                         aCTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATg                                                                                                      >  1:194893/1‑150 (MQ=255)
                                                    aGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTc                                                                                           >  2:314701/1‑150 (MQ=255)
                                                     gCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCa                                                                                          >  2:87907/1‑150 (MQ=255)
                                                      cTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCaa                                                                                         <  2:331407/150‑1 (MQ=255)
                                                       tCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAg                                                                                        >  2:51398/1‑150 (MQ=255)
                                                           cggcgATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATg                                                                                    <  2:139287/150‑1 (MQ=255)
                                                           cggcgATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACGCGGAAAACCGGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGTGACTACTTCAAGGATg                                                                                    <  2:215788/150‑1 (MQ=255)
                                                                  cATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGca                                                                             <  2:169183/150‑2 (MQ=255)
                                                                        cATCCTGCTGAACCTGATCTTCCACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGa                                                                       >  1:383819/1‑150 (MQ=255)
                                                                               ctgAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAg                                                                >  1:164069/1‑150 (MQ=255)
                                                                                tgAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAgg                                                               <  1:324032/150‑1 (MQ=255)
                                                                                               aCCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCGGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCaa                                                                                         <  2:54071/109‑1 (MQ=255)
                                                                                               aCCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCGGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCaa                                                                                         >  1:54071/1‑109 (MQ=255)
                                                                                               aCCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCGGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCaa                                                                                         >  1:300410/1‑109 (MQ=255)
                                                                                               aCCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCGGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCaa                                                                                         <  2:300410/109‑1 (MQ=255)
                                                                                               aCCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGa                                                <  1:314701/150‑1 (MQ=255)
                                                                                               aCCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGa                                                >  1:197998/1‑150 (MQ=255)
                                                                                                 cATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAg                                                                                                                                  >  2:221415/1‑66 (MQ=255)
                                                                                                 cATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAg                                                                                                                                  <  1:221415/66‑1 (MQ=255)
                                                                                                  aTTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGc                                             <  1:51398/150‑1 (MQ=255)
                                                                                                    ttCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTg                                           <  2:197998/150‑1 (MQ=255)
                                                                                                          accGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTg                                     >  2:454658/1‑150 (MQ=255)
                                                                                                              gCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcg                                                                                                                                            >  1:281908/1‑43 (MQ=255)
                                                                                                              gCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAgcg                                                                                                                                            <  2:281908/43‑1 (MQ=255)
                                                                                                              gCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAgg                                                               <  1:207662/120‑1 (MQ=255)
                                                                                                              gCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAgg                                                               >  2:207662/1‑120 (MQ=255)
                                                                                                                  cTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCa                             <  2:93470/150‑1 (MQ=255)
                                                                                                                   tCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCAc                            >  1:417546/1‑150 (MQ=255)
                                                                                                                    cGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACg                           <  2:194893/150‑1 (MQ=255)
                                                                                                                         aaCCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAAt                      <  1:87907/150‑1 (MQ=255)
                                                                                                                              gtcagtGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCACTTATGGTCCGAGGGTAACGAAGTGACGATGCTGAAACTGGACGAGCACGGCAATGAAAc                 >  1:5126/1‑150 (MQ=255)
                                                                                                                               tcagtGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACg                <  2:164069/150‑1 (MQ=255)
                                                                                                                                  gtgtTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTc             >  2:107860/1‑150 (MQ=255)
                                                                                                                                             aGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGc  >  2:140033/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CTTCTACCACCACTTCCCGAACTGGTTCGAGACCATCTTCCACTCCGGCATCAGCTCGGCGGCGATCATGGCCATCCTGCTGAACCTGATCTTCAACCATTTCACCACCGGCAACTCGGAAAACCAGTCAGTGTTCGCTGCAGCCTACGAGCGCACCATCCAGTATTCGGACATTTCCGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGC  >  NC_002947/5268674‑5268964

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