Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,351,944 C→A 10.5% L87L (CTG→CTT)  PP_3824 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,351,9440CA10.5% 152.2 / 5.3 52L87L (CTG→CTT) PP_3824hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base A (4/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAG  >  NC_002947/4351805‑4352091
                                                                                                                                           |                                                                                                                                                   
gaggaCCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCAGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCGTGTCGATATCCAACAGCGCCGCt                                                                                                                                            >  4:8708/1‑149 (MQ=255)
   gaaccGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGAAGCTGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACAAGCGCCGCTgtg                                                                                                                                         >  3:42690/4‑147 (MQ=255)
          tGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCAc                                                                                                                                  >  6:118791/1‑149 (MQ=255)
                  aTCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGt                                                                                                                           >  3:301398/1‑148 (MQ=255)
                  aTCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGACAAGGCGCAGCGGGAAGTTATTCAGCCGCTCTTGAAGATTGATGCTGCTGACGCTGCCGAGCTTGTCGATATCCACAAGCGCCTCTGGGCCGGCACCTGAGGt                                                                                                                           >  1:208668/1‑148 (MQ=255)
                    cAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGAt                                                                                                                        >  4:18638/1‑149 (MQ=255)
                    cAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGAt                                                                                                                        >  4:22866/1‑149 (MQ=255)
                     aaCTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCt                                                                                                                                            >  1:184963/1‑128 (MQ=255)
                     aaCTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCt                                                                                                                                            <  2:184963/128‑1 (MQ=255)
                     aaCTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGAt                                                                                                                        <  2:151686/148‑1 (MQ=255)
                          ccGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCt                                                                                                                   <  4:245378/148‑1 (MQ=255)
                           cGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTg                                                                                                                  <  4:301398/148‑1 (MQ=255)
                                aTGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTAc                                                                                                            >  5:90211/1‑149 (MQ=255)
                                  gCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGa                                                                                                          >  5:89462/1‑149 (MQ=255)
                                             ccGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAgg                                                                                               <  3:18638/149‑1 (MQ=255)
                                             ccGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAgg                                                                                               <  3:22866/149‑1 (MQ=255)
                                                      gcACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAg                                                                                                      >  5:15131/1‑133 (MQ=255)
                                                      gcACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAg                                                                                                      <  6:15131/133‑1 (MQ=255)
                                                        acaTCGCCAAGGCGCAGCGGGAAGTTATACAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACg                                                                                    >  1:179059/1‑149 (MQ=255)
                                                          aTCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCg                                                                                                                <  1:188627/119‑1 (MQ=255)
                                                          aTCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCg                                                                                                                >  2:188627/1‑119 (MQ=255)
                                                           tCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACgctg                                                                                 <  6:133843/149‑1 (MQ=255)
                                                                        gaggtaaGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCGGCTGACGCGCTTGTCGATATCCTCCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGATGCTCTTTGGTCAGc                                                                    <  2:208668/144‑1 (MQ=255)
                                                                                  aTTCCGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCAcgc                                                           <  8:58486/148‑1 (MQ=255)
                                                                                      aGCCGCTCGTGTAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTc                                                                                          >  3:228778/1‑113 (MQ=255)
                                                                                      aGCCGCTCGTGTAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTc                                                                                          <  4:228778/113‑1 (MQ=255)
                                                                                      aGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCTGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATc                                                      <  1:60009/149‑1 (MQ=255)
                                                                                       gCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCg                                                     >  7:50153/1‑149 (MQ=255)
                                                                                       gCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCg                                                     <  6:90211/149‑1 (MQ=255)
                                                                                              gTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCt                                              >  7:182098/1‑149 (MQ=255)
                                                                                              gTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCt                                              >  1:47166/1‑149 (MQ=255)
                                                                                                 cAGTTTGAAGCTGCTGCCGCTGCCGCGCGTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCtt                                             <  3:180433/147‑1 (MQ=255)
                                                                                                  aGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTtctc                                          >  5:326747/1‑149 (MQ=255)
                                                                                                   gTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTcc                                         >  2:118525/1‑149 (MQ=255)
                                                                                                       gATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCttt                                      <  5:118791/148‑1 (MQ=255)
                                                                                                        aTGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGctc                                                                              <  2:222254/107‑1 (MQ=255)
                                                                                                        aTGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGctc                                                                              >  1:222254/1‑107 (MQ=255)
                                                                                                        aTGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTAc                                    <  1:118525/149‑1 (MQ=255)
                                                                                                        aTGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTAc                                    >  5:156613/1‑149 (MQ=255)
                                                                                                        aTGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACACTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTAc                                    <  6:89462/149‑1 (MQ=255)
                                                                                                                    ctgccgcGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATa                        <  6:156613/149‑1 (MQ=255)
                                                                                                                     tgccgcGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAg                       >  6:152536/1‑149 (MQ=255)
                                                                                                                     agcggcgcTGGTGGATATCGACAAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCg                             <  1:11491/139‑1 (MQ=255)
                                                                                                                     agcggcgcTGGTGGATATCGACAAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCg                             >  2:6398/5‑143 (MQ=255)
                                                                                                                     agcggcgcTGGTGGATATCGACAAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCg                             <  1:6398/139‑1 (MQ=255)
                                                                                                                     agcggcgcTGGTGGATATCGACAAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCg                             >  2:11491/5‑143 (MQ=255)
                                                                                                                               gTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGc             >  5:246173/1‑149 (MQ=255)
                                                                                                                               gTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGc             >  2:201663/1‑149 (MQ=255)
                                                                                                                                  gATATCCACCACCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGctct          >  8:11680/1‑149 (MQ=255)
                                                                                                                                   atatCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCtt         <  1:162328/149‑1 (MQ=255)
                                                                                                                                   atatCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCtt         >  2:160835/1‑149 (MQ=255)
                                                                                                                                   atatCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCtt         <  2:179059/149‑1 (MQ=255)
                                                                                                                                       ccaccaGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGGCGGCATAGGCGCGTGGGCTCTTGCCa     <  8:50153/149‑1 (MQ=255)
                                                                                                                                          ccaGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCagag  <  5:152536/149‑1 (MQ=255)
                                                                                                                                           caGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCa     >  1:153314/1‑145 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GAGGACCCGTTGCTTCAGATCAACTTCCGGAGATGCATTGCGTTGCCGCAGTGCGCACATCGCCAAGGCGCAGCGGGAAGTTATTCAGCCGCTCGTGCAGTTTGATGCTGCTGCCGCTGCCGCGCTTGTCGATATCCACCAGCGCCGCTGGGCCGGCACCTGAGGTGATTTTCTGCGGTACGATCAGCCGCAGGCCGTCGCCACGCTGCTCTTCGGTCAGCGGGTCACGCGTATCGGCAAGCTTCTCCTTTACGCAGTCGGCATAGGCGCGTGGGCTCTTGCCAGAG  >  NC_002947/4351805‑4352091

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