Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,847,352 C→A 8.1% A686D (GCT→GAT)  PP_4261 → cation‑transporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,847,3520CA 8.1% 151.7 / 2.8 49A686D (GCT→GAT) PP_4261cation‑transporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base A (2/2);  total (24/25)
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.68e-01

TGCAAGGCGCGGGGCTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTGCTGTCCAACCGCCTGCAAGC  >  NC_002947/4847207‑4847499
                                                                                                                                                 |                                                                                                                                                   
tgcAAGGCGCGGGGCTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGAc                                                                                                                                                 <  2:99873/150‑1 (MQ=255)
 gcAAGGCGCGGGGCTGGCACACCATGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGATCGTCTGAAAGCCCTGCAGGCTGACg                                                                                                                                                >  2:555492/1‑150 (MQ=255)
        gcgGGGCTGGCACTCCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGATGGCTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCaa                                                                                                                                         <  2:367310/150‑1 (MQ=255)
           gggCTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGt                                                                                                                                      <  1:752207/150‑1 (MQ=255)
            ggCTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTg                                                                                                                                     >  2:1090402/1‑150 (MQ=255)
              cTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGATGACGGGCGCAAGGTGCt                                                                                                                                   <  1:825135/150‑1 (MQ=255)
              cTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGATGACGGGCGCAAGGTGCt                                                                                                                                   <  1:825134/150‑1 (MQ=255)
              cTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGATGACGGGCGCAAGGTGCt                                                                                                                                   >  2:825134/1‑150 (MQ=255)
              cTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGATGACGGGCGCAAGGTGCt                                                                                                                                   >  2:825135/1‑150 (MQ=255)
               tGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTg                                                                                                                                  <  2:903053/150‑1 (MQ=255)
                ggCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGa                                                                                                                                 >  2:410068/1‑150 (MQ=255)
                 gCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGAt                                                                                                                                >  1:1105591/1‑150 (MQ=255)
                     acCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCCGGCTGACGGGCGCAAGGTGCTGCTGCTc                                                                                                                            >  2:197010/1‑150 (MQ=255)
                       ccTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTAATGCTCgg                                                                                                                          <  2:645908/150‑1 (MQ=255)
                       ccTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGAACAATATACGGGACCGCCTGAAAGCCCTTCAGGCTGACGGGCGCAAGGTGCTGATGCTCgg                                                                                                                          >  1:182767/1‑150 (MQ=255)
                               tgTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACgggg                                                                                                                  >  1:819199/1‑150 (MQ=255)
                                       gACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGAt                                                                                                          <  2:420902/150‑1 (MQ=255)
                                           gCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTAc                                                                                                      <  2:265726/150‑1 (MQ=255)
                                                 cGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTAc                                                                                                <  1:1057959/150‑1 (MQ=255)
                                                                   ccGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCa                                                                              <  1:410068/150‑1 (MQ=255)
                                                                             gggTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCAt                                                                    <  1:323823/150‑1 (MQ=255)
                                                                               gTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACACACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGTCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTATCGGTACTGGCCGCCGCCGATATCAGCCTCGCCATgg                                                                  >  2:265930/1‑150 (MQ=255)
                                                                               gTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATgg                                                                  <  2:1041605/150‑1 (MQ=255)
                                                                                 ctgcaCCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGc                                                                >  2:212530/5‑150 (MQ=255)
                                                                                  tCGACCAGGCCATTGGCGGCCTTCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCa                                                               <  1:56575/150‑1 (MQ=255)
                                                                                  tCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCa                                                               <  2:756851/150‑1 (MQ=255)
                                                                                        aGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAgc                                                             >  2:747772/1‑146 (MQ=255)
                                                                                        aGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAgc                                                             <  1:747772/146‑1 (MQ=255)
                                                                                               tGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCt                                                  <  2:310727/150‑1 (MQ=255)
                                                                                                          gcCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCa                                       <  2:524668/150‑1 (MQ=255)
                                                                                                            ccGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGgcc                                                                                           <  1:313781/96‑1 (MQ=255)
                                                                                                            ccGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGgcc                                                                                           >  2:313781/1‑96 (MQ=255)
                                                                                                                  gACAAACTGGACCGCCTGAAATCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCCGGCCAAGACCAGTGCCggc                               >  2:369047/1‑148 (MQ=255)
                                                                                                                    cAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGc                             >  1:338295/1‑150 (MQ=255)
                                                                                                                     aaaCTGGACCGCCTGAAAGCCCTGCAGGCTTACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAgc                                                             <  2:563881/117‑1 (MQ=255)
                                                                                                                     aaaCTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAgc                                                             >  1:563881/1‑117 (MQ=255)
                                                                                                                        cTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTc                         <  2:249470/150‑1 (MQ=255)
                                                                                                                          ggACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGTGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCtt                       <  1:800430/150‑1 (MQ=255)
                                                                                                                           gACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTg                      >  2:839008/1‑150 (MQ=255)
                                                                                                                                ccTGAAAGCCCTGCCGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTAc                                                                                                >  2:887809/1‑71 (MQ=255)
                                                                                                                                ccTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTAc                                                                                                <  1:887809/71‑1 (MQ=255)
                                                                                                                                 cTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCAATCTCGCTggac                >  2:612728/1‑147 (MQ=255)
                                                                                                                                   gAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTGCTGTCCaa              <  2:736758/150‑1 (MQ=255)
                                                                                                                                      aGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTGCTGTCCAACCg           >  1:1033144/1‑150 (MQ=255)
                                                                                                                                          cTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTGCTGTCCAACCGCCTg       >  1:642663/1‑150 (MQ=255)
                                                                                                                                           tGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTGCTGTCCAACCGCCtac      >  1:796602/1‑148 (MQ=255)
                                                                                                                                               ggCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCGGTGGCGACGCCGTCTTGCTGTCCAACCGCCTGCAAgc  >  2:932307/1‑150 (MQ=255)
                                                                                                                                                gCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATc                                                                               >  1:321860/1‑72 (MQ=255)
                                                                                                                                                gCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATc                                                                               <  2:321860/72‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGCAAGGCGCGGGGCTGGCACACCCTGCTGCTGTCGGGCGACAGCTCGCCGATGGTTGCCGAAGTGGCCGCGCAGCTGGGTATCGACCAGGCCATTGGCGGCCTGCGCCCGGATGACAAACTGGACCGCCTGAAAGCCCTGCAGGCTGACGGGCGCAAGGTGCTGATGCTCGGCGACGGGGTCAACGATGTACCGGTACTGGCCGCCGCCGATATCAGCATCGCCATGGGCAGCGCCACCGACCTGGCCAAGACCAGTGCCGACGCCGTCTTGCTGTCCAACCGCCTGCAAGC  >  NC_002947/4847207‑4847499

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