Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,388,117 C→T 7.8% G127S (GGC→AGC)  PP_1209 ← cold‑shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,388,1170CT 7.8% 144.5 / 3.0 51G127S (GGC→AGC) PP_1209cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base T (2/2);  total (28/23)
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.08e-01

TCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGTCAGGGTCTGCA  >  NC_002947/1387969‑1388266
                                                                                                                                                    |                                                                                                                                                     
tCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCg                                                                                                                                                      >  1:1075426/1‑150 (MQ=255)
    cgcgCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCAt                                                                                                                                                  <  2:1075426/150‑1 (MQ=255)
             cTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCa                                                                                                                                         >  1:740123/1‑150 (MQ=255)
              ttCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCAc                                                                                                                                        <  1:258297/150‑1 (MQ=255)
                 gACCAGAATGCGGTGGCCATCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACgtg                                                                                                                                     >  2:498538/1‑150 (MQ=255)
                   ccAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTTTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGtt                                                                                                                                   >  1:328006/1‑150 (MQ=255)
                   ccAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGtt                                                                                                                                   >  1:328010/1‑150 (MQ=255)
                            ggTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCcgc                                                                                                                          <  2:629602/150‑1 (MQ=255)
                                  ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCTGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCt                                                                                                                    >  2:813112/1‑150 (MQ=255)
                                  ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCTGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCt                                                                                                                    >  2:813123/1‑150 (MQ=255)
                                  ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCt                                                                                                                    >  1:22231/1‑150 (MQ=255)
                                  ccTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCt                                                                                                                    >  1:700945/1‑150 (MQ=255)
                                     cACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                                                                                                 <  2:22231/150‑1 (MQ=255)
                                                 cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGa                                                                                                                                 <  2:954384/122‑1 (MQ=255)
                                                 cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGa                                                                                                                                 >  1:954384/1‑122 (MQ=255)
                                                 cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                                                                                                     >  2:172047/1‑150 (MQ=255)
                                                 cccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                                                                                                     <  1:205690/150‑1 (MQ=255)
                                                  ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGc                                                                                                    >  2:1090434/1‑150 (MQ=255)
                                                           cACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCAc                                                                                           <  1:787652/150‑1 (MQ=255)
                                                            aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGcc                                                                                                                            >  2:422105/1‑116 (MQ=255)
                                                            aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGcc                                                                                                                            <  1:422105/116‑1 (MQ=255)
                                                            aCGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACg                                                                                          >  2:665816/1‑150 (MQ=255)
                                                                     tCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                                                                                                     <  1:1057377/130‑1 (MQ=255)
                                                                     tCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                                                                                                     >  2:1057377/1‑130 (MQ=255)
                                                                       aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc                                                                                        <  1:360232/141‑1 (MQ=255)
                                                                       aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc                                                                                        >  2:360232/1‑141 (MQ=255)
                                                                        tCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacg                                                                              <  1:521200/150‑1 (MQ=255)
                                                                         cGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacgg                                                                             <  2:740123/150‑1 (MQ=255)
                                                                          gCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggc                                                                            >  1:13937/1‑150 (MQ=255)
                                                                              gAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggccacg                                                                        <  1:1090434/150‑1 (MQ=255)
                                                                                 tctcGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGCGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCacggccacggcc                                                                     >  1:733503/1‑150 (MQ=255)
                                                                                   tagggAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGACCAg                                                                   >  1:43271/3‑150 (MQ=255)
                                                                                     gggAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAgcg                                                                 >  1:1020906/1‑150 (MQ=255)
                                                                                           tGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGcac                                                           >  2:39418/1‑150 (MQ=255)
                                                                                                 cAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCa                                                     <  2:13937/150‑1 (MQ=255)
                                                                                                  aaaGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCTGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAg                                                    <  1:813123/150‑1 (MQ=255)
                                                                                                  aaaGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCTGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAg                                                    <  1:813112/150‑1 (MQ=255)
                                                                                                  aaaGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAg                                                    >  2:434849/1‑150 (MQ=255)
                                                                                                   aaGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAgg                                                   <  1:536410/150‑1 (MQ=255)
                                                                                                      ccTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggccacggc                                                                      >  2:990468/1‑128 (MQ=255)
                                                                                                      ccTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggccacggc                                                                      <  1:990468/128‑1 (MQ=255)
                                                                                                             aCGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGt                                         <  1:792531/150‑1 (MQ=255)
                                                                                                                  ttGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGccc                                    >  2:490698/1‑150 (MQ=255)
                                                                                                                    gAACCAGTTGACGGTGCCGGGAGCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCAt                                  >  2:742704/1‑150 (MQ=255)
                                                                                                                      aCCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGAACCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCAttt                                <  1:851625/150‑1 (MQ=255)
                                                                                                                          tttGACGGTGCCGGTATCACGGTTGCCGGCATCTTGCGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCggg                            <  2:733503/150‑1 (MQ=255)
                                                                                                                                          tCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGcaca                                                          >  2:39292/1‑104 (MQ=255)
                                                                                                                                          tCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGcaca                                                          <  1:39292/104‑1 (MQ=255)
                                                                                                                                            aCGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGGCAg          >  1:205428/1‑150 (MQ=255)
                                                                                                                                             cgatctccGGCATCTTGTCCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACCACCGGGCGGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGTCAgg         <  2:588215/144‑1 (MQ=255)
                                                                                                                                             cGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGTCAgg         <  2:1056534/150‑1 (MQ=255)
                                                                                                                                              ggTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGGCAggg        >  2:649279/1‑150 (MQ=255)
                                                                                                                                                    cGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGTCAGGGTCTgca  >  2:810428/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCACGCGCTGGCCTTCGACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGCCACGGCCAGCGCCGCACACACCAGGGCTGCCAGGTTGCCCATTTCGGGGCGGGCCAGTAGCGTCAGGGTCTGCA  >  NC_002947/1387969‑1388266

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