Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,486,260 A→C 7.5% S30R (AGC→CGC)  PP_3978 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,486,2600AC 7.5% 144.2 / 2.8 53S30R (AGC→CGC) PP_3978hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base C (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

ACGGCAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAT  >  NC_002947/4486111‑4486400
                                                                                                                                                     |                                                                                                                                            
aCGGCAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGa                                                                                                                                              <  1:1066669/150‑1 (MQ=255)
           aaTGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGaa                                                                                                                                   >  1:781523/1‑150 (MQ=255)
           aaTGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGaa                                                                                                                                   <  2:781523/150‑1 (MQ=255)
                   gcCACTACAATTACATAGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCt                                                                                                                           >  1:253812/1‑150 (MQ=255)
                    ccACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTg                                                                                                                          >  2:468043/1‑150 (MQ=255)
                       cTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACGGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATg                                                                                                                       >  2:1327812/1‑150 (MQ=255)
                        tACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATgg                                                                                                                      <  1:163285/150‑1 (MQ=255)
                         aCAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTAGCTGCCATGAAGCCGCCCTGATGGa                                                                                                                     <  1:398293/150‑1 (MQ=255)
                               aCATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCCGCTATGAGCTCCACAGGCACACAAAGCATGACCCTGCAGCGAGCTGCCATGAAGCCGGCCTGATGGAGCAttt                                                                                                               <  2:651044/150‑1 (MQ=255)
                                cATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCCCAAAAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTg                                                                                                              >  2:1112555/1‑150 (MQ=255)
                                 aTCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCCCCCTGATGGAGCATTTGc                                                                                                             >  2:650436/1‑150 (MQ=255)
                                          gacCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                      >  1:129700/1‑114 (MQ=255)
                                          gacCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                      <  2:129700/116‑3 (MQ=255)
                                                  cATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACg                                                                                            <  1:1102430/150‑1 (MQ=255)
                                                  cATTCCATAGGTATGGGCTTCATGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACg                                                                                            <  2:253812/150‑1 (MQ=255)
                                                       cATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGCAGGAGTTGCACCGTGGCTTTAAGACGGGATc                                                                                       <  1:1327812/150‑1 (MQ=255)
                                                        aTAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGt                                                                                                        <  2:849713/132‑1 (MQ=255)
                                                        aTAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGt                                                                                                        >  1:849713/1‑132 (MQ=255)
                                                        aTAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGATATGAGCTCCACAGGCACACAAAGCATGACCCGACAGCGAGCTGCCATGAAGCCGCCCTGATGGGGCAATTGCACCGTGGCTTTAAGTCGGGATCg                                                                                      >  1:5286/1‑150 (MQ=255)
                                                            gtctgtgCCACCCCTGGGGAAGCCAGCCACGCACACTCGATCGCAGCAGATCAGAACTCCACAAGCACACATAGCATGACCCGGCAGCGAGCTGCCAGGAAGCCGCCCCGATGGAGCATTTCCACCGTAGGTTTACCACGGCATCGAtcc                                                                                  >  2:22183/4‑148 (MQ=255)
                                                            gTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                      <  1:740316/98‑3 (MQ=255)
                                                            gTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCTCGCTGccgg                                                                                                                                      >  2:740316/1‑96 (MQ=255)
                                                            gTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCt                                                                                                                                           <  1:1310715/93‑1 (MQ=255)
                                                            gTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCt                                                                                                                                           >  2:1310715/1‑93 (MQ=255)
                                                                  gCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAg                                                                                                                                  >  1:202064/1‑96 (MQ=255)
                                                                  gCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAg                                                                                                                                  <  2:202064/96‑1 (MQ=255)
                                                                  gCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTAAGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAg                                                                            <  1:320370/150‑1 (MQ=255)
                                                                      cccGTGGTGTAGCCATAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGc                                                                                                                   >  2:60020/1‑107 (MQ=255)
                                                                      cccGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGc                                                                                                                   <  1:60020/107‑1 (MQ=255)
                                                                      cccGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACACGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATATGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTggg                                                                        <  1:1362528/150‑1 (MQ=255)
                                                                       ccGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGc                                                                       >  1:1225008/1‑150 (MQ=255)
                                                                       ccGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGCATCGGTGCCGATATTGGGc                                                                       >  2:36242/1‑150 (MQ=255)
                                                                       ccGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAACTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACTGACTGTgggg                                                                       >  1:283059/1‑149 (MQ=255)
                                                                              gtAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATg                                                                 <  2:1138261/149‑1 (MQ=255)
                                                                              gtAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATg                                                                 >  1:1138261/1‑149 (MQ=255)
                                                                                   cacaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTaa                                                           >  1:1282083/1‑150 (MQ=255)
                                                                                   cacaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTaa                                                           <  2:251055/150‑1 (MQ=255)
                                                                                   cacaAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCAATGAAGCCGCCCTGATGGAGCATTTGCATCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTaa                                                           <  1:22183/150‑1 (MQ=255)
                                                                                     caAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGAACCATGACTTTAAGACGGGATCGATACCTACAGTGGGCTACATGGCGTAAGc                                                         <  2:557562/150‑1 (MQ=255)
                                                                                        aCGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTc                                                      >  2:392012/1‑150 (MQ=255)
                                                                                               aTCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                               <  2:952984/150‑1 (MQ=255)
                                                                                               aTCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                               <  2:1034821/150‑1 (MQ=255)
                                                                                               aTCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGaa                                               >  1:316631/1‑150 (MQ=255)
                                                                                                   aTGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCa                                           <  2:316631/150‑1 (MQ=255)
                                                                                                          cagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTggg                                                                        >  2:365842/1‑114 (MQ=255)
                                                                                                          cagcTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTggg                                                                        <  1:365842/114‑1 (MQ=255)
                                                                                                                   gCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCAc                                                                                                           <  1:714130/70‑1 (MQ=255)
                                                                                                                   gCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCAc                                                                                                           >  2:714130/1‑70 (MQ=255)
                                                                                                                         caGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAg                     <  2:1282083/150‑1 (MQ=255)
                                                                                                                          aGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAgg                    >  2:1067964/1‑150 (MQ=255)
                                                                                                                            gCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGTGTGATGCGAACTCACCTTCCGAACTCCAGCATTAAGGGc                  >  2:473570/1‑150 (MQ=255)
                                                                                                                            gCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGc                  >  1:333458/1‑150 (MQ=255)
                                                                                                                                        tGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCcgc      <  1:1067964/150‑1 (MQ=255)
                                                                                                                                            ccGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAt  >  1:1035889/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
ACGGCAAGATCAATGCTGCGCCACTACAATTACATCGTTGACGACCCGATCATTCCATAGGTATGTGCCACCCGTGGTGTAGCCACAAACGCAACATCGATGGCAGCAGCTATGAGCTCCACAGGCACACAAAGCATGACCCGGCAGCGAGCTGCCATGAAGCCGCCCTGATGGAGCATTTGCACCGTGGCTTTAAGACGGGATCGATACCGACAGTGGGCTACATGGCGTAAGCGTCCGGCGAACTCACCTTCCGAACTCCAGCATTAAGGGCGATGGTGTCCGCGTAT  >  NC_002947/4486111‑4486400

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