Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,959,653 G→A 8.0% L244L (CTC→CTT)  mnxG ← manganese‑oxidizing multicopper oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,959,6530GA 8.0% 144.0 / 2.6 50L244L (CTC→CTT) mnxGmanganese‑oxidizing multicopper oxidase
Reads supporting (aligned to +/- strand):  ref base G (26/20);  new base A (2/2);  total (28/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 = 9.46e-01

TTGCCGTTGACCATGGCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGCCCCTTCGCGCTCG  >  NC_002947/3959507‑3959791
                                                                                                                                                  |                                                                                                                                           
ttGCCGTTGACCATGGCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAg                                                                                                                                           >  6:167764/1‑149 (MQ=255)
         aCCATGGCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCg                                                                                                                                  >  7:227213/1‑149 (MQ=255)
           catcccAATCACCGGATCGCCGGCCTCGCCCTCGCCGACCCACGCTTCGACCTGCCGATAACGGGCTTCCTACTCCATCACGGGCTGGCCGGTGGCGGATCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCACTGAGGGTGTTGCggt                                                                                                                                <  4:322510/144‑1 (MQ=255)
           cATGGCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCgg                                                                                                                                 >  3:322510/1‑148 (MQ=255)
           cATGGCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCgg                                                                                                                                 >  1:336872/1‑148 (MQ=255)
              ggCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCggtagg                                                                                                                             >  5:139796/1‑149 (MQ=255)
              ggCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCggtagg                                                                                                                             >  2:311391/1‑149 (MQ=255)
                 aaTGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTAc                                                                                                                          >  2:75884/1‑149 (MQ=255)
                         gTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGt                                                                                                                  <  1:98095/149‑1 (MQ=255)
                         gTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGt                                                                                                                  >  2:98095/1‑149 (MQ=255)
                              ccGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTggg                                                                                                             >  2:128459/1‑149 (MQ=255)
                                   cttgccttCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑Ga                                                                                                       <  8:227213/149‑1 (MQ=255)
                                         ttCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑Gaccacca                                                                                                 >  4:255630/1‑149 (MQ=255)
                                             gCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTg                                                                                                               >  5:291631/1‑132 (MQ=255)
                                             gCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTg                                                                                                               <  6:291631/132‑1 (MQ=255)
                                                 tCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGcc                                                                                          >  3:307662/1‑148 (MQ=255)
                                                 tCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCg                                                                                         <  3:255630/149‑1 (MQ=255)
                                                   cACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTg                                                                                             >  8:52852/1‑143 (MQ=255)
                                                   cACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTg                                                                                             <  7:52852/143‑1 (MQ=255)
                                                     cGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACa                                                                                     <  6:139796/149‑1 (MQ=255)
                                                      ggTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAg                                                                                    <  1:128459/149‑1 (MQ=255)
                                                           gACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑Gaccacc                                                                                                  <  8:197998/130‑1 (MQ=255)
                                                           gACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑Gaccacc                                                                                                  >  7:197998/1‑130 (MQ=255)
                                                             ccTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTg                                                                             <  5:167764/149‑1 (MQ=255)
                                                                 cGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGga                                                                           <  1:75884/147‑1 (MQ=255)
                                                                       aCGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTg                                                                    <  2:336872/148‑1 (MQ=255)
                                                                          ggCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCccg                                                                >  8:309447/1‑149 (MQ=255)
                                                                             ttCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCccgccg                                                             >  5:320333/1‑149 (MQ=255)
                                                                                      gATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGcccaaac                                                    >  5:32427/1‑145 (MQ=255)
                                                                                      gATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTcgc                                                    >  3:329222/1‑149 (MQ=255)
                                                                                                tGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACg                                           >  5:111856/1‑148 (MQ=255)
                                                                                                  gCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTc                                                      >  2:344763/1‑135 (MQ=255)
                                                                                                  gCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTc                                                      <  1:344763/135‑1 (MQ=255)
                                                                                                  gCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTa                                         <  6:111856/148‑1 (MQ=255)
                                                                                                      gTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGc                                    >  4:82814/1‑149 (MQ=255)
                                                                                                          cGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTg                                <  1:311391/149‑1 (MQ=255)
                                                                                                                gCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCtggtg                          <  4:307662/149‑1 (MQ=255)
                                                                                                                gCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCtggtg                          <  7:309447/149‑1 (MQ=255)
                                                                                                                 ccccGTGGTGGCCAGGCGCATTTCTTCTTCAGTAAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTcgcc                                                   >  5:160868/1‑123 (MQ=255)
                                                                                                                 ccccGTGGTGGCCAGGCGCATTTCTTCTTCAGTAAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTcgcc                                                   <  6:160868/123‑1 (MQ=255)
                                                                                                                 ccccGTGGTGGCCAGGCGCATTTCTTCTTCAGTAAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTcgcc                                                   >  7:172325/1‑123 (MQ=255)
                                                                                                                 ccccGTGGTGGCCAGGCGCATTTCTTCTTCAGTAAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTcgcc                                                   <  8:172325/123‑1 (MQ=255)
                                                                                                                  cccGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTggc                        >  8:313772/1‑149 (MQ=255)
                                                                                                                   ccGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGACGAACGAAGCGGCTTGCCTGGTggcg                       >  4:302541/1‑149 (MQ=255)
                                                                                                                                  gcATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGCCCCTTcg        >  5:128442/1‑149 (MQ=255)
                                                                                                                                    aTTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGCCCCTTcgcg      <  1:253364/149‑1 (MQ=255)
                                                                                                                                       tcttcttcAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGc               <  4:21637/137‑1 (MQ=255)
                                                                                                                                       tcttcttcAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGc               >  3:21637/1‑137 (MQ=255)
                                                                                                                                        cttcttcAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGCCCCTTCGCGCTc   <  3:82814/148‑1 (MQ=255)
                                                                                                                                         ttcttcAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCACGGACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGCCCCTTCGCGCTCg  >  1:76775/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
TTGCCGTTGACCATGGCAATGATCGGTTTGCCGGCCTTGCCTTCGGCGATCCACGGTTCGACCTGCGGATAACGGGCTTCGTAGTCGATCACGGGCTGGCCGGTGGCGGTTCGCCCCGTGGTGGCCAGGCGCATTTCTTCTTCAGTGAGGGTGTTGCGGTAGGTACGGCCCAGTTTGGGCAC‑GACCACCACCTGGCCGAACAGGCCGTTGGAGACGTTGCCCGCCGTGCCCTCGCCGCCGAACGTAGCGGCTTTGCTGGTGGCGGCGAAGGCCCCTTCGCGCTCG  >  NC_002947/3959507‑3959791

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: ‑