Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,568,739 T→G 18.1% E284D (GAA→GAC PP_2252 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,568,7390TG18.1% 48.3 / 9.4 29E284D (GAA→GACPP_2252hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (8/15);  new base G (2/3);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

ACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCT  >  NC_002947/2568596‑2568888
                                                                                                                                               |                                                                                                                                                     
aCCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGATGGTCGaac                                                                                                                                                 <  1:100446/150‑3 (MQ=255)
        cggcgCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGACCACCCGCAGATGGTCGAACAACTGCgg                                                                                                                                         <  2:321653/150‑1 (MQ=255)
           cgcACTGGCCCGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGggtg                                                                                                                                      <  1:389777/150‑1 (MQ=255)
                  gCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACgg                                                                                                                               <  2:750029/150‑1 (MQ=255)
                       cAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCCTCTGCGGGTGGTCACGGTCATg                                                                                                                          <  2:348280/150‑1 (MQ=255)
                             gTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcat                                                                                                                    <  2:112967/150‑1 (MQ=255)
                              tCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGtcatc                                                                                                                   >  1:346661/1‑150 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGTTGGTCACGGTCATGGGTCATCAc                                                                                                                 >  2:412265/1‑150 (MQ=255)
                                cTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCAc                                                                                                                 >  1:168102/1‑150 (MQ=255)
                                              aCAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTg                                                                                                   >  1:342434/1‑150 (MQ=255)
                                                    aaTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGTTGGTCACGTTCATGGGTCATCACCCCGTAACCGGCTGCGGCaa                                                                                             <  1:787843/150‑1 (MQ=255)
                                                     aTAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAAt                                                                                            <  1:704560/150‑1 (MQ=255)
                                                         aCCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGaa                                                                                        >  2:647974/1‑150 (MQ=255)
                                                                      ccccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTc                                                                                                                                                     >  2:212833/1‑76 (MQ=255)
                                                                      ccccGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTc                                                                                                                                                     <  1:212833/76‑1 (MQ=255)
                                                                                        ccATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCg                                                         <  2:346661/150‑1 (MQ=255)
                                                                                         cATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGa                                                        <  1:785699/150‑1 (MQ=255)
                                                                                              aCTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCCGTTGTTCGACCCTCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCa                                                   <  1:235351/150‑1 (MQ=255)
                                                                                                    gCCATTCCTGTTTCAACGCCACGTTCAGGATAACCCGCAGTTGTTCCACCACCTGCGGGTGGTCAATGTCATGGGTCATCACTCCGTAGCCGGCTCCGGCAATCGAAAACATCGCCTCTGCCATCTTGTCCTCTGCCGATGCCACGGCAc                                             <  2:89622/150‑1 (MQ=255)
                                                                                                    gCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCAc                                             >  2:344088/1‑150 (MQ=255)
                                                                                                       aTTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTgg                                          <  2:843015/150‑1 (MQ=255)
                                                                                                            tGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATc                                     <  2:780630/150‑1 (MQ=255)
                                                                                                                    cGCCACGTTCAGGATCGCCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTgg                             <  1:647974/150‑1 (MQ=255)
                                                                                                                       cACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                          <  2:168102/150‑1 (MQ=255)
                                                                                                                       cACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGc                          >  2:374266/1‑150 (MQ=255)
                                                                                                                                    aCCCGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg                <  1:522490/147‑1 (MQ=255)
                                                                                                                                    aCCCGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGgcg                >  2:522490/1‑147 (MQ=255)
                                                                                                                                     cccGCAGTTGCTTGACCCTCTGCGTTTGATCAAGGTCACGGGTCATCCCCCCGTAACCGCCTGCGGCAATCGAAATCATCCCCCCTGCCCCCTCGCCCCCTCCCCCTGCCACGGCACTCGCGATCCCCCGTGGTCCCAGGGTGTCGCCaa            <  2:32013/150‑1 (MQ=255)
                                                                                                                                      ccGCAGATGGTCGAACAACTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGTCAAg           >  1:591066/1‑150 (MQ=255)
                                                                                                                                            ttgttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                       <  2:586611/100‑1 (MQ=255)
                                                                                                                                            ttgttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGAt                                                       >  1:586611/1‑100 (MQ=255)
                                                                                                                                              gttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTCGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGc   >  2:599371/1‑150 (MQ=255)
                                                                                                                                               ttCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCGCCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCATCTACACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCt  >  1:692427/1‑150 (MQ=255)
                                                                                                                                               ttCGACCATCTGCGGGTGGTCACGGTCATGCGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCGTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCt  <  2:518919/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ACCGGCCGCGGCGCACTGGCCAGCAAGCTGTCCTCAACCTGCCCCGACAGGTAATAGACCCCCGCCATGGCCCCGCTCTGGCGGTTCTCCATCAACTCCTGCCATTCCTGGTTCAACGCCACGTTCAGGATCACCCGCAGTTGTTCGACCATCTGCGGGTGGTCACGGTCATGGGTCATCACCCCGTAACCGGCTGCGGCAATCGAAATCATCGCCCCTGCCACCTTGCCCACTGCCGATGCCACGGCACTGGCGATCCCGCGTGGGGCCAGGGTGGCGCCAAGCTTTTCGCT  >  NC_002947/2568596‑2568888

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