Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,624,187 G→C 12.2% intergenic (‑92/+81) PP_3194 ← / ← PP_3195 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,624,1870GC12.2% 160.8 / 12.8 57intergenic (‑92/+81)PP_3194/PP_3195hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/19);  new base C (4/3);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGC  >  NC_002947/3624045‑3624332
                                                                                                                                              |                                                                                                                                                  
cAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGAt                                                                                                                                             >  2:598588/1‑150 (MQ=255)
   gCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATAcc                                                                                                                                          >  2:1027661/1‑150 (MQ=255)
     cAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATc                                                                                                                                                <  1:357365/142‑1 (MQ=255)
     cAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATc                                                                                                                                                >  2:357365/1‑142 (MQ=255)
      agggCCGTCGAACAGGGTGGAGACGTGAGCATGCAAGCCACTGATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTgg                                                                                                                                       <  2:800698/147‑1 (MQ=255)
             tCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGctctc                                                                                                                                >  2:154150/1‑150 (MQ=255)
               gAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTg                                                                                                                              >  1:418030/1‑150 (MQ=255)
                aaCAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTgg                                                                                                                             <  1:154150/150‑1 (MQ=255)
                               tgcaggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACga                                                                                                                  <  2:764593/146‑1 (MQ=255)
                               tgcaggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACga                                                                                                                  >  1:764593/1‑146 (MQ=255)
                                 caggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTgg                                                                                                                                       <  1:1033426/123‑1 (MQ=255)
                                 caggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGCCATCGATACCTgg                                                                                                                                       >  2:1033426/1‑123 (MQ=255)
                                  aggtgcaggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGgc                                                                                                           <  1:598588/150‑1 (MQ=255)
                                        aggGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGATGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAAt                                                                                                     <  2:58351/150‑1 (MQ=255)
                                                          ctctctCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAg                                                                                   <  2:418030/150‑1 (MQ=255)
                                                            ctctcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAggg                                                                                 >  1:693414/1‑150 (MQ=255)
                                                               tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCg                                                                              >  2:684860/1‑150 (MQ=255)
                                                               tcATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCg                                                                              >  1:965958/1‑150 (MQ=255)
                                                                cATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGCTCCAGGCTTGATCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCgg                                                                             >  1:812008/1‑150 (MQ=255)
                                                                cATATTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGTTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAGCGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCgg                                                                             <  2:491680/150‑1 (MQ=255)
                                                                 aTCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGc                                                                            >  2:288745/1‑150 (MQ=255)
                                                                  tCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCa                                                                           >  1:212953/1‑150 (MQ=255)
                                                                    ttGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCGATCGGCGCCGCCCAAGCCAGGGCCGGCAgc                                                                         >  1:882331/1‑150 (MQ=255)
                                                                         ggtggtGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCCGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGcgccgc                                                                                            >  1:171720/1‑126 (MQ=255)
                                                                         ggtggtGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGcgccgc                                                                                            <  2:171720/126‑1 (MQ=255)
                                                                               gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGc                                                              >  2:647980/1‑150 (MQ=255)
                                                                               gCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGc                                                              >  2:104915/1‑150 (MQ=255)
                                                                                 aaGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCg                                                            >  2:521749/1‑150 (MQ=255)
                                                                                  aggctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCgt                                                          <  1:648329/150‑1 (MQ=255)
                                                                                   ggctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCgag                                                          >  1:721070/1‑148 (MQ=255)
                                                                                    gctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCgtgt                                                        >  1:764886/1‑150 (MQ=255)
                                                                                    gctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCgtgt                                                        <  1:647980/150‑1 (MQ=255)
                                                                                    gctaggctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCACATCCGGAAGGCCTCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCgtgt                                                        >  2:595216/1‑150 (MQ=255)
                                                                                         gctGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGgccg                                                   <  2:486580/150‑1 (MQ=255)
                                                                                               gCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAgg                                             <  1:104915/150‑1 (MQ=255)
                                                                                                 gACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCAACTCGCTGCAGCGCCGATCAGGCCGTGTGGCCGCCAAgggg                                           >  1:968476/1‑149 (MQ=255)
                                                                                                 gACTGCCCGGCACCAGCAGATGCTATTGATGACCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAAAACTGCGCCAATCTGCGCCGCCCAAGCCAGGGCCGCGAGCGCGACACGTCTGATGTGGCCGCCAAGGGt                                           >  1:968112/1‑150 (MQ=255)
                                                                                                      cccGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCgg                                      >  2:750787/1‑150 (MQ=255)
                                                                                                      cccGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCgg                                      <  1:684860/150‑1 (MQ=255)
                                                                                                          gCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCGATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGAt                                  <  2:882331/150‑1 (MQ=255)
                                                                                                            aCCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATcc                                >  1:1021439/1‑150 (MQ=255)
                                                                                                             ccAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAgc                                                                         <  1:874535/109‑1 (MQ=255)
                                                                                                             ccAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAgc                                                                         >  2:874535/1‑109 (MQ=255)
                                                                                                                       cTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATg                     >  2:542644/1‑150 (MQ=255)
                                                                                                                          ttGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAg                  >  2:445217/1‑150 (MQ=255)
                                                                                                                               cgCCAGGTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTgg                                                      >  1:80884/2‑109 (MQ=255)
                                                                                                                               cgCCAGGTATCGATGCCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTgg                                                      <  2:80884/108‑1 (MQ=255)
                                                                                                                                 ccAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCgg           >  2:1136828/1‑150 (MQ=255)
                                                                                                                                  cAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGgccgc                                                  <  1:684020/110‑1 (MQ=255)
                                                                                                                                  cAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGgccgc                                                  >  2:684020/1‑110 (MQ=255)
                                                                                                                                     gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt       >  2:669156/1‑150 (MQ=255)
                                                                                                                                     gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt       <  1:288745/150‑1 (MQ=255)
                                                                                                                                     gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTAGAGGCGATCCAGGCGTCGATGTAGTGCCCGGCAtt       >  2:917654/1‑150 (MQ=255)
                                                                                                                                     gTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCATGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCAtt       >  2:823623/1‑150 (MQ=255)
                                                                                                                                        tCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGAt    <  2:764886/150‑1 (MQ=255)
                                                                                                                                         cGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATg   <  2:812008/150‑1 (MQ=255)
                                                                                                                                          gATGGCATCGATACGTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCAAAGCCAGGGCCGGCAGAGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGACCGGCATTGATGc  >  1:508424/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CAGGCCAGTGCCGTCGAACAGGGTGGAGCAGTGCAGGTGCAGGGCACTCATTGGTGAGCTCTCTCATCTTGGGGGTGGTGCAAGGCTAGGCTGTGGCGACTGCCCGGCACCAGCAGATGCTATTGATGGCCAGGTATCGATGGCATCGATACCTGGCGCTCTCTGGCGCAGAAACGAGAACGGCGCCAATCGGCGCCGCCCAAGCCAGGGCCGGCAGCGCCGATCAGGCCGTGTGGCCGCCAAGGGTTTCGGCGATCCAGTCGGCGATGTAGTGCCCGGCATTGATGC  >  NC_002947/3624045‑3624332

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