Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 226,661 A→T 8.0% intergenic (+68/‑65) PP_0174 → / → PP_0175 hypothetical protein/transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947226,6610AT 8.0% 142.6 / 3.4 50intergenic (+68/‑65)PP_0174/PP_0175hypothetical protein/transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base T (2/2);  total (24/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 = 6.56e-01

AGGAGTGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACC  >  NC_002947/226514‑226806
                                                                                                                                                   |                                                                                                                                                 
aGGAGTGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTgac                                                                                                                                                  >  3:249687/1‑149 (MQ=255)
       gCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAggg                                                                                                                                           >  4:283988/1‑149 (MQ=255)
        cGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAgggg                                                                                                                                          >  8:187741/1‑149 (MQ=255)
          ggACCGGGCGGTGCAGCTGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTCGTCAACAgggggc                                                                                                                                        <  8:307142/149‑2 (MQ=255)
           gACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGAt                                                                                                                                       >  2:157280/1‑149 (MQ=255)
           gACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGAt                                                                                                                                       <  3:3608/149‑1 (MQ=255)
              cGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTc                                                                                                                                    <  8:245462/149‑1 (MQ=255)
                 gCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCggg                                                                                                                                 <  7:273004/149‑1 (MQ=255)
                 gCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCggg                                                                                                                                 <  4:94769/149‑1 (MQ=255)
                     tGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGCCg                                                                                                                              >  3:56329/1‑148 (MQ=255)
                         gATGTGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTatat                                                                                                                         <  2:242077/149‑1 (MQ=255)
                         gATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGGGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTatat                                                                                                                         >  2:308652/1‑149 (MQ=255)
                          aTGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATAtt                                                                                                                        >  3:13287/1‑149 (MQ=255)
                              gcagcaCATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATa                                                                                                                    <  3:283988/149‑1 (MQ=255)
                                agcaCATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGt                                                                                                                             <  2:52466/138‑1 (MQ=255)
                                agcaCATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGt                                                                                                                             >  1:52466/1‑138 (MQ=255)
                                  cacaTTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTa                                                                                                                <  2:202548/149‑1 (MQ=255)
                                                    gAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACt                                                                                              >  4:160870/1‑149 (MQ=255)
                                                    gAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTCGTCAACAggggg                                                                                                                                         <  8:144850/106‑1 (MQ=255)
                                                    gAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTCGTCAACAggggg                                                                                                                                         >  7:144850/1‑106 (MQ=255)
                                                       gCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTgg                                                                                           <  1:157280/149‑1 (MQ=255)
                                                              gACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCtt                                                                                    <  8:104919/149‑1 (MQ=255)
                                                                 cAGATCCATCCCTGGGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCAt                                                                                 <  1:308652/149‑1 (MQ=255)
                                                                   gATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGt                                                                               <  7:187741/149‑1 (MQ=255)
                                                                   gATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGt                                                                               >  5:52602/1‑149 (MQ=255)
                                                                     tccatccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTcc                                                                             <  3:160870/149‑1 (MQ=255)
                                                                       catccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCt                                                                           >  5:241927/1‑149 (MQ=255)
                                                                        atccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCtt                                                                          <  4:249687/149‑1 (MQ=255)
                                                                         tccCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGTATATACTTGGTTCTCTTCATGTCCCtt                                                                          >  7:168964/1‑148 (MQ=255)
                                                                           ccTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGa                                                                        <  6:241927/148‑1 (MQ=255)
                                                                              gAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGAttt                                                                     >  1:276585/1‑148 (MQ=255)
                                                                               aGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTc                                                                    <  5:104998/148‑1 (MQ=255)
                                                                                     tctcTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGa                                                                                                     >  8:302959/1‑109 (MQ=255)
                                                                                     tctcTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGa                                                                                                     <  7:302959/109‑1 (MQ=255)
                                                                                                 tgactgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTcag                                                   >  6:325289/1‑147 (MQ=255)
                                                                                                     tgaAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTgg                                             >  8:137039/1‑149 (MQ=255)
                                                                                                        aaaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGcc                                           <  6:52602/148‑1 (MQ=255)
                                                                                                         aaaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGATTATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTg                                         <  4:56329/149‑1 (MQ=255)
                                                                                                          aaTGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTgg                                        <  7:137039/149‑1 (MQ=255)
                                                                                                            tGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTa                                      <  4:13287/149‑1 (MQ=255)
                                                                                                               acaacaTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTc                                   >  1:137437/1‑149 (MQ=255)
                                                                                                                caacaTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCg                                  >  8:89979/1‑149 (MQ=255)
                                                                                                                  acaTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCt                                <  3:96194/149‑1 (MQ=255)
                                                                                                                    aTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTgg                                        >  8:111437/1‑139 (MQ=255)
                                                                                                                    aTTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTgg                                        <  7:111437/139‑1 (MQ=255)
                                                                                                                                 aTCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTcagca                                                 >  4:98498/1‑117 (MQ=255)
                                                                                                                                 aTCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTcagca                                                 <  3:98498/117‑1 (MQ=255)
                                                                                                                                    cGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGgcggcg               >  6:14968/1‑148 (MQ=255)
                                                                                                                                        tcccccTGTCGTCAACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCt          >  5:7682/2‑149 (MQ=255)
                                                                                                                                            ccTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCa      >  5:347137/1‑149 (MQ=255)
                                                                                                                                                ttGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATAcc  >  7:121154/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGGAGTGGCGGGACCGGGCGGTGCAGATGGGCAGCACATTGGCCTTGACCATGAAGCTGGAAGACCAGATCCATCCCTGAGTATTTCTCTTCAGTGCTGACTGAAAAATGAACAACATTTAACCAGGCTATCCGTCGCCCCCTGTTGACGACAGGGGGATGTCGGGGCGTATATTGATAGTTAGCAAACTATGAATATACTTGGTTCTCTTCATGTCCCTTGATTTCCTGCGCCTGCAAGTCAGCAGTGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACC  >  NC_002947/226514‑226806

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: