Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,256,978 C→T 8.6% R59H (CGC→CAC)  PP_2853 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,256,9780CT 8.6% 125.4 / 2.9 47R59H (CGC→CAC) PP_2853hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/20);  new base T (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

GGAATCGGCTGTCAGCGCCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCGGTGGCGGGGGCCAAGGTCGTCAAGCA  >  NC_002947/3256833‑3257120
                                                                                                                                                  |                                                                                                                                              
ggAATCGGCTGTCAGCGCCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCGCTAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGgcgc                                                                                                                                              <  1:267849/149‑1 (MQ=255)
 gAATCGGCTGTCAGCGCCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCGCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCGCGACCTGGCGca                                                                                                                                             >  5:275170/1‑148 (MQ=255)
              gcgcCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTc                                                                                                                                   <  1:943/145‑1 (MQ=255)
              gcgcCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTc                                                                                                                                   >  2:943/1‑145 (MQ=255)
                gcCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCAc                                                                                                                             >  6:36140/1‑149 (MQ=255)
                  cTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCa                                                                                                                           >  3:112710/1‑149 (MQ=255)
                        gTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGTGCACACGTTGCTCCAGCACCATCAGTg                                                                                                                     >  5:333561/1‑149 (MQ=255)
                           ggAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCTGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCAGCCTCAGTGCCt                                                                                                                  >  4:80987/1‑149 (MQ=255)
                            gAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTg                                                                                                                 >  5:89945/1‑149 (MQ=255)
                                 cccTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTAGCt                                                                                                            >  6:82162/1‑149 (MQ=255)
                                  ccTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGTGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCtg                                                                                                           >  3:202472/1‑149 (MQ=255)
                                      tAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGTGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGc                                                                                                        <  6:333561/148‑1 (MQ=255)
                                        gTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCt                                                                                                      >  4:263738/1‑148 (MQ=255)
                                         tGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTc                                                                                                     >  2:336711/1‑148 (MQ=255)
                                         tGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCa                                                                                                    <  1:64514/149‑1 (MQ=255)
                                            ccGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAg                                                                                                 <  2:187529/149‑1 (MQ=255)
                                                tGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAgcggc                                                                                             <  5:82162/149‑1 (MQ=255)
                                                   gcCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCa                                                                                                                        >  4:217563/1‑119 (MQ=255)
                                                   gcCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCa                                                                                                                        <  3:217563/119‑1 (MQ=255)
                                                    ccAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGttt                                                                                         <  2:158524/149‑1 (MQ=255)
                                                      aGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGTGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGa                                                                                       <  4:202472/149‑1 (MQ=255)
                                                      aGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGa                                                                                       <  3:263738/149‑1 (MQ=255)
                                                             aCGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCaa                                                                                 >  7:148888/1‑148 (MQ=255)
                                                             aCGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAAt                                                                                <  1:336711/149‑1 (MQ=255)
                                                              cGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATc                                                                               <  6:280320/149‑1 (MQ=255)
                                                                   tGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTc                                                                          >  7:105426/1‑149 (MQ=255)
                                                                      gCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCt                                                                        <  8:300921/148‑1 (MQ=255)
                                                                       cACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGc                                                                      >  2:317464/1‑149 (MQ=255)
                                                                               aCCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGttt                                                                                         >  4:240626/1‑122 (MQ=255)
                                                                                   ccTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTc                                                           >  4:308294/1‑148 (MQ=255)
                                                                                    cTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAAt                                                                                <  4:185609/126‑1 (MQ=255)
                                                                                    cTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAAt                                                                                >  3:185609/1‑126 (MQ=255)
                                                                                        gCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGAACACGATGCTCCAGCAAAATCAGTGACTGTTGCTGTCCCTCATAGCGGAGATTGAGTTAAATCAGTTCCTGCCTCAGCGCTTCGACtag                                                     >  7:323171/1‑147 (MQ=255)
                                                                                         ccG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTtgt                                                    <  6:89945/149‑1 (MQ=255)
                                                                                          cG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTtgtg                                                   >  5:3939/1‑149 (MQ=255)
                                                                                                    gACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTcggcgg                                          >  1:320817/1‑149 (MQ=255)
                                                                                                     aCCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCggcggc                                         >  4:175845/1‑149 (MQ=255)
                                                                                                        aaGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGGCGGCGGCTGc                                      <  8:105426/149‑1 (MQ=255)
                                                                                                               tGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAg                                <  6:3939/148‑1 (MQ=255)
                                                                                                                gggTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGt                               <  2:320817/148‑1 (MQ=255)
                                                                                                                 ggTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTGGCTCCAGCACCATCAGTGCCTGTGGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTg                              <  6:275170/148‑1 (MQ=255)
                                                                                                                  gTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCg                            <  3:175845/149‑1 (MQ=255)
                                                                                                                            agcagGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCGGTGGCggggg                  <  1:317464/149‑1 (MQ=255)
                                                                                                                              gtggggCCTCCTCGACCTGGAGCACACGTTGCTCGAGCACCATCAGTGGCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTGAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCGGTGGCGGGGGcc                <  8:323171/147‑1 (MQ=255)
                                                                                                                                       cctcGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCgg                           <  8:105109/129‑1 (MQ=255)
                                                                                                                                       cctcGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCgg                           >  7:105109/1‑129 (MQ=255)
                                                                                                                                       cctcGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCGGTGGCGGGGGCCAAGgtcgtc       <  5:36140/149‑1 (MQ=255)
                                                                                                                                            aCCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCGGTGGCGGGGGCCAAGGTCGTCAAGca  >  8:207565/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGAATCGGCTGTCAGCGCCTGGCCGTAGGAGCTCCCTGTAGTGCCCGGTGCGCCAGAAGCGACGGTCTGAGCACCTTTCACCCCCTCGGCCG‑GAGACTTGACCAAGCGTTTGGGTTGCGCGGCAGCAGGGGCCTCCTCGACCTGGCGCACACGTTGCTCCAGCACCATCAGTGCCTGTTGCTGTGCCTCATAGCGGCGTTTGAGTTCAATCAGTTCCTGCCTCAGCGCTTCGACTTGTGGGTCGGCGGCTGCGTAAAGTGCGGTGGCGGGGGCCAAGGTCGTCAAGCA  >  NC_002947/3256833‑3257120

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