Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,388,043 G→T 12.5% G151G (GGC→GGA PP_1209 ← cold‑shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,388,0430GT12.5% 127.9 / 7.2 48G151G (GGC→GGAPP_1209cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base G (23/19);  new base T (3/3);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

GTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCAC  >  NC_002947/1387897‑1388189
                                                                                                                                                    |                                                                                                                                                  
gTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGACCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGccc                                                                                                                                                 <  1:183915/150‑1 (MQ=255)
       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)
                                                                                                                                                    |                                                                                                                                                  
GTTGGTCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCAC  >  NC_002947/1387897‑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: