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_000913	1,638,571	0	G	A	100.0%	96.7 / NA	26	H33Y (CAT→TAT)	cspI	cold shock‑like protein CspI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/16);  total (10/16)

TGATTTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACC  >  NC_000913/1638328‑1638810                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                aatacactggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAac                                                                                                                                                                                                                                          >  1:214466‑M2/140‑251 (MQ=255)  ataccctggaaatgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAac                                                                                                                                                                                                                                          >  2:53521‑M2/139‑250 (MQ=255)                       ttacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCGGCCATCTTTCg                                                                                                                                                                                                                    >  2:335448‑M2/118‑251 (MQ=255)                        tacccaaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCgg                                                                                                                                                                                                                   >  2:234451‑M2/117‑251 (MQ=255)                              aagctatgccggactgaatgtccgccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATGCCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCGTTGCTGCCATCTTTCGGCGtgat                                                                                                                                                                                                             <  2:114694‑M2/141‑1 (MQ=255)                                         gactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCTGCGTGATGAAACCAAAAc                                                                                                                                                                                                  <  1:87422‑M2/152‑1 (MQ=255)                                               atgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                                                                                                                                            <  1:294012‑M2/158‑1 (MQ=255)                                                  gctaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtatactcatgcagttaacccactcacTACATGAACAGCGGCAGGAACTTTAGGTTCGTTCTCAATACCAAATTCAACTTCCTGATTCTAAGGTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACCTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAATCCTTTTTCAg                                                                                                                                                                                         <  2:230125‑M2/161‑1 (MQ=255)                                                       cgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa                                                                                                                                                                                    <  1:32169‑M2/166‑1 (MQ=255)                                                                tttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                                                                                                                                                           >  2:251392‑M2/77‑251 (MQ=255)                                                                                    gccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                                                                                                                                                         <  2:335582‑M2/193‑1 (MQ=255)                                                                                                     cacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAcc                                                                                                                                       <  1:210837‑M2/211‑1 (MQ=255)                                                                                                        aaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt                                                                                                                                   <  1:288768‑M2/215‑1 (MQ=255)                                                                                                          aaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATCGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTg                                                                                                                                  >  1:45268‑M2/35‑250 (MQ=255)                                                                                                           aaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTgaa                                                                                                                                <  1:53521‑M2/218‑1 (MQ=255)                                                                                                                     actcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGcc                                                                                                                      >  1:51575‑M2/24‑251 (MQ=255)                                                                                                                                 aacccacttacTACATGAACAGCGGCAGGACCTTTAGTACAGTTCTCAATGCCATATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACAGCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACTTAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCCGAAt                                                                                                          <  2:240311‑M2/240‑1 (MQ=255)                                                                                                                                    ccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTACCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGt                                                                                                       <  1:262594‑M2/243‑1 (MQ=255)                                                                                                                                         cACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTTTCAATACCAAATTCAACGTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCg                                                                                                    <  1:205894/249‑1 (MQ=255)                                                                                                                                                                             ctcAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAAtct                                                              >  2:268196/1‑251 (MQ=255)                                                                                                                                                                                                        cgcAGTTAATGGTTTGAAATCGTTGTGCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATtata                                    <  2:306799/248‑1 (MQ=255)                                                                                                                                                                                                           aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACa                                 >  1:299501/1‑250 (MQ=255)                                                                                                                                                                                                            gTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACCCATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTGTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAAt                               <  2:58934/251‑1 (MQ=255)                                                                                                                                                                                                                                gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAAttt            <  2:307108/250‑1 (MQ=255)                                                                                                                                                                                                                                    ggAATTGCTGAAAAATAGACAAAAACATCTTTGATGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTACAGGGTTAAACCATTTCACTAAACCAGGCATTTAGATAGACATAACTAGTGACTTTTGAAGAAATTAGCCCTTGGGAAGAATGGTCCGAAAAAAAATATCAGAGATAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAAc       <  2:184169/250‑1 (MQ=255)                                                                                                                                                                                                                                         tGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAAcc  >  1:364563/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                TGATTTACAAAAGATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACC  >  NC_000913/1638328‑1638810

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