Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 1,051,638 C→T 9.4% R60C (CGT→TGT)  cspG → cold shock protein homolog, cold‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,051,6380CT 9.4% 118.7 / 2.5 42R60C (CGT→TGT) cspGcold shock protein homolog, cold‑inducible
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (2/2);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTTAACGCAGATAAAGGTTTTGGCTTTATCACTCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAA  >  NC_000913_3_pae_tpiA/1051494‑1051781
                                                                                                                                                |                                                                                                                                               
tttAACGCAGATAAAGGTTTTGGCTTTATCACTCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCAATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTACTGCGCAGGGGAAACCTga                                                                                                                                             >  2:376094/1‑148 (MQ=255)
  tAACGCAGATAAAGGTTTTGGCTTTATCACTCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGcc                                                                                                                                           >  1:769208/1‑149 (MQ=255)
          aTAAAGGTTTTGGCTTTATCACTCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCa                                                                                                                                   >  1:300744/1‑149 (MQ=255)
                               cTCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCCCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAATTGATTTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTaa                                                                                                              <  1:321683/149‑1 (MQ=255)
                                tCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTacg                                                                                                             >  2:890657/1‑147 (MQ=255)
                                tCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAg                                                                                                             >  2:890648/1‑149 (MQ=255)
                                tCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCGATGGAGCAGGGCCAACGTGGCCCCGGGGCGGCGAACTGTGTTACGGTCTact                                                                                                             >  1:96188/1‑147 (MQ=255)
                                 ccTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAgg                                                                                                            >  2:174268/1‑149 (MQ=255)
                                     atgatgGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGc                                                                                                        <  1:552660/149‑1 (MQ=255)
                                       gatgGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCa                                                                                                      <  1:566635/149‑1 (MQ=255)
                                        atgGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCat                                                                                                     <  2:175075/149‑1 (MQ=255)
                                          ggCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCatta                                                                                                   >  2:61970/1‑149 (MQ=255)
                                             agcaAAGACGTTTTCGTCCATTTCACCTCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTGTCTATTGAGCAAGAGTAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCAttatta                                                                                                <  1:881652/149‑1 (MQ=255)
                                                aaaGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTc                                                                                             <  1:407911/149‑1 (MQ=255)
                                                      gTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACAtc                                                                                       >  1:589516/1‑149 (MQ=255)
                                                      gTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATAAGATCTGTGAATTTTATATTTAGAAGGGTCAACGTGGCGCCGCGGCAGGGGACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACAtc                                                                                       >  1:225298/1‑149 (MQ=255)
                                                        tttCGTCCATTTCACCGCCATCCAGAGCAATTAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGTCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTGGCCATTATTACTCAACAtctc                                                                                     <  2:209869/149‑1 (MQ=255)
                                                         ttCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTcc                                                                                    <  1:805501/149‑1 (MQ=255)
                                                             tCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCAttt                                                                                >  1:501134/1‑149 (MQ=255)
                                                                  ttCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCt                                                                           <  2:389051/149‑1 (MQ=255)
                                                                  ttCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGGTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTTTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCt                                                                           <  2:673513/149‑1 (MQ=255)
                                                                      ccGCCATCCAGAGCAATGAATTCCGCACGCTTAACGAAATTCAGAAAGTTGATTTGTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCATCGAACTTTTTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTcc                                                                       <  1:598319/149‑1 (MQ=255)
                                                                            tCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAATGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCAtgttg                                                                 >  1:816037/1‑149 (MQ=255)
                                                                            tCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAATGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCAtgttg                                                                 >  1:294627/1‑149 (MQ=255)
                                                                            tCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCAtgttg                                                                 >  1:108494/1‑149 (MQ=255)
                                                                               agagCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATTTTGTCa                                                              >  2:936415/1‑149 (MQ=255)
                                                                                    aaTGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGtt                                                         >  1:368312/1‑149 (MQ=255)
                                                                                    aaTGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGtt                                                         <  1:502798/149‑1 (MQ=255)
                                                                                     aTGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTGATAGTTc                                                        >  2:43460/1‑149 (MQ=255)
                                                                                      tGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTca                                                       >  1:44398/1‑149 (MQ=255)
                                                                                        aaTTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTcaca                                                     <  1:893709/149‑1 (MQ=255)
                                                                                           tCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTa                                                  <  2:380295/149‑1 (MQ=255)
                                                                                             cGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTAcc                                                >  2:783306/1‑149 (MQ=255)
                                                                                                         gAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCCTGTTGTCATGGTTCACAGTACCGCACATCGGCAt                                    >  1:414248/1‑149 (MQ=255)
                                                                                                            aaTCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATAGGCAtcat                                 >  2:698232/1‑146 (MQ=255)
                                                                                                                  aaaGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACACCTCCATTTCCGCTGTCCATGTTTTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGa                           <  1:642800/149‑1 (MQ=255)
                                                                                                                      ttGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACGCCACATCGGCATTCGATGTGACGGa                       >  2:72430/1‑149 (MQ=255)
                                                                                                                      ttGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGa                       >  1:829247/1‑149 (MQ=255)
                                                                                                                                     aGCAGGGGCAATGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTggg        <  2:816037/149‑1 (MQ=255)
                                                                                                                                     aGCAGGGGCAATGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTggg        <  2:294627/149‑1 (MQ=255)
                                                                                                                                        aGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACGGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGgcgc     >  1:764056/1‑149 (MQ=255)
                                                                                                                                           tgCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGCTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTaa  <  2:300744/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TTTAACGCAGATAAAGGTTTTGGCTTTATCACTCCTGATGATGGCAGCAAAGACGTTTTCGTCCATTTCACCGCCATCCAGAGCAATGAATTCCGCACGCTGAACGAAAATCAGAAAGTTGAATTTTCTATTGAGCAGGGGCAACGTGGCCCCGCGGCAGCGAACGTTGTTACGCTCTAAGGTTGCCATTATTACTCAACATCTCCATTTCCGCTGTCCATGTTGTCATGGTTCACAGTACCGCACATCGGCATTCGATGTGACGGAGCGAAACCCTTTGGGCGCTAA  >  NC_000913_3_pae_tpiA/1051494‑1051781

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