Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,305,651 G→A 8.3% Q150* (CAG→TAG)  PP_3777 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,305,6510GA 8.3% 136.1 / 2.4 48Q150* (CAG→TAG) PP_3777hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base A (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.10e-01

ACCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCA  >  NC_002947/4305502‑4305794
                                                                                                                                                     |                                                                                                                                               
acCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTg                                                                                                                                                 >  2:737814/1‑150 (MQ=255)
acCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTg                                                                                                                                                 <  1:698131/150‑1 (MQ=255)
acCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAACGACAGCCCTTGGGCAGCCTCGATGCGCTg                                                                                                                                                 >  1:211116/1‑150 (MQ=255)
  cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                               >  1:1028908/1‑150 (MQ=255)
  cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                               <  2:1023579/150‑1 (MQ=255)
  cGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTa                                                                                                                                               <  1:492734/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                                                                                                                                     >  2:573592/1‑124 (MQ=255)
                                      gCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAg                                                                                                                                     <  1:573592/124‑1 (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                                                                            <  1:729280/144‑1 (MQ=255)
                                                                           gAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATcc                                                                            >  2:729280/1‑144 (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                                                                                        >  1:641331/1‑68 (MQ=255)
                                                                                                                                           tCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGAt                                                                                        <  2:641331/68‑1 (MQ=255)
                                                                                                                                           tCGATGCACTACAGCGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTg      >  2:350813/1‑150 (MQ=255)
                                                                                                                                               tGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCa  <  2:343108/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
ACCGGCAATGCCGACCAACGAATTCAGCCCGATGATCAGCTTGCCTCCGGGCTTCAGGTAGCTCGGCACCCGTTCGAAAAACCTGCGCAAGGTCCGCTCGCCATAGAACATCGCTTCCAGCGACAGCCCTTGGGCAGCCTCGATGCGCTGTAGTGCATCGAGGAAGCCATCCGAGGGCGTCTGCCACCCTGGCGGATTGAATAGGATCACATCGTATCCGTCCTCTGCCGCAGGCATTGCGTCGAACAGATCTCCCTGTACGAATTCGATACAGGATTCGCCCAGGTTGACCA  >  NC_002947/4305502‑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: