Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 333,184 A→G 7.4% Y108C (TAC→TGC)  PP_5435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947333,1840AG 7.4% 155.9 / 2.4 54Y108C (TAC→TGC) PP_5435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base G (2/2);  total (25/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 = 9.97e-01

CCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAT  >  NC_002947/333037‑333325
                                                                                                                                                   |                                                                                                                                             
cccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                              >  6:260919/1‑149 (MQ=255)
 ccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                              >  7:18560/1‑148 (MQ=255)
 ccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                              >  3:40668/1‑148 (MQ=255)
     tCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACa                                                                                                                                          >  6:21725/1‑148 (MQ=255)
       aaaaCTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAAt                                                                                                                                       >  8:4815/1‑149 (MQ=255)
        aaaCTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATg                                                                                                                                      <  1:4754/149‑1 (MQ=255)
         aaCTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATgg                                                                                                                                     <  5:154341/149‑1 (MQ=255)
            tCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTg                                                                                                                                   >  5:18938/1‑148 (MQ=255)
            tCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCCCCTGATCGCGTGGCTCAATGCACGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCCCGCCGAGTACGACAAATGGTGc                                                                                                                                  <  6:284674/149‑1 (MQ=255)
              aGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGc                                                                                                                                  <  7:138220/147‑1 (MQ=255)
              aGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGc                                                                                                                                  <  4:203498/147‑1 (MQ=255)
              aGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCa                                                                                                                                 <  8:18560/148‑1 (MQ=255)
                 cccGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCGGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGaa                                                                                                                              <  4:335585/148‑1 (MQ=255)
                   cGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCg                                                                                                                           <  6:18938/149‑1 (MQ=255)
                    gCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGa                                                                                                                          >  7:235435/1‑149 (MQ=255)
                             aaaGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTa                                                                                                                 <  7:301960/149‑1 (MQ=255)
                             aaaGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGACTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTa                                                                                                                 <  7:311488/149‑1 (MQ=255)
                                                 cGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGa                                                                                               <  8:235435/147‑1 (MQ=255)
                                                   tgtTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATAtt                                                                                           <  6:233458/149‑1 (MQ=255)
                                                    gtTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGCACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATAttt                                                                                          <  3:249301/149‑1 (MQ=255)
                                                         atatCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAAtt                                                                                     <  7:4815/149‑1 (MQ=255)
                                                          tatCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTc                                                                                    <  5:21725/149‑1 (MQ=255)
                                                              aaGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                                  >  1:380422/1‑147 (MQ=255)
                                                              aaGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACa                                                                                <  2:380422/149‑1 (MQ=255)
                                                               aGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCcgc                                                                                                   <  6:354516/129‑1 (MQ=255)
                                                               aGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAg                                                                               <  5:190082/149‑1 (MQ=255)
                                                                gAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAgg                                                                              >  7:233350/1‑149 (MQ=255)
                                                                   aaaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCa                                                                                        <  4:71211/136‑1 (MQ=255)
                                                                   aaaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCa                                                                                        >  3:71211/1‑136 (MQ=255)
                                                                    aaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGa                                                                           >  8:180448/1‑148 (MQ=255)
                                                                      cGAATGCAGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGacaa                                                                        <  4:40668/149‑1 (MQ=255)
                                                                                 tGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                                  <  2:204160/128‑1 (MQ=255)
                                                                                 tGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                                  >  1:204160/1‑128 (MQ=255)
                                                                                                aTGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                <  2:264261/147‑1 (MQ=255)
                                                                                                aTGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                >  1:264261/1‑147 (MQ=255)
                                                                                                          gggTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGa                                    <  5:260919/149‑1 (MQ=255)
                                                                                                            gTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCGACTGGACTAGATg                                  >  8:105210/1‑149 (MQ=255)
                                                                                                              aCTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                <  1:175019/133‑1 (MQ=255)
                                                                                                              aCTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                                >  2:175019/1‑133 (MQ=255)
                                                                                                                  tctACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAAAGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCg                                                               >  6:335947/1‑114 (MQ=255)
                                                                                                                       aTCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTg                                                                                                         >  8:107035/1‑67 (MQ=255)
                                                                                                                          aaGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTaa                     <  7:180448/148‑1 (MQ=255)
                                                                                                                              ccATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTaa                     >  2:124753/1‑144 (MQ=255)
                                                                                                                              ccATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTaa                     <  1:124753/144‑1 (MQ=255)
                                                                                                                              accatttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATg                                  <  4:261201/127‑1 (MQ=255)
                                                                                                                              accatttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATg                                  >  3:261201/5‑131 (MQ=255)
                                                                                                                                  tttGTCGTACTCGGCGTGCGACAAATGGTTCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAAt                                                                                      <  1:268393/75‑1 (MQ=25)
                                                                                                                                  tttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAAt                                                                                      >  2:268393/1‑75 (MQ=37)
                                                                                                                                       cGCACGCCGAGTACGACAAATTGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGc       >  6:138131/1‑149 (MQ=255)
                                                                                                                                        gCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATATCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCAATGCt      >  4:47143/1‑149 (MQ=255)
                                                                                                                                            gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTggg                                                                                                       >  4:273385/1‑48 (MQ=255)
                                                                                                                                            gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTggg                                                                                                       <  3:273385/48‑1 (MQ=255)
                                                                                                                                            gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAt  >  6:380132/1‑149 (MQ=255)
                                                                                                                                            gCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAt  >  8:53262/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCTCGAT  >  NC_002947/333037‑333325

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