Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	1,427,990	G→A	intergenic (+296/‑178)	ydaY → / → ynaA	hypothetical protein/predicted tail protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,427,990	0	G	A	100.0%	184.3 / NA	62	intergenic (+296/‑178)	ydaY/ynaA	hypothetical protein/predicted tail protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/33);  total (29/33)

CGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCG  >  W3110S.gb/1427923‑1428057                                                                    |                                                                    cGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAg                                                                  >  1:462759/1‑71 (MQ=255)     attgattgAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtg                                                              >  1:1995765/1‑71 (MQ=255)     attgattgAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtctg                                                              >  1:2728444/1‑71 (MQ=255)      ttgattgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAg                                                                  <  1:2722230/66‑1 (MQ=255)        gattgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtgtg                                                            <  1:3011401/70‑1 (MQ=255)         attgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAgtgtgtgt                                                           <  1:4281722/70‑1 (MQ=255)         attgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTa                                                          <  1:1307171/71‑1 (MQ=255)          ttgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAg                                                         >  1:174123/1‑71 (MQ=255)           tgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgt                                                        >  1:84188/1‑71 (MQ=255)           tgaGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgt                                                        >  1:180969/1‑71 (MQ=255)             aGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgtt                                                       <  1:4166201/70‑1 (MQ=255)               cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgttgtt                                                    <  1:997256/71‑1 (MQ=255)               cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgttgtt                                                    <  1:2254810/71‑1 (MQ=255)               cgcgATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgttgtt                                                    <  1:354421/71‑1 (MQ=255)                  gATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAAt                                                 <  1:1720457/71‑1 (MQ=255)                  gATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGATAAt                                                 <  1:3112108/71‑1 (MQ=255)                      aaTTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACg                                             <  1:3184771/71‑1 (MQ=255)                      aaTTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACg                                             <  1:2311004/71‑1 (MQ=255)                      aaTTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACg                                             <  1:2981091/71‑1 (MQ=255)                        tttGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTa                                           <  1:1357393/71‑1 (MQ=255)                        tttGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTa                                           <  1:1682147/71‑1 (MQ=255)                          tGAAGCAGTGATCTTTGCACTGGATGACCAGTATCTTCAGGAACAgtgtgtgt                                                           >  1:1196655/1‑53 (MQ=255)                                 gTGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                   >  1:4584415/1‑70 (MQ=255)                                  tGATCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                   <  1:3853755/69‑1 (MQ=255)                                    aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:1199288/1‑69 (MQ=255)                                    aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:2114432/1‑69 (MQ=255)                                    aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:4359962/1‑69 (MQ=255)                                    aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:3353017/1‑69 (MQ=255)                                    aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:1429922/1‑69 (MQ=255)                                    aTCTTTGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:2497778/1‑69 (MQ=255)                                    aTCTTTGAACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                 >  1:2230744/1‑69 (MQ=255)                                         tGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACagag                            <  1:285084/69‑1 (MQ=255)                                         tGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                             <  1:4543497/68‑1 (MQ=255)                                         tGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                             <  1:694276/68‑1 (MQ=255)                                         tGCACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                             <  1:2385942/68‑1 (MQ=255)                                          gCACTGGATGACCAGTATTTTCTGGAACAGTGTGTTTAGTTGTTAATTACGTACACTCTGt                                   <  1:1244430/61‑1 (MQ=255)                                           cACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgt                        <  1:4193536/71‑1 (MQ=255)                                           cACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgt                        <  1:1328909/71‑1 (MQ=255)                                           cACTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgt                        <  1:3571406/71‑1 (MQ=255)                                            aCTGGATGACCAGTATTTTCAGGAACAGTGTGTGTAgttgtt                                                    <  1:1859472/42‑1 (MQ=255)                                               ggATGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATg                    >  1:130164/1‑71 (MQ=255)                                                 aTGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGgt                  <  1:1860029/71‑1 (MQ=255)                                                 aTGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGgt                  >  1:2013786/1‑71 (MQ=255)                                                 aTGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGgt                  >  1:499605/1‑71 (MQ=255)                                                  tGACCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGgtg                 >  1:3468929/1‑71 (MQ=255)                                                    aCCAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGgtgt                <  1:4321061/70‑1 (MQ=255)                                                     ccAGTATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCt              >  1:3324574/1‑71 (MQ=255)                                                         tATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTaca                                         <  1:846393/40‑1 (MQ=255)                                                         tATTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                             <  1:1368799/52‑1 (MQ=255)                                                          aTTTTCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaat         >  1:2442958/1‑71 (MQ=255)                                                           ttttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaat         >  1:416898/1‑70 (MQ=255)                                                           ttttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaat         >  1:4343462/1‑70 (MQ=255)                                                           ttttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaat         >  1:3618690/1‑70 (MQ=255)                                                           ttttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaat         >  1:257214/1‑70 (MQ=255)                                                             ttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                   >  1:851650/1‑42 (MQ=255)                                                             ttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGa                      >  1:1398350/1‑55 (MQ=255)                                                             ttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaatta       >  1:969300/1‑70 (MQ=255)                                                             ttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaatta       >  1:2790494/1‑70 (MQ=255)                                                             ttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaatta       <  1:2998623/70‑1 (MQ=255)                                                             ttCAGGAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTttaatta       <  1:4425763/70‑1 (MQ=255)                                                                 ggAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCg  <  1:1317261/71‑1 (MQ=255)                                                                 ggAACAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCg  <  1:3501590/71‑1 (MQ=255)                                                                    |                                                                    CGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCG  >  W3110S.gb/1427923‑1428057

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