Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,800,129 T→G 8.2% L155R (CTG→CGG)  PP_3360 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,800,1290TG 8.2% 147.3 / 3.7 49L155R (CTG→CGG) PP_3360membrane protein
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base G (2/2);  total (24/25)
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.68e-01

TGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGTGGACAAGCCGTTC  >  NC_002947/3799986‑3800271
                                                                                                                                               |                                                                                                                                              
tccccTGAGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTACGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGAGGTCAGTCGTAATGCTGCCGa                                                                                                                                           >  4:57847/4‑149 (MQ=255)
  gCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCGGGCGATg                                                                                                                                         >  4:297035/1‑149 (MQ=255)
  gCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCGGGCGATg                                                                                                                                         <  3:297035/149‑1 (MQ=255)
  gCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGCGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCGGGCGATg                                                                                                                                         >  4:302042/1‑149 (MQ=255)
  cccTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTATCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCGGGCGATg                                                                                                                                         <  3:302042/148‑1 (MQ=255)
    cTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCt                                                                                                                                       <  4:144581/149‑1 (MQ=255)
     atctggCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGGCAGTCGTAGTGCTGGCGATGCTg                                                                                                                                      <  7:326451/145‑1 (MQ=255)
       tgGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTc                                                                                                                                    <  4:255116/149‑1 (MQ=255)
          gCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTg                                                                                                                                      >  8:326451/1‑144 (MQ=255)
               ccGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTgg                                                                                                                            <  8:50235/149‑1 (MQ=255)
                cgcgCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTggg                                                                                                                           <  8:319493/149‑1 (MQ=255)
                 gcgcCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTgggg                                                                                                                          <  5:359828/149‑1 (MQ=255)
                  cgcCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTCGCAACCCTTTCAGCCTGGGTCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTggggg                                                                                                                         >  8:325917/1‑149 (MQ=255)
                    ccATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTg                                                                                                                       <  2:46125/149‑1 (MQ=255)
                    ccATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTg                                                                                                                       <  8:259768/149‑1 (MQ=255)
                                         gTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTTGTGCTGGCGATGCTGTCGAACCTgggg                                                                                                                          >  5:188519/1‑125 (MQ=255)
                                         gTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTgggg                                                                                                                          <  6:188519/125‑1 (MQ=255)
                                                      ccACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCtg                                                                                            <  7:72696/142‑1 (MQ=255)
                                                      ccACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCtg                                                                                            >  8:72696/1‑142 (MQ=255)
                                                                gACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGc                                                                           <  3:57847/149‑1 (MQ=255)
                                                                    gCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGaa                                                                                                                     >  6:307704/1‑103 (MQ=255)
                                                                      ctcaaCCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCTCCAGCCtgggcgtgg                                                                                     <  8:191920/130‑1 (MQ=255)
                                                                            cAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCgctggc                                                               >  5:192937/1‑149 (MQ=255)
                                                                                     gcACTGGCAACCCTTTCAGACTGGTGCGCGCGTGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGaa                                                      <  6:192937/149‑1 (MQ=255)
                                                                                      cACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCgct                                                                  <  5:136359/136‑1 (MQ=255)
                                                                                      cACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCgct                                                                  >  6:136359/1‑136 (MQ=255)
                                                                                          ggCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATcc                                                 >  8:210748/1‑149 (MQ=255)
                                                                                          ggCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATcc                                                 >  7:316371/1‑149 (MQ=255)
                                                                                           gCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTtc                                                                                                        <  1:165749/93‑1 (MQ=255)
                                                                                           gCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTtc                                                                                                        >  2:165749/1‑93 (MQ=255)
                                                                                           gCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGcc                                                                                                  >  6:124694/1‑99 (MQ=255)
                                                                                           gCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGcc                                                                                                  <  5:124694/99‑1 (MQ=255)
                                                                                           gCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGa                                                       <  5:52333/142‑1 (MQ=255)
                                                                                           gCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGa                                                       >  6:52333/1‑142 (MQ=255)
                                                                                             accccTTTCAGCGTGGGGGGCGGGGGTCGGGTTGGGGTGAGTCGTAGTGCTGGCGATGCTGGCGAAGCTGGGGGTGAAGATCCGCGCGTTCGTCGCCCGCCTGGGCGGGGGCGGCATCGCAGTGGCGCTGGCGGTGCTGAATATGCTTg                                              >  6:172068/3‑149 (MQ=255)
                                                                                              aCCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTgg                                             <  5:140096/149‑1 (MQ=255)
                                                                                                  tttCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGAc                                         >  7:206908/1‑149 (MQ=255)
                                                                                                      aGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGt                                     >  6:155375/1‑149 (MQ=255)
                                                                                                       gCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGtt                                    >  5:362279/1‑149 (MQ=255)
                                                                                                        ccTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTc                                   <  7:210748/149‑1 (MQ=255)
                                                                                                            gggCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGc                                 <  7:325917/147‑1 (MQ=255)
                                                                                                                          tgatgtgGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGTCGTGGGCGGCATCGCAGTGTCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTACGTTGTCGATCGcc                 >  4:13191/4‑149 (MQ=255)
                                                                                                                               tgGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGt               >  2:215834/1‑146 (MQ=255)
                                                                                                                               tgGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGt               <  1:215834/146‑1 (MQ=255)
                                                                                                                                    aGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGTGGACAAGc       <  8:316371/149‑1 (MQ=255)
                                                                                                                                         tAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGTGGACAAGCCGTTc  >  1:14018/1‑149 (MQ=255)
                                                                                                                                              cTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTc                        <  8:196015/122‑1 (MQ=255)
                                                                                                                                              cTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTc                        >  7:196015/1‑122 (MQ=255)
                                                                                                                                              cTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCATCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATACTTCGCGACCTCTTCGCTTCGTTGTCGATCGCCGTGGACAAGCCGtta  >  5:155853/1‑143 (MQ=255)
                                                                                                                                               tGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGTGGACAAGCCg     >  4:126577/1‑140 (MQ=255)
                                                                                                                                               tGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGTGGACAAGCCg     <  3:126577/140‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGGCCTGTGGGCCAACCGCGCCATCGGCCTGGGCCTGAGCCGTTACTTCGCCCGCCACCGTACCGACGGCCTGAACCAGGGCAGCGCACTGGCAACCCTTTCAGCCTGGGGCGCGCGGGTCCTGTTGTGGTCAGTCGTAGTGCTGGCGATGCTGTCGAACCTGGGGGTGAACATCACCGCCTTCGTCGCCAGCCTGGGCGTGGGCGGCATCGCAGTGGCGCTGGCTGTGCAGAATATCCTTGGCGACCTGTTCGCTTCGTTGTCGATCGCCGTGGACAAGCCGTTC  >  NC_002947/3799986‑3800271

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