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%	114.2 / NA	33	H33Y (CAT→TAT)	cspI	cold shock‑like protein CspI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/21);  total (12/21)

GATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGA  >  NC_000913/1638340‑1638813                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    atgatgccggttacccaaagctatgccggactgaatgtccaccgtcagtaatttttataccctgcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATTCCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTtgctg                                                                                                                                                                                                                                 <  2:445745‑M2/124‑1 (MQ=255)       ccggttacccaaagctatgccggactgaatgtccaccggcagtaatttttataccccgcgtaactgccgggtgattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTGTGAAATCGTTGCTCTTAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCt                                                                                                                                                                                                                           <  1:213512‑M2/130‑1 (MQ=255)                 aaagctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGGTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAC‑AACACCATCTTTGCTGCCATCTTTCGGCgtgg                                                                                                                                                                                                                 >  2:257016‑M2/112‑250 (MQ=255)                    gctatgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTACTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGtgatg                                                                                                                                                                                                               >  2:87377‑M2/109‑250 (MQ=255)                       atgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGaa                                                                                                                                                                                                             <  1:428537‑M2/144‑1 (MQ=255)                       atgccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGGCCGTTCTCAATACCAAATTAAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAAAAACA‑CATCTTTTCTGCCATCTTTCGGCGTGATGAAAc                                                                                                                                                                                                           >  2:322888‑M2/106‑251 (MQ=255)                          ccggactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCCACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                                                                                                                                                                                                         >  1:323413‑M2/103‑250 (MQ=255)                             gactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAAccaaagc                                                                                                                                                                                                     <  1:388762‑M2/152‑3 (MQ=255)                             gactgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                                                                                                                                                     <  1:388773‑M2/152‑1 (MQ=255)                                tgaatgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                                                                                                                                                                  >  1:111617‑M2/97‑251 (MQ=255)                                   atgtccaccgtcagtaatttttatacccggcgtaactgccgggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                                                                                                                                                               <  2:111617‑M2/158‑1 (MQ=255)                                                                     actgcctggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGTACCTTTAGGTCCGTTCTCAATTCCAAATTCAACTTCCTGAATCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTGTTCAGGGTTAAACCATTTCACTAAACCAGTCAtttt                                                                                                                                                             <  2:431872‑M2/192‑1 (MQ=255)                                                                           gggttattgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa                                                                                                                                                       <  2:319049‑M2/198‑1 (MQ=255)                                                                            ggttattgcttgtcaaaaaaaagtggtagagncatgcagataaccaaatcacTACATGAACAGCGGCAGGACCTTTATGTCCGTTCTCAATGCTAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGCCAACCA‑CATCGTTGCTGTCATCTTTCGGCGTGATGAAACCAAATCCTTTTTCAGGGTGAAACCATTTCCCTAAACCAGTCATTTTGTTAGAc                                                                                                                                                      <  2:215208‑M2/199‑1 (MQ=255)                                                                                attgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAATCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATa                                                                                                                                                   <  1:53881‑M2/202‑1 (MQ=255)                                                                                attgcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATa                                                                                                                                                   >  2:435923‑M2/49‑250 (MQ=255)                                                                                   gcttgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAATCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAGTTTGTTAGACATAatta                                                                                                                                               <  2:36923‑M2/206‑1 (MQ=255)                                                                                      tgtcacaaaaaagtggtagactcatgcagttaactcactcacTACATGAACAGCGGCAGGACCTTTAGGTTCGTTCTCAATACCAAATTCAACTGCCTGATTCTCAGTGAATGTTTTGAAATCGTTGCTCGGAATTGCTGAGAAATAGACAATCA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAATCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattat                                                                                                                                              <  2:166739‑M2/207‑1 (MQ=255)                                                                                      tgtcacaaaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATTCCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCTGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAattatt                                                                                                                                             <  1:353485‑M2/208‑1 (MQ=255)                                                                                          acaaaaaagtggtagactcatgcagttaacccactcacTACAAGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATCCCAAATTCAACTTCCTTATTCTCAGTTAATGTTTTGAAATCGTTTCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTTCCATCTTTCGGCGTGATGAAACCAAATCCTGTTTCATGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt                                                                                                                                        <  2:340220‑M2/213‑1 (MQ=255)                                                                                             aaaaagtggtagactcatgcagttaacccactcacTACATGAACAGCGGCAGGACCTTTAGGCCCGTTCTCAATACAAAATTCAAATTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGGAATGAAACCAAAACATTTTTCCGGGTTAAACCATTTCACTAAACCAGTCATTTTTTTAGCCATAATTATTACCTTTTg                                                                                                                                     >  2:250215‑M2/36‑251 (MQ=255)                                                                                                                          actcacTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTcc                                                                                                        <  2:322350‑M2/245‑1 (MQ=255)                                                                                                                                                      ttAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACca                                                                             <  2:80896/250‑1 (MQ=255)                                                                                                                                                                        ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCANAGAGAAAAACCAACAAGGAACTCTCAAGAgg                                                          >  2:229693/1‑251 (MQ=255)                                                                                                                                                                          aaaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAgga                                                         >  2:347302/1‑249 (MQ=255)                                                                                                                                                                               cAACTTCCTTATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat                                                   <  1:359802/251‑1 (MQ=255)                                                                                                                                                                               cAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaat                                                   <  1:194053/251‑1 (MQ=255)                                                                                                                                                                               cAACTGCCAGATTCGCAGTTAATGTTTTGAAATCGTGGCTCTGAATTGCTGAGAAATAGACAAACA‑TATCTTTGCTGCCATCTTTCGGCGTGATGAAACGAAAACATTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAAATAATTATTACCTTTTGAAGAAATTAGCCCTTTGGTAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAAACAACAATGATTTCTCAAGATGTATAaat                                                   <  2:393391/251‑1 (MQ=255)                                                                                                                                                                                 aCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaataa                                                 <  1:435923/251‑1 (MQ=255)                                                                                                                                                                                 aCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAaata                                                  >  2:456173/1‑250 (MQ=255)                                                                                                                                                                                       tGATTATCAGCTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGACATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAACTCTCAAGAGGTACAAATAATAAAAt                                           >  2:360608/1‑251 (MQ=255)                                                                                                                                                                                                                       ctGAATTGCTGAGAAATAGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTa            >  1:21920/1‑250 (MQ=255)                                                                                                                                                                                                                                 gagaAATAGACAAACC‑CTTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAATCCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAACTATTACCATTTGAAAAAATTAGCCCTTAGGCAGAATGGTCCGAAAAAAAATAACAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGa  <  2:115809/250‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    GATACATTTTGCCCCGAAAGGATTCATGTCAGAAACATCAAAAGATGATGTTCTATACTGGGAACCATGACAACACGGGCATCTAAAGTGAATATGGTTTGTAATATTGTCTACCTCAAAGCGCCACTACATGAACAGCGGCAGGACCTTTAGGTCCGTTCTCAATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACA‑CATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTTGAAGAAATTAGCCCTTGGGCAGAATGGTCCGAAAAAAAATATCAGAGAGAAAAACCAACAAGGAAATCTCAAGAGGTACAAATAATAAAATTATAACAATGACTGCTTCAGATAAATTTGTAACAAACCAGA  >  NC_000913/1638340‑1638813

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