Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,146,373 G→A 8.0% G326G (GGC→GGT)  PP_5392 ← WD40/YVTN repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,146,3730GA 8.0% 140.3 / 2.3 50G326G (GGC→GGT) PP_5392WD40/YVTN repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base A (2/2);  total (22/28)
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.86e-01

CTTGGTCTTGGCGCCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCGGTGG  >  NC_002947/6146228‑6146506
                                                                                                                                                 |                                                                                                                                     
cTTGGTCTTGGCGCCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGcc                                                                                                                                     >  7:243264/1‑148 (MQ=255)
 ttGGTCTTGGCGCCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCaa                                                                                                                                   >  4:118218/1‑149 (MQ=255)
         ggCGCCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTAt                                                                                                                             <  2:310248/147‑1 (MQ=255)
         ggCGCCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATg                                                                                                                            <  3:118218/148‑1 (MQ=255)
           cgcCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGt                                                                                                                         <  2:54403/149‑1 (MQ=255)
            gcCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTa                                                                                                                        <  1:155280/149‑1 (MQ=255)
             ccAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCAGTTTTCGTCTCGTGCTCAGCTTGCACATCACCGGGCTATTGGCCACTCGGCGGCTCGCAGCTGATCCCGTTTGGCCCGGCGCCAACCTGATGCGTAg                                                                                                                       >  2:128808/1‑149 (MQ=255)
                      aCCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTTGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGAc                                                                                                                <  3:107776/147‑1 (MQ=255)
                            cACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCACCGACCTTACGgc                                                                                                        >  4:181588/1‑149 (MQ=255)
                              cGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACATTACTgct                                                                                                       >  6:62606/1‑148 (MQ=255)
                               gACAAACAAATGACTGCGCCGGAGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATAAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTATCAACGACCTTACGgctgc                                                                                                     <  5:62606/149‑1 (MQ=255)
                                   aacaaaTGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGgct                                                                                                       >  7:294631/1‑143 (MQ=255)
                                   aacaaaTGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGgct                                                                                                       <  8:294631/143‑1 (MQ=255)
                                   aacaaaTGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCg                                                                                                  >  3:287866/1‑148 (MQ=255)
                                      aaaTGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATc                                                                                              >  6:238924/1‑149 (MQ=255)
                                      aaaTGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATc                                                                                              >  1:32236/1‑149 (MQ=255)
                                          gACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATg                                                                                          >  2:99616/1‑149 (MQ=255)
                                               cgcCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACa                                                                                       <  3:105814/147‑1 (MQ=255)
                                               cgcCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAg                                                                                      <  8:243264/148‑1 (MQ=255)
                                                          aTGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGTCTg                                                                          <  5:124604/149‑1 (MQ=255)
                                                             cTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACa                                                                                       <  6:133773/133‑1 (MQ=255)
                                                             cTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACa                                                                                       >  5:133773/1‑133 (MQ=255)
                                                              tCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCtt                                                                       >  7:34808/1‑148 (MQ=255)
                                                              tCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCtt                                                                       >  7:302992/1‑148 (MQ=255)
                                                                aTTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCaa                                                                    <  3:250554/149‑1 (MQ=255)
                                                                 tttttCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCaa                                                                    >  8:192502/1‑148 (MQ=255)
                                                                    ttCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAatat                                                                 <  8:34808/148‑1 (MQ=255)
                                                                       gTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCt                                                             <  5:246773/149‑1 (MQ=255)
                                                                       gTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCt                                                             <  7:192502/149‑1 (MQ=255)
                                                                                  aGGTTGCACATCACCCCGCGATTGTGCACGCGGCGGCTCTCACCTGAGCCCGTTTGGCCCGGCGCCAACCTTATCGGTAGCAACGACGTTACGGCTGCTCGTATCTATGACAGAAACGGTGCCTCCTTCAATATTGCTGACAAACACAt                                                  <  4:99994/149‑1 (MQ=255)
                                                                                  aGGTTGCACATCAACGGGCTATTGGGCACTCCGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAacac                                                    <  8:276220/147‑1 (MQ=255)
                                                                                    gTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATaa                                                <  2:175792/149‑1 (MQ=255)
                                                                                      tGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGCCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAgc                                              <  8:302992/149‑1 (MQ=255)
                                                                                           aTCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGt                                         >  3:162362/1‑149 (MQ=255)
                                                                                              aaCGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAgc                                              <  1:138625/141‑1 (MQ=255)
                                                                                              aaCGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAgc                                              >  2:138625/1‑141 (MQ=255)
                                                                                                           gCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCcac                         >  2:103466/1‑149 (MQ=255)
                                                                                                               tCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCACCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATaa                                                <  5:65225/122‑1 (MQ=255)
                                                                                                               tCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCACCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATaa                                                <  5:60587/122‑1 (MQ=255)
                                                                                                               tCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCACCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATaa                                                >  6:65225/1‑122 (MQ=255)
                                                                                                               tCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCACCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATaa                                                >  6:60587/1‑122 (MQ=255)
                                                                                                                 ggcggcTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATg                   <  4:287866/149‑1 (MQ=255)
                                                                                                                      cTCGAAGCTTATCCCGTTTGGCCCGGCGCAAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCAc              <  3:181588/149‑1 (MQ=255)
                                                                                                                          aaGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTACTTCAATATTGCTGACAAACACATAAgcgc                                            >  7:129525/1‑115 (MQ=255)
                                                                                                                          aaGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTACTTCAATATTGCTGACAAACACATAAgcgc                                            <  8:129525/115‑1 (MQ=255)
                                                                                                                               gATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCgg     >  5:234902/1‑149 (MQ=255)
                                                                                                                                aTCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCggt    >  3:244331/1‑149 (MQ=255)
                                                                                                                                 tCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCggtg   <  5:238924/149‑1 (MQ=255)
                                                                                                                                 tCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACCACCTTACGGCTGCTCGTATCGATGACACATACGGTGCCTGCTTCATTATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCggt    <  1:99616/148‑1 (MQ=255)
                                                                                                                                   ccGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCggtgg  <  2:32236/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                     
CTTGGTCTTGGCGCCAAGCCTCACCTGGCACGACAAACAAATGACTGCGCCGGCGTGAATGCTCATTTTTCGTCTCGTGCTCAGGTTGCACATCAACGGGCTATTGGGCACTCGGCGGCTCGAAGCTGATCCCGTTTGGCCCGGCGCCAACCTTATGGGTAGCAACGACCTTACGGCTGCTCGTATCGATGACAGATACGGTGCCTGCTTCAATATTGCTGACAAACACATAAGCGCCGTCGCTGCTGATGGCCACACCATGGGCACCCTTGCCGGTGG  >  NC_002947/6146228‑6146506

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