Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,956,687 T→G 7.0% D1233A (GAC→GCC)  mnxG ← manganese‑oxidizing multicopper oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,956,6870TG 7.0% 160.5 / 2.4 54D1233A (GAC→GCC) mnxGmanganese‑oxidizing multicopper oxidase
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base G (2/2);  total (27/27)
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.14e-01

GCCAGCCGTGTTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGCTTGGGCAGCCCG  >  NC_002947/3956539‑3956834
                                                                                                                                                     |                                                                                                                                                   
gcCAGCCGTGTTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGt                                                                                                                                                     <  6:38001/149‑1 (MQ=255)
  cAGCCGTGTTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCg                                                                                                                                                   <  4:322884/149‑1 (MQ=255)
    gCCGTGTTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGc                                                                                                                                                 <  6:206713/149‑1 (MQ=255)
        tgtTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCTGGCCGGCGACCt                                                                                                                                               >  5:366853/1‑147 (MQ=255)
        tgtTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGCCGACGGCCTGGCCGGCGACCtca                                                                                                                                             <  7:109543/149‑3 (MQ=255)
         gtTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCa                                                                                                                                            >  3:248986/1‑149 (MQ=255)
            gATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATg                                                                                                                                          <  4:43006/148‑1 (MQ=255)
                  cggcggGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGc                                                                                                                                    <  3:115935/148‑1 (MQ=255)
                    gcggGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTg                                                                                                                                 <  4:248986/149‑1 (MQ=255)
                               gTAGTGCCGCCGGTAGGGCCGACGGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCCGGCCGTCGGCCTGCATGCCCCGCTGGCCGTCCGGGcc                                                                                                                      >  6:83078/1‑149 (MQ=255)
                                   tGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGAcc                                                                                                                  >  4:319430/1‑149 (MQ=255)
                                     ccgccgGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCtt                                                                                                                <  3:290764/149‑1 (MQ=255)
                                      cgccgGTAGGGCCGACTGCCGAGGGGATCGCCGATGGCGACAGGGTCCTCGCCGGGATGGCGCGGTCGGTACGGGTGTGCAGGGCATAGCTGAGGGCGCCGGCCAGGCCGGCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTg                                                                                                               <  5:83078/149‑1 (MQ=255)
                                         cgGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTc                                                                                                            <  4:119211/149‑1 (MQ=255)
                                          ggTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTTCAGGGCAATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGc                                                                                                                                         <  1:63542/120‑1 (MQ=255)
                                          ggTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGc                                                                                                                                         >  2:63542/1‑119 (MQ=255)
                                           gTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCgg                                                                                                          <  5:366723/149‑1 (MQ=255)
                                                gCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCa                                                                                                      >  7:118461/1‑148 (MQ=255)
                                                gCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCAt                                                                                                     <  8:118461/149‑1 (MQ=255)
                                                         cGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACAcgcg                                                                                              >  7:40666/1‑147 (MQ=255)
                                                               gATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGc                                                                                        >  6:34816/1‑147 (MQ=255)
                                                                 tCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGc                                                                                    >  7:10426/1‑149 (MQ=255)
                                                                  cGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTg                                                                                                                                 >  6:363435/1‑103 (MQ=255)
                                                                  cGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTg                                                                                                                                 <  5:363435/103‑1 (MQ=255)
                                                                                 gTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGt                                                                    >  2:62169/1‑149 (MQ=255)
                                                                                   tCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGc                                                                                                                                    <  3:51179/83‑1 (MQ=255)
                                                                                   tCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGc                                                                                                                                    >  4:51179/1‑83 (MQ=255)
                                                                                        gCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTgacga                                                              <  1:62169/148‑1 (MQ=255)
                                                                                        gCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTga                                                                 >  8:295854/1‑145 (MQ=255)
                                                                                        gCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTga                                                                 <  7:295854/145‑1 (MQ=255)
                                                                                              aTGGCGCGGTCGGTACGGGTGTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCCGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACAt                                                       <  7:226717/149‑1 (MQ=255)
                                                                                                ggCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATcc                                                     >  1:40848/1‑149 (MQ=255)
                                                                                                           tACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATcc                                          <  8:40666/149‑1 (MQ=255)
                                                                                                              gggTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCAcc                                       >  8:245383/1‑149 (MQ=255)
                                                                                                              gggTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCAcc                                       <  2:40848/149‑1 (MQ=255)
                                                                                                                gTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCTTCCACCga                                     <  8:10426/149‑1 (MQ=255)
                                                                                                                      aGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCAc                                        <  2:70811/140‑1 (MQ=255)
                                                                                                                      aGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCAc                                        >  1:70811/1‑140 (MQ=255)
                                                                                                                       gggCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGcc                              >  4:34472/1‑149 (MQ=255)
                                                                                                                         gCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTg                            >  2:193135/1‑149 (MQ=255)
                                                                                                                            aTAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCg                          >  4:312381/1‑149 (MQ=255)
                                                                                                                             tAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCgg                         <  5:34816/149‑1 (MQ=255)
                                                                                                                                cTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGCCCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCa                      >  7:276569/1‑149 (MQ=255)
                                                                                                                                     gccgacGGCCTGGCCGGCGACCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGcc                 >  3:316862/6‑149 (MQ=255)
                                                                                                                                     gTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATcccc                                                   <  7:74868/115‑1 (MQ=255)
                                                                                                                                     gTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATcccc                                                   >  8:74868/1‑115 (MQ=255)
                                                                                                                                        gCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGtt                                                                   <  7:122259/96‑1 (MQ=255)
                                                                                                                                        gCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGtt                                                                   >  8:122259/1‑96 (MQ=255)
                                                                                                                                        gCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGc              >  8:73659/1‑149 (MQ=255)
                                                                                                                                         ccGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCg                                <  1:192024/130‑1 (MQ=255)
                                                                                                                                         ccGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCg                                >  2:192024/1‑130 (MQ=255)
                                                                                                                                         ccGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGCt             >  4:215660/1‑149 (MQ=255)
                                                                                                                                            gCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGCTTgg          <  3:34472/149‑1 (MQ=255)
                                                                                                                                               aGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGCTTggg         >  7:69734/1‑147 (MQ=255)
                                                                                                                                                    gTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGCTTGGGCAGCCCg  <  3:319430/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCCAGCCGTGTTGATATGCGGCGGGAACAGCGTAGTGCCGCCGGTAGGGCCGACTGCCGAGGTGATCGCCGATGGCGACAGGTTCATCGCCGGGATGGCGCGGTCGGTACGGGTTTGCAGGGCATAGCTGAGGTCGCCGGCCAGGCCGTCGGCCTGCATGCCACGCTTGCCGTCCGGGCCGACCTTGTTCGGGTCATACACGCGCAGGGCCAGCGGCTCGTTGCGGTAGTTGACGACGAACATCCCCGGGTCATCCACCGAGATCGCCTGCGGGCATGGCCGGCTTGGGCAGCCCG  >  NC_002947/3956539‑3956834

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