Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 333,174 C→T 7.8% H105Y (CAC→TAC)  PP_5435 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947333,1740CT 7.8% 145.0 / 2.3 51H105Y (CAC→TAC) PP_5435hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.39e-01

GATCGCATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCT  >  NC_002947/333027‑333321
                                                                                                                                                   |                                                                                                                                                   
gATCGCATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGc                                                                                                                                                     >  8:311488/1‑148 (MQ=255)
     cATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGcc                                                                                                                                                >  6:190082/1‑148 (MQ=255)
          cccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                          >  6:260919/1‑149 (MQ=255)
           ccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                          >  3:40668/1‑148 (MQ=255)
           ccACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTAc                                                                                                                                          >  7:18560/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                                                                                    >  3:71211/1‑136 (MQ=255)
                                                                             aaaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCa                                                                                    <  4:71211/136‑1 (MQ=255)
                                                                              aaCGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGa                                                                       >  8:180448/1‑148 (MQ=255)
                                                                                cGAATGCAGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGacaa                                                                    <  4:40668/149‑1 (MQ=255)
                                                                                           tGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                              >  1:204160/1‑128 (MQ=255)
                                                                                           tGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGa                                                                              <  2:204160/128‑1 (MQ=255)
                                                                                                          aTGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                            >  1:264261/1‑147 (MQ=255)
                                                                                                          aTGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                            <  2:264261/147‑1 (MQ=255)
                                                                                                                    gggTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGa                                <  5:260919/149‑1 (MQ=255)
                                                                                                                      gTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCGACTGGACTAGATg                              >  8:105210/1‑149 (MQ=255)
                                                                                                                        aCTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                            >  2:175019/1‑133 (MQ=255)
                                                                                                                        aCTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGc                                            <  1:175019/133‑1 (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                              >  3:261201/5‑131 (MQ=255)
                                                                                                                                        accatttGTCGTACTCGGCGTGCGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATg                              <  4:261201/127‑1 (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)
                                                                                                                                                   |                                                                                                                                                   
GATCGCATCCCCCACTCAAAACTCAGGCCCGCCAGCCGCAAAGGGCCGCTTGATATCTACGTGTTCGATATCAAGAAAAACGAATGCCGCCTGATCGCGTGGCTCAATGCGCGAACGGGTACTCTCTACATCAAGGCCATTTTGTCGCACGCCGAGTACGACAAATGGTGCAGAAGCGATATCAGGTAGCAACCTGGGCCGCGCGATATTTCAATTCGACAGGTGACAACATGAAGCGGGATAAAATCGAAGCCAATGGACTAGATGCCTTCATCGCTAATATTTCACCCATGCT  >  NC_002947/333027‑333321

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