Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,884,028	G→A	M1M (ATG→ATA) †	ysdD →	protein YsdD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,884,028	0	G	A	100.0%	142.5 / NA	41	M1M (ATG→ATA) †	ysdD	protein YsdD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/26);  total (15/26)

CCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATGACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCC  >  NC_000913/3883893‑3884161                                                                                                                                        |                                                                                                                                      ccTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCa                                                                                                                                   <  2:369207/140‑1 (MQ=255) ccTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCa                                                                                                                                   <  2:350062/140‑1 (MQ=255)    aaaCTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAg                                                                                                                                >  1:523172/1‑140 (MQ=255)    aaaCTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAg                                                                                                                                >  1:714723/1‑140 (MQ=255)         gcgcAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACaaa                                                                                                                           >  1:293365/1‑140 (MQ=255)          cgcAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACaaaa                                                                                                                          >  2:71529/1‑140 (MQ=255)          cgcAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACaaaa                                                                                                                          >  1:544642/1‑140 (MQ=255)          cgcAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACaaaa                                                                                                                          >  2:58181/1‑140 (MQ=255)               gagaTCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTgg                                                                                                                     >  2:350374/1‑140 (MQ=255)                     ttctGTTTTTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  1:296233/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  1:71529/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  2:523172/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  2:293365/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  2:556375/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  2:688926/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  2:714723/140‑1 (MQ=255)                     ttctGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  1:473872/140‑1 (MQ=255)                         gTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGAt                                                                                                           <  1:466158/140‑1 (MQ=255)                          tttCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATc                                                                                                          >  1:347048/1‑140 (MQ=255)                              tcACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCtaat                                                                                                      <  1:78478/140‑1 (MQ=255)                                    aTTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               >  2:306647/1‑125 (MQ=255)                                    aTTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAAt                                                                                                               <  1:306647/125‑1 (MQ=255)                                     ttttCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATc                                                                                               <  1:350374/140‑1 (MQ=255)                                                     acaGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTgcg                                                                               <  2:217153/140‑1 (MQ=255)                                                       aGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCt                                                                             >  1:264439/1‑140 (MQ=255)                                                       aGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGACGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAGTAAAGATCCAGGACGATCCTAGCGCt                                                                             >  2:710922/1‑140 (MQ=255)                                                        ggACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCtt                                                                            >  2:435310/1‑140 (MQ=255)                                                                cAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATc                                                                    <  1:710922/140‑1 (MQ=255)                                                                  gAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAg                                                                  <  2:254678/140‑1 (MQ=255)                                                                   aaCTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGc                                                                 <  2:708586/140‑1 (MQ=255)                                                                     cTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGccc                                                               <  2:154200/140‑1 (MQ=255)                                                                        gCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGta                                                            >  1:364291/1‑140 (MQ=255)                                                                                       tCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCgg                                             >  2:43539/1‑140 (MQ=255)                                                                                           gCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGgcgcg                                         <  2:264439/140‑1 (MQ=255)                                                                                                gAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATg                                    <  2:544642/140‑1 (MQ=255)                                                                                                 aGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGc                                   >  1:54104/1‑140 (MQ=255)                                                                                                     gggCACGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGt                               <  2:54104/140‑1 (MQ=255)                                                                                                          cGATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCa                          <  2:387437/140‑1 (MQ=255)                                                                                                           gATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCAc                         <  2:19278/140‑1 (MQ=255)                                                                                                           gATTTACGCCGCATATTGGAAAATTTAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCAc                         <  2:795747/140‑1 (MQ=255)                                                                                                                                  tttAATAACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTcc  <  2:149416/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      CCTAAACTGCGCAAGAGATCTTCTGTTTCTCACAGATTTTCCCGATTTATCCACAGGACTTTCCAGAACTCGCTAAGTGTAAACGATCCTGCCGCGAGGCGGGCACGATTTACGCCGCATATTGGAAAATTTAATGACCATAGACAAAAATTGGCTTAATCGATCTAATAAAGATCCAGGACGATCCTTGCGCTTTACCCATCAGCCCGTATAATCCTCCACCCGGCGCGCCATGCTGGTTTCCACTGGTGTGAGGTCGTACATTTTCC  >  NC_000913/3883893‑3884161

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