Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,007,107 T→C 7.4% I14T (ATC→ACC)  PP_5674 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,007,1070TC 7.4% 150.9 / 2.5 54I14T (ATC→ACC) PP_5674hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (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 = 6.53e-01

ACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGAACAACATGCAGAGCT  >  NC_002947/5006967‑5007254
                                                                                                                                            |                                                                                                                                                   
aCCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCg                                                                                                                                            >  2:570076/1‑150 (MQ=255)
      tCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAg                                                                                                                                                 <  2:145407/139‑1 (MQ=255)
      tCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAg                                                                                                                                                 >  1:145407/1‑139 (MQ=255)
        gatgTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGCAGCCGCCTGATGGCACAAAACTAAAAATCCTTGTGGCAACATAGCAAGATCCAGAAGCGTTTAGAcc                                                                                                                                    >  2:1477043/4‑149 (MQ=255)
        gACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACa                                                                                                                                    >  2:612947/1‑150 (MQ=255)
             cTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCTTTTCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATc                                                                                                                               <  1:1477043/150‑1 (MQ=255)
                aTGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGaa                                                                                                                            >  2:1713366/1‑150 (MQ=255)
                     tCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAAc                                                                                                                       >  1:1891569/1‑150 (MQ=255)
                         ccTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTa                                                                                                                   <  2:123626/150‑1 (MQ=255)
                            aaaTTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATACTCGAAAAAACGCTAgcg                                                                                                                <  2:995411/150‑1 (MQ=255)
                               ttCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCAt                                                                                                             >  1:1147821/1‑150 (MQ=255)
                                    aaTCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGa                                                                                                        <  2:1891569/150‑1 (MQ=255)
                                     aTCAAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAACGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGAc                                                                                                       <  2:597963/150‑1 (MQ=255)
                                         aGCATATGTGTCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATACCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCaa                                                                                                   >  2:1449815/1‑150 (MQ=255)
                                                ggggCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGt                                                                                            >  2:450632/1‑150 (MQ=255)
                                                 gggCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGtt                                                                                           <  1:1713366/150‑1 (MQ=255)
                                                   gCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTaa                                                                                         <  2:1147821/150‑1 (MQ=255)
                                                    cATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAAc                                                                                        <  2:1039638/150‑1 (MQ=255)
                                                     aTTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGAAAAAATAAAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACt                                                                                       <  2:12303/150‑1 (MQ=255)
                                                          tGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACaaa                                                                                  >  1:1604574/1‑150 (MQ=255)
                                                             cAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATAGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACaaaaaa                                                                               <  2:406905/150‑1 (MQ=255)
                                                              aaTCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAg                                                                              >  2:749520/1‑150 (MQ=255)
                                                                tCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAgc                                                                                                                 <  2:747250/113‑1 (MQ=255)
                                                                tCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAgc                                                                                                                 >  1:747250/1‑113 (MQ=255)
                                                                 cTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATACCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTaa                                                                           <  1:1449815/150‑1 (MQ=255)
                                                                       ccGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACt                                                                     <  1:749520/150‑1 (MQ=255)
                                                                                 aaTCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTa                                                                       >  2:1073603/1‑138 (MQ=255)
                                                                                 aaTCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTa                                                                       <  1:1073603/138‑1 (MQ=255)
                                                                                       cGGAGGCCGCCTGATGAAAGTCTCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGc                                                     >  1:1232647/1‑150 (MQ=255)
                                                                                               gccTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACaaaaaa                                                                               <  2:170373/116‑1 (MQ=255)
                                                                                               gccTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACaaaaaa                                                                               >  1:170373/1‑116 (MQ=255)
                                                                                                    atgaAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGt                                                                                                                                           >  2:1137459/1‑51 (MQ=255)
                                                                                                    atgaAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGt                                                                                                                                           <  1:1137459/51‑1 (MQ=255)
                                                                                                      taaaGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACt                                      >  1:280297/2‑150 (MQ=255)
                                                                                                            aTCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCCACTGCAGCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGa                                <  1:612947/150‑1 (MQ=255)
                                                                                                            aTCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTCTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGa                                >  2:625203/1‑150 (MQ=255)
                                                                                                               ttCCGAACGTCAGTGGCAACATATCAATACCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGaaa                                             >  2:1601027/1‑134 (MQ=255)
                                                                                                               ttCCGAACGTCAGTGGCAACATATCAATACCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGaaa                                             <  1:1601027/134‑1 (MQ=255)
                                                                                                                  cGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATc                                                                                                                               >  2:1207654/1‑49 (MQ=255)
                                                                                                                  cGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATc                                                                                                                               <  1:1207654/49‑1 (MQ=255)
                                                                                                                         cAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATg                   <  1:450632/150‑1 (MQ=255)
                                                                                                                            tGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGaa                <  2:238969/150‑1 (MQ=255)
                                                                                                                            tGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGaa                <  2:1232647/150‑1 (MQ=255)
                                                                                                                              gCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGaaca              >  1:163105/1‑150 (MQ=255)
                                                                                                                               cAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGaacaa             <  2:1604574/150‑1 (MQ=255)
                                                                                                                               cAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGaacaa             >  1:1472406/1‑150 (MQ=255)
                                                                                                                               cAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGaacaa             >  2:1039940/1‑150 (MQ=255)
                                                                                                                                aaCATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTaa                                                                                         <  1:370042/73‑1 (MQ=255)
                                                                                                                                aaCATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTaa                                                                                         >  2:370042/1‑73 (MQ=255)
                                                                                                                                aaCATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTaa                                                                                         <  1:685932/73‑1 (MQ=255)
                                                                                                                                aaCATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTaa                                                                                         >  2:685932/1‑73 (MQ=255)
                                                                                                                                     atCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGAacaaca           >  2:1861992/1‑146 (MQ=255)
                                                                                                                                     atCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGAacaaca           <  1:1861992/146‑1 (MQ=255)
                                                                                                                                     atCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGATGTCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGAACAACATGCa       >  1:1089346/1‑150 (MQ=255)
                                                                                                                                          tatCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGAACAACATGCAGAGCt  <  2:280297/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
ACCAAGTCGACGTCTGATGGCTCGCCCTAAATTCAAAATCCAGCATATGGGGCATTGCTGACAATCTTGCCCCGCCAACTAAATCTCCGGAGGCCGCCTGATGAAAGTATCTTCCGAACGTCAGTGGCAACATATCAATATCCAGAAGCGTTTAGACATAATCGAAAAAACGCTAGCGCATCTCGACTCAAAAAATGTTAACTACAAAAAAGTAACCTACTTAGCTAGACACCTAGCGAATGAAATCACACTCGCGGAGCGCCGGGAGATGGAACAACATGCAGAGCT  >  NC_002947/5006967‑5007254

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