Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,225,015 T→C 7.9% I540V (ATT→GTT)  PP_1963 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,225,0150TC 7.9% 143.8 / 2.2 52I540V (ATT→GTT) PP_1963hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/25);  new base C (2/2);  total (25/27)
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.87e-01

TTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCATACGGTCCGCCATTGCAGCGCCTG  >  NC_002947/2224872‑2225158
                                                                                                                                               |                                                                                                                                               
ttCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCa                                                                                                                                             >  1:209338/1‑148 (MQ=255)
   gTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCa                                                                                                                                             >  1:173219/1‑145 (MQ=255)
        ttCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCCTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTc                                                                                                                                    <  3:39916/149‑1 (MQ=255)
                 gggAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCtgctgc                                                                                                                           >  6:122647/1‑149 (MQ=255)
                 gggAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAACGGCAAAGATTTTCACCtgctgc                                                                                                                           >  4:144596/1‑149 (MQ=255)
                 gggAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAACGGCAAAGATTTTCACCtgctgc                                                                                                                           >  4:149460/1‑149 (MQ=255)
                          cgctGCCCTTGTCTATGCAGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGa                                                                                                                  <  8:162336/146‑1 (MQ=255)
                           gATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGAt                                                                                                                 <  6:71975/149‑1 (MQ=255)
                             tGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCtgctttgc                                                                                                                      >  6:120324/1‑142 (MQ=255)
                             tGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCtgctttgc                                                                                                                      <  5:120324/142‑1 (MQ=255)
                               cACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCAtt                                                                                                             <  2:57389/149‑1 (MQ=255)
                                aCTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTg                                                                                                            <  4:298996/149‑1 (MQ=255)
                                  ttGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATg                                                                                                                                                >  6:128364/1‑111 (MQ=255)
                                  ttGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATg                                                                                                                                                <  5:128364/111‑1 (MQ=255)
                                    gTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGc                                                                                                          <  2:209338/147‑1 (MQ=255)
                                     tCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACct                                                                                                                                <  4:106206/124‑1 (MQ=255)
                                     tCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACct                                                                                                                                >  3:106206/1‑124 (MQ=255)
                                       tATGCTGTAGCATATCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGACTGGTAGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACAAGGAAAGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCt                                                                                                     >  3:248077/1‑149 (MQ=255)
                                         tGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCtt                                                                                                    <  2:205156/148‑1 (MQ=255)
                                           cTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCttac                                                                                                  >  1:291295/1‑146 (MQ=255)
                                              tAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTACCCaa                                                                                              <  2:291295/149‑1 (MQ=255)
                                                  aTTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCTAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAgg                                                                                          >  8:293196/1‑149 (MQ=255)
                                                   tttctttcCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGc                                                                                         <  5:308285/149‑1 (MQ=255)
                                                       tttcCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGcc                                                                                                 <  4:164377/137‑1 (MQ=255)
                                                       tttcCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGcc                                                                                                 >  3:164377/1‑137 (MQ=255)
                                                           cATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGaa                                                                                  <  5:122647/148‑1 (MQ=255)
                                                            aTTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACt                                                                                >  2:137574/1‑149 (MQ=255)
                                                                 gAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACgg                                                                           >  7:104278/1‑149 (MQ=255)
                                                                 gAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGACCAAGAGGCCCGTGAACTTACgg                                                                           >  2:313688/1‑149 (MQ=255)
                                                                  aaCCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGt                                                                          >  8:172634/1‑149 (MQ=255)
                                                                     cGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTtaa                                                                       <  2:173219/149‑1 (MQ=255)
                                                                           cTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAg                                                                   >  4:196212/1‑147 (MQ=255)
                                                                           cTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTa                                                                 <  7:172634/149‑1 (MQ=255)
                                                                              gAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGaa                                                                                  >  6:224535/1‑129 (MQ=255)
                                                                              gAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGaa                                                                                  <  5:224535/129‑1 (MQ=255)
                                                                                    cTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTaa                                                        >  2:100294/1‑149 (MQ=255)
                                                                                     ttCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAAc                                                       >  5:24131/1‑149 (MQ=255)
                                                                                       ccAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAATCg                                                     >  3:290028/1‑149 (MQ=255)
                                                                                               gTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGAc                                              <  7:293196/148‑1 (MQ=255)
                                                                                                    tACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAACGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCAtt                                        <  3:144596/149‑1 (MQ=255)
                                                                                                    tACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAACGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCAtt                                        <  3:149460/149‑1 (MQ=255)
                                                                                                      cTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGt                                                                          <  8:163956/113‑1 (MQ=255)
                                                                                                      cTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGt                                                                          >  7:163956/1‑113 (MQ=255)
                                                                                                      cTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTcc                                      >  3:186325/1‑149 (MQ=255)
                                                                                                      cTGAATTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTcc                                      <  6:326195/149‑1 (MQ=255)
                                                                                                          ttttCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCa                                   <  2:39701/148‑1 (MQ=255)
                                                                                                                 aTCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATc                           <  1:137574/149‑1 (MQ=255)
                                                                                                                  tCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCa                          <  1:100294/149‑1 (MQ=255)
                                                                                                                    cccAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCATa                        <  6:24131/149‑1 (MQ=255)
                                                                                                                                   gATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCATACGGTCCGCCATTGCa         <  3:196212/149‑1 (MQ=255)
                                                                                                                                           ggAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACGTCATTCCGCAGTCTTATCATACGGTCCGCCATTGCAGCGCCTg  >  6:217631/1‑148 (MQ=255)
                                                                                                                                           ggAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTCAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCATACGCTCCGCCATTGCAGCGCCTg  >  6:212983/1‑148 (MQ=255)
                                                                                                                                               tGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCATACGGTCCGCCATTgta         >  6:270106/1‑135 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TTCGTAAATTCGAACCTGGGAGCAGCCGATGCACTTGTCTATGCTGTAGCATTTCTTTCCATTAGGAACCGATTGCTCGAATCCCTTCCAGTCTGGTTGATACTGATTTTCGCATCCCCAGAGCGGGCCGGGATGACCAGGAATGGCAAAGATTTTCACCTGCTGCTTTGCTTGATCATTGGCACGCTTGCCCAAGAGGCCCGTGAACTTACGGTTAATAAGTAGGGCTGTAACCGCCTCGACCTCATTCCGCAGTCTTATCATACGGTCCGCCATTGCAGCGCCTG  >  NC_002947/2224872‑2225158

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