Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,178,240 A→T 7.0% Q159L (CAG→CTG)  PP_1931 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,178,2400AT 7.0% 169.7 / 2.9 57Q159L (CAG→CTG) PP_1931hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base T (2/2);  total (28/29)
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.39e-01

CAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGCTGAAACAA  >  NC_002947/2178097‑2178385
                                                                                                                                               |                                                                                                                                                 
cAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATAGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATAGCGACCGTAACTTCACCCAGGcc                                                                                                                                               >  5:98829/1‑148 (MQ=255)
 aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCaa                                                                                                                                             <  5:49158/149‑1 (MQ=255)
 aaGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCaa                                                                                                                                             <  5:43847/149‑1 (MQ=255)
        aCGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAg                                                                                                                                      <  6:140916/149‑1 (MQ=255)
        aCGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAg                                                                                                                                      <  7:126589/149‑1 (MQ=255)
          ggTAAAGCCATTCGTCGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCTTGGCCTGGGTGAAgtt                                                                                                                                    <  1:285721/149‑3 (MQ=255)
          ggTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGc                                                                                                                                    >  1:26036/1‑149 (MQ=255)
                   aTTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCTCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATgca                                                                                                                           <  7:50723/149‑1 (MQ=255)
                             tCTGCCGTCAGTTCTTGATCCTCTGCCGTAACCTAAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgag                                                                                                                 <  2:142447/149‑1 (MQ=255)
                                 ccGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcaga                                                                                                             <  6:123947/149‑1 (MQ=255)
                                 ccGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagca                                                                                                               >  5:378461/1‑147 (MQ=255)
                                   gTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcagag                                                                                                            >  1:139106/1‑148 (MQ=255)
                                    tCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcagagca                                                                                                          >  3:95911/1‑149 (MQ=255)
                                    tCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCgagcagagca                                                                                                          >  4:66476/1‑149 (MQ=255)
                                      aGTTCGTGATCCTCTGCCGTAACGTCAACCTTTTCTCTAAGTCGATCATTGCCATCGACGGCAGCATGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCAtt                                                                                                        <  7:156531/149‑1 (MQ=255)
                                         tCGTGATCCTCTGCCTTAACCTCAAACTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGa                                                                                                      <  8:326857/148‑1 (MQ=255)
                                          cgggctCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTCAACAAACGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAg                                                                                                     <  2:357486/143‑1 (MQ=255)
                                              aTCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCt                                                                                                 <  2:26036/148‑1 (MQ=255)
                                                            ccTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATgg                                                                                  <  6:378461/149‑1 (MQ=255)
                                                              tCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATg                                                                                   <  6:172300/146‑1 (MQ=255)
                                                              tCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATg                                                                                   >  5:172300/1‑146 (MQ=255)
                                                                  ccTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg                                                                            <  2:139106/149‑1 (MQ=255)
                                                                   cTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATc                                                                                                                    >  6:219919/1‑108 (MQ=255)
                                                                   cTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATc                                                                                                                    <  5:219919/108‑1 (MQ=255)
                                                                    ttttCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg                                                                            >  5:139415/1‑147 (MQ=255)
                                                                        ctctGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg                                                                            >  8:65596/1‑143 (MQ=255)
                                                                        ctctGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCgg                                                                            <  7:65596/143‑1 (MQ=255)
                                                                           tGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCggg                                                                   >  3:325303/1‑149 (MQ=255)
                                                                             aGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCgg                                                                    <  3:66476/146‑1 (MQ=255)
                                                                               tCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAAt                                                               >  1:15185/1‑149 (MQ=255)
                                                                                 gATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGc                                                             >  5:28296/1‑149 (MQ=255)
                                                                                  aTCATTGCCATCGACGGCAGCAAGATCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGATACCTTACAGCGATGGATTCGGCAGACCGGGAAATGcc                                                            >  5:58705/1‑149 (MQ=255)
                                                                                   tCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCg                                                           >  7:215589/1‑149 (MQ=255)
                                                                                         ccATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTg                                                     <  4:325303/149‑1 (MQ=255)
                                                                                         ccATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTg                                                     >  7:148982/1‑149 (MQ=255)
                                                                                         ccATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGAAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCATGAGGtt                                                     >  3:80161/1‑148 (MQ=255)
                                                                                           aTCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCg                                                                             >  1:169588/1‑123 (MQ=255)
                                                                                           aTCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCg                                                                             <  2:169588/123‑1 (MQ=255)
                                                                                           aTCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGc                                                   >  7:326216/1‑149 (MQ=255)
                                                                                              gACGGCAGCAAGTTGAAACCAGTGTACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGGCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGTTGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCAt                                                <  6:98829/149‑1 (MQ=255)
                                                                                                    agcaAGATAAAAGCATTTAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGaa                                          <  6:58705/149‑1 (MQ=255)
                                                                                                        aGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGcc                                      <  8:63105/149‑1 (MQ=255)
                                                                                                                      aaCAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAaga                        <  1:29868/149‑1 (MQ=255)
                                                                                                                       aCAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGaaa                          >  7:353735/1‑146 (MQ=255)
                                                                                                                       aCAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGaaa                          >  4:232621/1‑146 (MQ=255)
                                                                                                                       aCAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAaga                        <  8:353735/148‑1 (MQ=255)
                                                                                                                           tCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAagag                       <  1:188358/145‑1 (MQ=255)
                                                                                                                           tCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAagag                       >  2:188358/1‑145 (MQ=255)
                                                                                                                           tCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAagag                       >  1:128820/1‑145 (MQ=255)
                                                                                                                                  cGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAAc            >  6:239810/1‑149 (MQ=255)
                                                                                                                                    tAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGc          <  8:148982/149‑1 (MQ=255)
                                                                                                                                       cTTCACCTTGGCCTGGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGc                                  >  3:298521/1‑122 (MQ=255)
                                                                                                                                       cTTCACCTTGGCCTGGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGc                                  <  4:298521/122‑1 (MQ=255)
                                                                                                                                       cTTCACCTTGGCCTGGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGCTg        >  6:332622/1‑148 (MQ=255)
                                                                                                                                       cTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGCTg        >  5:141759/1‑148 (MQ=255)
                                                                                                                                            cccAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGCTGAaacaa  <  6:28296/149‑1 (MQ=255)
                                                                                                                                             ccAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGCTGAaacaa  >  1:151268/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CAAGGACAACGGTAAAGCCATTCGTAGCGTCTGCCGTCAGTTCGTGATCCTCTGCCGTAACCTCAACCTTTTCTCTGAGTCGATCATTGCCATCGACGGCAGCAAGTTCAAAGCAGTGAACAATCGCGACCGTAACTTCACCCAGGCCAAGGTGAAGGCCCGCATGCAGCAAATCGAGCAGAGCATTGAGCGCTACCTTACAGCGATGGATTCGGCAGACCGGGCAATGCCGGAGGTGGCCATGGCGAAGGCCGAGCGGCTGAAAGAGAAGATCGAAACGCTGAAACAA  >  NC_002947/2178097‑2178385

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