Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,305,652 T→C 8.3% L149L (CTA→CTG PP_3777 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,305,6520TC 8.3% 139.3 / 2.4 48L149L (CTA→CTGPP_3777hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base C (2/2);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCA  >  NC_002947/4305504‑4305794
                                                                                                                                                    |                                                                                                                                              
cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                               >  1:1028908/1‑150 (MQ=255)
cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                               <  1:492734/150‑1 (MQ=255)
cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                               <  2:1023579/150‑1 (MQ=255)
       gCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATc                                                                                                                                        <  2:713429/150‑1 (MQ=255)
         cGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGa                                                                                                                                      >  1:343108/1‑150 (MQ=255)
         cGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGa                                                                                                                                      <  2:33789/150‑1 (MQ=255)
          gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAg                                                                                                                                     <  1:206463/150‑1 (MQ=255)
          gACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCACTACAGCGCATCGAg                                                                                                                                     <  1:620530/150‑1 (MQ=255)
             cAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCGAGAGTCGCTCGCCATAGAACATCACTTCAGGCAAACATGCTTGGGCAGCCTCGATGGGCTGTAGTGCATCGAGGaa                                                                                                                                  >  1:1172972/1‑150 (MQ=255)
              aaCGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAg                                                                                                                                 >  2:490174/1‑150 (MQ=255)
              aaCGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAg                                                                                                                                 <  2:981917/150‑1 (MQ=255)
               aCGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGc                                                                                                                                <  1:490174/150‑1 (MQ=255)
                  aaTTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCAt                                                                                                                             <  1:737814/150‑1 (MQ=255)
                        gCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAgg                                                                                                                       >  1:695241/1‑150 (MQ=255)
                          ccGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGc                                                                                                                     >  1:1141023/1‑150 (MQ=255)
                                    gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAg                                                                                                                                     <  1:573592/124‑1 (MQ=255)
                                    gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAg                                                                                                                                     >  2:573592/1‑124 (MQ=255)
                                    gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTccgc                                                                                                               >  1:101530/1‑143 (MQ=255)
                                    gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCCGc                                                                                                               <  2:101530/146‑1 (MQ=255)
                                             cGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTg                                                                                                        <  1:527071/144‑1 (MQ=255)
                                             cGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTg                                                                                                        >  2:527071/1‑144 (MQ=255)
                                              gggCTTCAGGTAGCTCGGCACCCGTTTGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGAtt                                                                                                 <  2:620340/150‑1 (MQ=255)
                                                 cTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGaa                                                                                              <  2:1172972/150‑1 (MQ=255)
                                                                   ccGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATcc                                                                            >  2:1094682/1‑150 (MQ=255)
                                                                    cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAg                                                                                                                                     <  1:706613/92‑1 (MQ=255)
                                                                    cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAg                                                                                                                                     >  2:706613/1‑92 (MQ=255)
                                                                    cGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCg                                                                           >  1:38549/1‑150 (MQ=255)
                                                                        cGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCct                                                                       <  2:9892/150‑1 (MQ=255)
                                                                         gAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATcc                                                                            >  2:729280/1‑144 (MQ=255)
                                                                         gAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATcc                                                                            <  1:729280/144‑1 (MQ=255)
                                                                            aaaCCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGc                                                                   >  2:663637/1‑150 (MQ=255)
                                                                                          tCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGt                                                     <  2:1028908/150‑1 (MQ=255)
                                                                                           ccGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGCATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTc                                                    <  2:695241/150‑1 (MQ=255)
                                                                                              cTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTAt                                                                              <  2:32164/121‑1 (MQ=255)
                                                                                              cTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTAt                                                                              >  1:32164/1‑121 (MQ=255)
                                                                                                   cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCa                                                               <  2:453769/131‑1 (MQ=255)
                                                                                                   cATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCa                                                               >  1:453769/1‑131 (MQ=255)
                                                                                                        aaCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                     >  2:252126/1‑88 (MQ=255)
                                                                                                        aaCATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCg                                                                                                     <  1:252126/88‑1 (MQ=255)
                                                                                                                            cccTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTAc                                  <  2:1040606/135‑1 (MQ=255)
                                                                                                                            cccTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTAc                                  >  1:1040606/1‑135 (MQ=255)
                                                                                                                             ccTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCAGCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGAt                  >  2:978602/1‑150 (MQ=255)
                                                                                                                                    cttccTCGATGCACTACAGCGCATCGAGGCTGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCa           >  1:213460/4‑150 (MQ=255)
                                                                                                                                       ccTCGATGTGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGt        >  2:220004/1‑150 (MQ=255)
                                                                                                                                         tCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGAt                                                                                        <  2:641331/68‑1 (MQ=255)
                                                                                                                                         tCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGAt                                                                                        >  1:641331/1‑68 (MQ=255)
                                                                                                                                         tCGATGCACTACAGCGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTg      >  2:350813/1‑150 (MQ=255)
                                                                                                                                             tGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCa  <  2:343108/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCA  >  NC_002947/4305504‑4305794

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: