Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 107,990 Δ1 bp 7.0% coding (451/642 nt) PP_5426 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947107,9900C. 7.0% 182.8 / 7.0 57coding (451/642 nt)PP_5426hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base . (2/2);  total (29/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 = 9.42e-01

CACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAG  >  NC_002947/107842‑108131
                                                                                                                                                     |                                                                                                                                             
cACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTc                                                                                                                                               <  5:175919/149‑1 (MQ=255)
cACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGAGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTc                                                                                                                                               <  5:170962/149‑1 (MQ=255)
     ccAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGTTGATTGGCGAAAGCTTTGCCaa                                                                                                                                          <  5:326596/149‑1 (MQ=255)
              ccACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACtg                                                                                                                                   <  1:106362/147‑1 (MQ=255)
              ccACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGATTGGCGAAAGCTTTGCCAACCa                                                                                                                                       <  6:332733/143‑1 (MQ=255)
              ccACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGATTGGCGAAAGCTTTGCCAACCa                                                                                                                                       >  5:332733/1‑143 (MQ=255)
               cACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACtggtg                                                                                                                                <  6:154482/149‑1 (MQ=255)
                 cgccgcGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTggg                                                                                                                              <  2:181454/149‑1 (MQ=255)
                                           gCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCaa                                                                                                                                      >  6:44554/1‑115 (MQ=255)
                                           gCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCaa                                                                                                                                      <  5:44554/115‑1 (MQ=255)
                                                tctcGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCg                                                                                                            <  8:12382/136‑1 (MQ=255)
                                                tctcGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCg                                                                                                            >  7:12382/1‑136 (MQ=255)
                                                  tcGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGa                                                                                             >  7:204334/1‑149 (MQ=255)
                                                    gATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCAACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGAt                                                                                            >  1:186979/1‑148 (MQ=255)
                                                            tCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgctg                                                                                   >  7:242261/1‑149 (MQ=255)
                                                              cACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgctgct                                                                                 >  8:135586/1‑149 (MQ=255)
                                                                ccGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAATGGctgctg                                                                                >  2:61216/1‑148 (MQ=255)
                                                                ccGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGgctgcttg                                                                               <  8:204334/149‑3 (MQ=255)
                                                                ccGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGAGCACCACCTCATCCAGCTGGTTGGCAAAGCTATCGCCAATCAACTGGTGGGCAATCAGGTTGGAGACCGTGGTGTCGACCGGTATCAATGGatgctg                                                                                >  2:276870/1‑148 (MQ=255)
                                                                   aTTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGtt                                                                                                                    >  3:161325/1‑109 (MQ=255)
                                                                   aTTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGtt                                                                                                                    <  4:161325/109‑1 (MQ=255)
                                                                    ttCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTTGTGTCGACCGGGATCAATGGCTGCTGGAAGt                                                                           <  1:276870/149‑1 (MQ=255)
                                                                    ttCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAATGGCTGCTGGAAGt                                                                           <  1:61216/149‑1 (MQ=255)
                                                                     tCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTg                                                                          >  7:74436/1‑149 (MQ=255)
                                                                      cAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTgg                                                                         >  7:131590/1‑149 (MQ=255)
                                                                       aaCGGCAAAGCTCTCATCATTAAGCATTTCGTCAACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGc                                                                        <  2:186979/149‑1 (MQ=255)
                                                                                   ctcATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGt                                                            >  3:208667/1‑149 (MQ=255)
                                                                                    tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGt                                                            >  6:41808/1‑148 (MQ=255)
                                                                                    tcatcaTTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGt                                                            >  7:206717/1‑148 (MQ=255)
                                                                                     catcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCa                                                          >  4:105552/1‑149 (MQ=255)
                                                                                      atcatTAAGCATTTGGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGtt                                                                                                                    >  7:299600/1‑90 (MQ=255)
                                                                                      atcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGtt                                                                                                                    <  8:299600/90‑1 (MQ=255)
                                                                                      atcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCAt                                                         >  7:26892/1‑149 (MQ=255)
                                                                                       tcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCAtt                                                        <  8:26892/149‑1 (MQ=255)
                                                                                       tcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCAtt                                                        <  1:241743/149‑1 (MQ=255)
                                                                                       tcatTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGAGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCAt                                                         >  3:306841/1‑148 (MQ=255)
                                                                                          ttAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCgtg                                                                                                          <  2:302898/96‑1 (MQ=255)
                                                                                          ttAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCgtg                                                                                                          >  1:302898/1‑96 (MQ=255)
                                                                                             aGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACct                                                   <  8:242261/148‑1 (MQ=255)
                                                                                                    cgtcgACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTccac                                            >  2:345269/1‑148 (MQ=255)
                                                                                                     gtcgACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGCCCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTccac                                            >  1:250985/1‑147 (MQ=255)
                                                                                                               tggtgCACCACCTCATCCAGCTGGTTGACAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTAc                                                                >  8:164693/1‑117 (MQ=255)
                                                                                                               tggtgCACCACCTCATCCAGCTGGTTGACAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTAc                                                                <  7:164693/117‑1 (MQ=255)
                                                                                                                 gtgCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGAt                                                                    <  3:325668/111‑1 (MQ=255)
                                                                                                                 gtgCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGAt                                                                    >  4:325668/1‑111 (MQ=255)
                                                                                                                  tgCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAg                             <  8:131590/149‑1 (MQ=255)
                                                                                                                  tgCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAg                             <  8:206717/149‑1 (MQ=255)
                                                                                                                        accTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGc                       <  7:135586/149‑1 (MQ=255)
                                                                                                                           tCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAgg                            <  4:237101/141‑1 (MQ=255)
                                                                                                                           tCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAgg                            >  3:237101/1‑141 (MQ=255)
                                                                                                                            ccaccAGTTGATTGGCGAAAGCTTTGCCAACCAGCTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCa                   >  3:166125/4‑149 (MQ=38)
                                                                                                                             aTCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGc                                                                        <  7:140559/95‑1 (MQ=255)
                                                                                                                             aTCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGc                                                                        >  8:140559/1‑95 (MQ=255)
                                                                                                                             aTCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGAGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCa                   <  4:306841/148‑1 (MQ=255)
                                                                                                                                cAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATcc                <  2:250985/148‑1 (MQ=255)
                                                                                                                                 aGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATc                 >  3:281479/1‑146 (MQ=255)
                                                                                                                                      gTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAg           <  4:281479/147‑1 (MQ=255)
                                                                                                                                              aaGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAg  <  2:7392/149‑1 (MQ=255)
                                                                                                                                              aaGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACAGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTAATACAGAAGGCAGTCCAg  >  5:287303/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CACGGCCAGCGCTTCCACGCCGCGCCAGTTGCCGGGCTGCTCAGCTTCTCTCGATAAGCCTCCACCGATTCAACGGCAAAGCTCTCATCATTAAGCATTTCGTCGACCGAGTGGTGCACCACCTCATCCAGCTGGTTGGCAAAGCTTTCGCCAATCAACTGGTGGGCAATCAGGTTGGCGACCGTGGTGTCGACCGGGATCAACGGCTGCTGGAAGTGGCGGATGTACAAGTCATTCACCTCTTCCACCAGGCGGTGCGCCAGGTAGGCTTCATCCAGAAGGCAGTCCAG  >  NC_002947/107842‑108131

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: