Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	2,236,590	G→T	D127Y (GAT→TAT)	sanA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,236,590	0	G	T	100.0%	180.2 / NA	60	D127Y (GAT→TAT)	sanA	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (30/30);  total (30/30)

GATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGGATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATCACCCAACGTTTC  >  W3110S.gb/2236524‑2236655                                                                   |                                                                  aaTTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTc                                                               <  1:2743324/70‑1 (MQ=255)   tttAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCa                                                             <  1:2908314/71‑1 (MQ=255)         cGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGc                                                       >  1:1149054/1‑71 (MQ=255)          gctgctGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCg                                                      <  1:3953840/71‑1 (MQ=255)          gctgctGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCg                                                      <  1:4017605/71‑1 (MQ=255)                  tgtCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCg                                                      <  1:1841569/63‑1 (MQ=255)                   gtCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCa                                              >  1:2145237/1‑70 (MQ=255)                       atccATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGtt                                         <  1:2135319/69‑1 (MQ=255)                       atccATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGtt                                         <  1:2263835/69‑1 (MQ=255)                         ccATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGtttt                                       >  1:1550003/1‑71 (MQ=255)                          cATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTc                                      <  1:1851090/71‑1 (MQ=255)                          cATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTc                                      <  1:4074307/71‑1 (MQ=255)                          cATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTc                                      <  1:2288185/71‑1 (MQ=255)                            tCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGtttttc                                     >  1:3244838/1‑68 (MQ=255)                            tCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGa                                    >  1:1534038/1‑71 (MQ=255)                            tCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGa                                    <  1:2353023/71‑1 (MQ=255)                            tCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGa                                    >  1:2281702/1‑71 (MQ=255)                            ccAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTa                                                    <  1:2677296/54‑1 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTaca                                                  >  1:1775763/1‑55 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTaca                                                  >  1:4082515/1‑55 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCggt                                   >  1:4390163/1‑68 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGAt                                   >  1:3330434/1‑70 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGAt                                   >  1:616429/1‑70 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATa                                  >  1:3843150/1‑71 (MQ=255)                              aGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATa                                  >  1:4488552/1‑71 (MQ=255)                               gATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGa                                    <  1:1744974/68‑1 (MQ=255)                               gATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGa                                    <  1:3704362/68‑1 (MQ=255)                                atatTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAg                                           <  1:4584247/60‑1 (MQ=255)                                atatTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAg                                           <  1:1399528/60‑1 (MQ=255)                                atatTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAg                                           <  1:3269927/60‑1 (MQ=255)                                atatTGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGa                                    <  1:1745513/67‑1 (MQ=255)                                   ttGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTGCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             <  1:826569/71‑1 (MQ=255)                                   ttGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             <  1:535994/71‑1 (MQ=255)                                   ttGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             <  1:88468/71‑1 (MQ=255)                                   ttGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             <  1:441548/71‑1 (MQ=255)                                   ttGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             >  1:157360/1‑71 (MQ=255)                                    tGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             >  1:1545818/1‑70 (MQ=255)                                    tGTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             >  1:3552205/1‑70 (MQ=255)                                     gTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGa                           >  1:1333084/1‑71 (MQ=255)                                     gTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGa                           >  1:3959533/1‑71 (MQ=255)                                     gTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGa                           >  1:4601956/1‑71 (MQ=255)                                     gTTCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGa                           >  1:1722665/1‑71 (MQ=255)                                      ttCTCGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             <  1:3313194/68‑1 (MQ=255)                                        ctcGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGAttt                        >  1:320958/1‑71 (MQ=255)                                        ctcGATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGAttt                        >  1:4601728/1‑71 (MQ=255)                                           gATTACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCat                     <  1:791628/71‑1 (MQ=255)                                              tACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTa                                                    >  1:2702640/1‑37 (MQ=255)                                              tACGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCattat                  >  1:1846676/1‑71 (MQ=255)                                               aCGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCattat                  <  1:1390487/70‑1 (MQ=255)                                               aCGCAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCattat                  >  1:673902/1‑70 (MQ=255)                                               aCGCAGGCTTTCGTACGCTGTATTCAATCGTGCGTacac                                                 >  1:2328480/1‑39 (MQ=38)                                                  cAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTc                                      <  1:887186/47‑1 (MQ=255)                                                  cAGGCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATc              <  1:2311412/71‑1 (MQ=255)                                                    ggCTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATCAc            >  1:3759745/1‑71 (MQ=255)                                                      cTTTCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTcta                                    <  1:1809519/45‑3 (MQ=255)                                                        ttCGTACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATCACCCaa        >  1:3026173/1‑71 (MQ=255)                                                            tACGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAAt                             <  1:4658131/46‑1 (MQ=255)                                                              cGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATCACCCAACGTTTc  >  1:90211/1‑71 (MQ=255)                                                              cGCTGTATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATCACCCAACGTTTc  >  1:1324630/1‑71 (MQ=255)                                                              cGCTGTATTCCATCGTGCGTAAACGCAAAGTTTTCGATACTAATGATTTCATTATTATCACCCAACGTTTc  <  1:3792271/71‑1 (MQ=255)                                                                   |                                                                  GATTTAATCGCTGCTGGTGTCGACCCATCAGATATTGTTCTCGATTACGCAGGCTTTCGTACGCTGGATTCCATCGTGCGTACACGCAAAGTTTTCGATACTAATGATTTCATTATTATCACCCAACGTTTC  >  W3110S.gb/2236524‑2236655

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