Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,388,050 T→A 8.3% D149V (GAC→GTC)  PP_1209 ← cold‑shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,388,0500TA 8.3% 136.5 / 4.9 49D149V (GAC→GTC) PP_1209cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base A (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

AACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCA  >  NC_002947/1387904‑1388194
                                                                                                                                                   |                                                                                                                                                
aaCGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGtctc                                                                                                                                               >  1:400630/1‑150 (MQ=255)
       gTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGAt                                                                                                                                        >  2:38875/1‑150 (MQ=255)
             tGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCAGGAAGATGTCATCGCCCGAGGCTCGGgagagatcgga                                                                                                                                  >  2:327048/1‑143 (MQ=255)
             tGCAACCACGTCTTCGGCCTGCAGGCCCTCGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGga                                                                                                                                           >  2:801556/1‑141 (MQ=255)
             tGCAACCACGTCTTCGGCCTGCAGGCCCTCGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGga                                                                                                                                           <  1:801556/141‑1 (MQ=255)
              gCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCa                                                                                                                                 <  1:13516/150‑1 (MQ=255)
              gCAACCACGTCTTCGGCCTGCAGGCCCTTGTACCGGCG‑AACACATGGGCATTCTGGGCTCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCa                                                                                                                                 >  1:150053/1‑150 (MQ=255)
               cAACCACGTCATCGGCCTGCAGGCCCTTGTCGCTGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCAACCAGAATGCGGTGGCCCTCACCGCGGATTTCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCGCGGGCGAGGAAACGaa                                                                                                                                >  1:688398/1‑150 (MQ=255)
                  ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                        <  1:803941/139‑1 (MQ=255)
                  ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                        >  2:803941/1‑139 (MQ=255)
                  ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGc                                                                                                                             >  1:688885/1‑150 (MQ=255)
                     cGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACTCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCtt                                                                                                                          <  2:400630/150‑1 (MQ=255)
                        cTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                        >  1:514333/1‑149 (MQ=255)
                        cTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                        <  2:514333/149‑1 (MQ=255)
                         ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGa                                                                                                                      >  1:38302/1‑150 (MQ=255)
                         ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGa                                                                                                                      <  1:6552/150‑1 (MQ=255)
                               cTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCCGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGgcga                                                                                                                                         <  1:700101/125‑1 (MQ=255)
                               cTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCgggcga                                                                                                                                         >  2:700101/1‑122 (MQ=255)
                                 gCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGa                                                                                                              <  1:749208/150‑1 (MQ=255)
                                     gCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCa                                                                                                          >  1:613629/1‑150 (MQ=255)
                                         ttGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGga                                                                                                                                           <  2:701135/113‑1 (MQ=255)
                                         ttGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGga                                                                                                                                           >  1:701135/1‑113 (MQ=255)
                                              gcgGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGt                                                                                                 <  1:712835/150‑1 (MQ=255)
                                                                 cTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGcc                                                                               <  1:38875/150‑1 (MQ=255)
                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCg                                                                              >  2:808813/1‑150 (MQ=255)
                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCg                                                                              <  2:144630/150‑1 (MQ=255)
                                                                   ccACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCgg                                                                             >  1:396271/1‑150 (MQ=255)
                                                                               cTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCa                                                                 >  1:512350/1‑150 (MQ=255)
                                                                                    aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgtgt                                                            <  2:613629/150‑1 (MQ=255)
                                                                                           tGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGc                                                     <  2:396271/150‑1 (MQ=255)
                                                                                                    ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                                     >  1:60613/1‑125 (MQ=255)
                                                                                                    ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                                     <  2:60613/125‑1 (MQ=255)
                                                                                                    ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCt                                            >  1:84289/1‑150 (MQ=255)
                                                                                                        aCCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                                     <  1:142457/121‑1 (MQ=255)
                                                                                                        aCCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                                     >  2:142457/1‑121 (MQ=255)
                                                                                                         ccGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgcgc                                       >  1:676035/1‑150 (MQ=255)
                                                                                                                tGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGtgctgc                                <  2:38302/150‑1 (MQ=255)
                                                                                                                 ggCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCg                               <  1:808813/150‑1 (MQ=255)
                                                                                                                  gaccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCTGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                              <  2:688885/148‑1 (MQ=255)
                                                                                                                  gCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                              <  2:512350/150‑1 (MQ=255)
                                                                                                                  gCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCgg                              <  1:497901/150‑1 (MQ=255)
                                                                                                                    ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGtt                                                                                  >  1:387587/1‑96 (MQ=255)
                                                                                                                    ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGtt                                                                                  <  2:387587/96‑1 (MQ=255)
                                                                                                                    ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGc                            >  1:753790/1‑150 (MQ=255)
                                                                                                                       aaaGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCaaa                         <  2:150053/150‑1 (MQ=255)
                                                                                                                            gCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCa                    <  1:71483/150‑1 (MQ=255)
                                                                                                                              aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACg                  >  2:209954/1‑150 (MQ=255)
                                                                                                                                         aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcac       >  1:1337/1‑150 (MQ=255)
                                                                                                                                              ccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggcca  >  2:447271/1‑150 (MQ=255)
                                                                                                                                              ccGAGTCTCGGGAGACGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCCAAGCCCACGGCCAGCTGcacggcca  >  2:286895/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCAC‑CGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCA  >  NC_002947/1387904‑1388194

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