Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,769 A→C 7.2% K71Q (AAG→CAG)  PP_4562 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,7690AC 7.2% 153.7 / 4.5 56K71Q (AAG→CAG) PP_4562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base C (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

CAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAC  >  NC_002947/5182624‑5182915
                                                                                                                                                 |                                                                                                                                                  
cagcagCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGt                                                                                                                                                >  2:616007/1‑150 (MQ=255)
  gcagcaACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGAGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGtcc                                                                                                                                              >  2:325281/1‑148 (MQ=255)
     gcaACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCcca                                                                                                                                           <  2:257771/150‑1 (MQ=255)
      caACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTTTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAATGTGCtcac                                                                                                                                          <  1:708004/150‑1 (MQ=255)
            ccTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAAttt                                                                                                                                    >  1:806937/1‑150 (MQ=255)
               gTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCt                                                                                                                                 >  1:390858/1‑150 (MQ=255)
               gTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCt                                                                                                                                 >  1:465452/1‑150 (MQ=255)
                          tGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCAcc                                                                                                                      <  1:281325/150‑1 (MQ=255)
                           ggACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCt                                                                                                                     >  1:1146809/1‑150 (MQ=255)
                               aacaaACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAg                                                                                                                 >  1:894292/1‑150 (MQ=255)
                               aacaaACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTACTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAg                                                                                                                 >  1:605492/1‑150 (MQ=255)
                                 caaACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCa                                                                                                               >  2:1102772/1‑150 (MQ=255)
                                  aaaCGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTTCAGGTGctgca                                                                                                                                          <  2:1121607/122‑5 (MQ=255)
                                  aaaCGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTTCAGGTGctgca                                                                                                                                          >  1:1121607/1‑118 (MQ=255)
                                  aaaCGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCaa                                                                                                              >  1:382597/1‑150 (MQ=255)
                                  aaaCGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCaa                                                                                                              >  1:33478/1‑150 (MQ=255)
                                   aaCGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCaaa                                                                                                             <  1:635966/150‑1 (MQ=255)
                                    aCGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAAt                                                                                                            >  1:934505/1‑150 (MQ=255)
                                        aaCCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGaa                                                                                                        <  2:147031/150‑1 (MQ=255)
                                        aaCCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGaa                                                                                                        <  2:33478/150‑1 (MQ=255)
                                          ccTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAAGTTCCTCAAGGCTCACCTTGAGCAAATGGAAGt                                                                                                      <  2:92953/150‑1 (MQ=255)
                                                gtggtTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCaa                                                                                                >  2:907496/1‑150 (MQ=255)
                                                 tggtTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAg                                                                                               >  1:523357/1‑150 (MQ=255)
                                                    ttCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCt                                                                                            >  1:268292/1‑150 (MQ=255)
                                                     tCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTg                                                                                           <  1:585940/150‑1 (MQ=255)
                                                                 tGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCt                                                                               <  2:268292/150‑1 (MQ=255)
                                                                     aCTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTccac                                                                           >  2:142209/1‑150 (MQ=255)
                                                                       tgctgGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTgc                                                                                                    <  1:1205955/123‑1 (MQ=255)
                                                                       tgctgGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTgc                                                                                                    >  2:1205955/1‑123 (MQ=255)
                                                                           gggCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCGACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAg                                                                     >  1:1221628/1‑150 (MQ=255)
                                                                                  ggCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGcc                                                              <  2:947299/150‑1 (MQ=255)
                                                                                          cGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCGACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAgcgc                                                      <  2:1221628/150‑1 (MQ=255)
                                                                                          cGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAgcgc                                                      <  2:894292/150‑1 (MQ=255)
                                                                                          cGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAgcgc                                                      >  2:1161683/1‑150 (MQ=255)
                                                                                           gACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTc                                                                                                                                >  2:1418709/1‑75 (MQ=255)
                                                                                           gACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTc                                                                                                                                <  1:1418709/75‑1 (MQ=255)
                                                                                           gACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGcc                                                     <  2:523357/150‑1 (MQ=255)
                                                                                                      ccTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTg                                                                 <  2:1000879/127‑1 (MQ=255)
                                                                                                      ccTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTg                                                                 >  1:1000879/1‑127 (MQ=255)
                                                                                                           tgctgcACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGaa                                                                                                        <  1:545467/83‑1 (MQ=255)
                                                                                                           tgctgcACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGaa                                                                                                        >  2:545467/1‑83 (MQ=255)
                                                                                                              tgcACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGcc                                                     <  2:624530/131‑1 (MQ=255)
                                                                                                              tgcACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGcc                                                     >  1:624530/1‑131 (MQ=255)
                                                                                                              tgcACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCtgctgc                                  <  2:382597/150‑1 (MQ=255)
                                                                                                                     gCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCGTCACAAGCCCCGATCGACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGacac                           >  1:526394/1‑150 (MQ=255)
                                                                                                                     gCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGacac                           <  1:1400361/150‑1 (MQ=255)
                                                                                                                             gCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                   <  2:605492/150‑1 (MQ=255)
                                                                                                                             gCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCgg                   <  1:1102772/150‑1 (MQ=255)
                                                                                                                                 tGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTgg               <  1:1161683/150‑1 (MQ=255)
                                                                                                                                  gAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGCCACCCGGACGGGTgtc              <  2:1364582/150‑3 (MQ=255)
                                                                                                                                        gggcaCCTTCAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGcc                                                              <  2:987637/94‑1 (MQ=255)
                                                                                                                                        gggcaCCTTCAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGcc                                                              >  1:987637/3‑96 (MQ=255)
                                                                                                                                        cagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagcag        <  1:907496/150‑1 (MQ=255)
                                                                                                                                         agcaCCTGAAGGTGCCCACGAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagcagc       <  2:934505/150‑1 (MQ=255)
                                                                                                                                          gcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCt      <  2:465452/150‑1 (MQ=255)
                                                                                                                                              cTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAc  >  1:886256/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAC  >  NC_002947/5182624‑5182915

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