Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,810,941	T→C	G76G (GGA→GGG)	b1730 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,810,941	0	T	C	100.0%	116.0 / NA	38	G76G (GGA→GGG)	b1730	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/10);  total (28/10)

TATCCGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCTCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTCACATTG  >  NC_000913/1810851‑1811030                                                                                           |                                                                                          tATCCGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGcccc                                                                                         >  1:3918138/1‑93 (MQ=255)  aTCCGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCt                                                                                        >  2:866077/1‑93 (MQ=255)   tCCGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTg                                                                                       >  2:384947/1‑93 (MQ=255)     cGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGc                                                                                     >  1:2929994/1‑93 (MQ=255)       aaTTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCtt                                                                                   <  2:1177316/93‑1 (MQ=255)        aTTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCttt                                                                                  <  1:4257400/93‑1 (MQ=255)          ttCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGa                                                                                <  2:4536402/93‑1 (MQ=255)                 ccATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGcatca                                                                         >  2:490150/1‑93 (MQ=255)                        gATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTaa                                                                  >  2:2530765/1‑93 (MQ=255)                           tattaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGt                                                               >  2:2384658/1‑93 (MQ=255)                           tattaattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGt                                                               <  2:1024900/93‑1 (MQ=255)                               aattattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGa                                                           <  1:3553579/93‑1 (MQ=255)                                attattaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCCTTGAGGCATCATAGTTAATGTTAGAc                                                          >  1:335168/1‑93 (MQ=255)                                   attaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACata                                                       <  2:1563657/93‑1 (MQ=255)                                   attaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACata                                                       >  2:3122363/1‑93 (MQ=255)                                   attaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACata                                                       <  1:3867357/93‑1 (MQ=255)                                    ttaCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACatat                                                      >  2:1968063/1‑93 (MQ=255)                                     taCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACAtata                                                     >  2:3893642/1‑93 (MQ=255)                                     taCTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACAtata                                                     >  1:2468660/1‑93 (MQ=255)                                         ctATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCa                                                 >  1:1439003/1‑93 (MQ=255)                                          tataTTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCACCTGGCTTTGAGGCATCATAGTTAATGTTAGACATCTAGGCaa                                                >  2:2335938/1‑93 (MQ=255)                                            tatTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATc                                              >  2:1998868/1‑93 (MQ=255)                                                tAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCatat                                          >  1:2177290/1‑93 (MQ=255)                                                tAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCatat                                          >  1:2415223/1‑93 (MQ=255)                                                    aTTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAAtt                                      >  2:181441/1‑93 (MQ=255)                                                    aTTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAAtt                                      >  2:4189686/1‑93 (MQ=255)                                                       gTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTaaa                                   <  1:230814/93‑1 (MQ=255)                                                        tCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAg                                  >  1:4464446/1‑93 (MQ=255)                                                        tCCTGCAGATATGCATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAg                                  <  1:2678947/93‑1 (MQ=255)                                                           tGCAGATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGttt                               >  1:1744911/1‑93 (MQ=255)                                                               gATATGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGAtt                           >  2:702557/1‑93 (MQ=255)                                                                 tatGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCt                         <  2:3033122/93‑1 (MQ=255)                                                                  atGTATCAAAACTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTg                        >  1:3133662/1‑93 (MQ=255)                                                                           aaCTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCAt               >  2:2365118/1‑93 (MQ=255)                                                                           aaCTAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCAt               >  2:475170/1‑93 (MQ=255)                                                                              tAAGTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGc            >  1:1085213/1‑93 (MQ=255)                                                                                 gTGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGt         >  2:4323517/1‑93 (MQ=255)                                                                                  tGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTGGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTc        >  1:3385363/1‑93 (MQ=255)                                                                                  tGCTGGCGCCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTc        >  1:4376035/1‑93 (MQ=255)                                                                                        cgcCCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTCACATTg  >  2:4319851/1‑93 (MQ=255)                                                                                           |                                                                                          TATCCGAATTTCACTACCATAAAGATTATTAATTATTACTCTATATTTAACATTGTCCTGCAGATATGTATCAAAACTAAGTGCTGGCGCTCCTGGCTTTGAGGCATCATAGTTAATGTTAGACATATAGGCAATCATATAATTAAAGTTTGATTCTGCTATTTCATCGCAGTCACATTG  >  NC_000913/1810851‑1811030

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