Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,388,050 T→A 13.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,0500TA13.3% 115.8 / 10.0 46D149V (GAC→GTC) PP_1209cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base A (3/3);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

AACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCAC  >  NC_002947/1387904‑1388189
                                                                                                                                                    |                                                                                                                                           
aaCGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGtctc                                                                                                                                          >  1:527141/1‑150 (MQ=255)
   ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCggg                                                                                                                                       >  1:164077/1‑150 (MQ=255)
   ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCggg                                                                                                                                       >  2:123789/1‑150 (MQ=255)
      ggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCTGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGgaga                                                                                                                                    <  1:269323/150‑1 (MQ=255)
             tGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAAcc                                                                                                                             >  2:118556/1‑150 (MQ=255)
                  ccACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGaaa                                                                                                                               <  1:596017/143‑1 (MQ=255)
                  caaCGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGaaa                                                                                                                               >  2:596017/3‑143 (MQ=255)
                   cACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGcc                                                                                                                       >  1:82757/1‑150 (MQ=255)
                     cGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCtt                                                                                                                     >  2:745855/1‑150 (MQ=255)
                     cGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGAATTGACTGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGGCTCGGGAGGTGAAACCAAAGCCtt                                                                                                                     >  1:535817/1‑150 (MQ=255)
                      gTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                    <  1:348214/150‑1 (MQ=255)
                       tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCCGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                   <  1:572898/150‑1 (MQ=255)
                         ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                   <  1:381232/132‑1 (MQ=255)
                         ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                   >  2:381232/1‑132 (MQ=255)
                         ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGa                                                                                                                 >  2:119874/1‑150 (MQ=255)
                               cTGCAGGCCCTTGTCGCGGTGCCCACCGTAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATTGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                           >  1:706763/1‑150 (MQ=255)
                                     gCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCa                                                                                                     <  2:164077/150‑1 (MQ=255)
                                           gTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTgcc                                                                                               <  2:610390/150‑3 (MQ=255)
                                                   gCATCACCGAAAACTCCCCGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTAGCGAGTTGTGGTCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCgg                                                                                       <  2:706763/150‑1 (MQ=255)
                                                        aCCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCa                                                                                  <  1:123789/150‑1 (MQ=255)
                                                          cGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCgggcga                                                                                                                                    >  2:387937/1‑95 (MQ=255)
                                                          cGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGgcga                                                                                                                                    <  1:387937/98‑1 (MQ=255)
                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCa                                                                                  <  1:388998/141‑1 (MQ=255)
                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCa                                                                                  >  2:388998/1‑141 (MQ=255)
                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCg                                                                         >  1:97597/1‑150 (MQ=255)
                                                                  tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCg                                                                         >  1:357149/1‑150 (MQ=255)
                                                                         gcTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                   >  1:726470/1‑85 (MQ=39)
                                                                         gcTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                   <  2:726470/85‑1 (MQ=39)
                                                                           tGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                                <  2:82757/150‑1 (MQ=255)
                                                                           tGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGAGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTtgt                                                                <  2:535817/150‑1 (MQ=255)
                                                                                ttCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCAc                                                           >  1:339087/1‑150 (MQ=255)
                                                                                     ccAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGtt                                                      <  1:709716/150‑1 (MQ=255)
                                                                                         aaTGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                  >  2:505636/1‑150 (MQ=255)
                                                                                                     cTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCtt                                      <  1:505636/150‑1 (MQ=255)
                                                                                                            cgGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGtg                               >  1:505857/1‑150 (MQ=255)
                                                                                                            cgGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGtg                               <  2:505857/150‑1 (MQ=255)
                                                                                                               gATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGtgct                             <  1:118556/150‑1 (MQ=255)
                                                                                                                    cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                    <  2:32066/58‑1 (MQ=255)
                                                                                                                    cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                    >  1:32066/1‑58 (MQ=255)
                                                                                                                     ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGc                       >  2:405444/1‑150 (MQ=255)
                                                                                                                        aaaGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                    <  1:244934/134‑1 (MQ=255)
                                                                                                                        aaaGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                    >  2:244934/1‑134 (MQ=255)
                                                                                                                               aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCtt                                                                  <  2:234647/97‑1 (MQ=255)
                                                                                                                               aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCtt                                                                  >  1:234647/1‑97 (MQ=255)
                                                                                                                                 gAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc           >  2:256277/1‑150 (MQ=255)
                                                                                                                                  aaGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGcc          >  2:242437/1‑150 (MQ=255)
                                                                                                                                          aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcac  >  2:413851/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
AACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCAC  >  NC_002947/1387904‑1388189

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