Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,638,571 G→A H33Y (CAT→TAT)  cspI ← cold shock‑like protein CspI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,638,5710GA96.2% 85.2 / NA 26H33Y (CAT→TAT) cspIcold shock‑like protein CspI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/13);  minor base C (1/0);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGATTTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTA  >  NC_000913/1638327‑1638803
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
caattccctggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                                                                                                                                                                                                    <  1:330438‑M2/111‑1 (MQ=255)
              tgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctgc                                                                                                                                                                                                                      >  1:13963‑M2/127‑251 (MQ=255)
               gatgccggttacccaaagctatgccggactgaatgtcaaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTACTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGcc                                                                                                                                                                                                                     >  2:473451‑M2/126‑251 (MQ=255)
                      gttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTc                                                                                                                                                                                                              <  1:365864‑M2/133‑1 (MQ=255)
                                          actgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                                                                                                                                                                                          <  1:452069‑M2/153‑1 (MQ=255)
                                              aatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtttt                                                                                                                                                                                      <  2:165077‑M2/157‑1 (MQ=255)
                                                                        cccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAcc                                                                                                                                                            <  1:230943‑M2/183‑1 (MQ=255)
                                                                                     ccgggttattgcttgtcacaaacaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAACTCAACTTCCTGATTTTTAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCCTTCGACGTGCTGAAACCAAAACCTTTTTCTGGGCTAAACCATTTCCCTATACAAGTCATTTTgttt                                                                                                                                               >  2:6121‑M2/56‑250 (MQ=255)
                                                                                          ttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTAGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                                                                                                                                          >  2:133769‑M2/51‑251 (MQ=255)
                                                                                                                       tcatgcagttaacccactcacTACATGAACAGCGGCCGGACCTTTAGGTCCGTTCTCCATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                                                                                                                                                                                              >  2:260812‑M2/22‑138 (MQ=255)
                                                                                                                       tcatgcagttaacccacgcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTGTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                                                                                                                                                                                              <  1:260812‑M2/117‑1 (MQ=255)
                                                                                                                        catgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCtt                                                                                                            <  2:293542‑M2/231‑1 (MQ=255)
                                                                                                                              gttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGc                                                                                                        >  2:12553‑M2/15‑249 (MQ=255)
                                                                                                                                 aacccacttacgaCATGAACAGCGGCAGGACCTTTATGTCCGTTCTCAATACCAAATGTAACTTCCTGATTCGCAGCTAATGTTTTGAAATCGTTGCTCTGAATTGATGAGAAATAGACAAACACAGCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAATCCTTTTTCATGGTTAAACCATTTCACTAAACCATTCATTTTGTTAGACAGAATTATTACCTTTTGAAGAAATTAGCATTTGGGCAGAAt                                                                                                   <  2:217746/239‑1 (MQ=255)
                                                                                                                                                 gAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAata                                                                                   <  2:13963/251‑1 (MQ=255)
                                                                                                                                                                 tttAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAAcc                                                                    >  1:100563/1‑250 (MQ=255)
                                                                                                                                                                    aGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGACTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGACACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACcaac                                                                 >  1:330552/1‑250 (MQ=255)
                                                                                                                                                                                     cAAATTCAACTTGTTGATTCCCAGTTACTGTTTTGACTTCCCTGCTCTGCATTGCTGAGAAATCGACAAACACATCTTTGCTGCAATATTTCGGCGTGATGAAACCACAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATACTTCTTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGt                                               >  2:259306/1‑251 (MQ=255)
                                                                                                                                                                                          tCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACaa                                           >  1:404386/1‑250 (MQ=255)
                                                                                                                                                                                            aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATCGCTGACAAATAGACCAACCCATCTTTGCTGCCATCTTTCGGCGTGATGACACCACAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat                                         >  1:39345/1‑250 (MQ=255)
                                                                                                                                                                                                tCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataat                                      <  2:330552/249‑1 (MQ=255)
                                                                                                                                                                                                       tctcAGTTAATGTTTTGAAATCGTTGCTCTGAAGTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTTCCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAGCCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtat                              <  1:259514/250‑1 (MQ=255)
                                                                                                                                                                                                             ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATg                       >  1:430972/1‑251 (MQ=255)
                                                                                                                                                                                                                          aTCGTTGCTCTGAATTGCTGAGAAATATACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATCTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTTGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATa          <  1:529563/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGCGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATaa         >  2:323522/1‑250 (MQ=255)
                                                                                                                                                                                                                                     gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAATTCACTAAACCAGTCATTTTGTTAGACATAATTATTGCCTTTTTAAGAAATTAGCCCTGGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAGCAATCACGGCTTTAGATAAAGTTGTa  <  2:485880/248‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
TTGATTTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTA  >  NC_000913/1638327‑1638803

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: