Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,743,076	A→G	N115S (AAC→AGC)	yiaK →	2,3‑diketo‑L‑gulonate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,743,076	0	A	G	100.0%	286.8 / NA	81	N115S (AAC→AGC)	yiaK	2,3‑diketo‑L‑gulonate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (46/35);  total (46/35)

TAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATC  >  NC_000913/3742941‑3743212                                                                                                                                        |                                                                                                                                         taACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCAc                                                                                                                                      >  2:1010434/1‑140 (MQ=255) taACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCAc                                                                                                                                      >  1:862476/1‑140 (MQ=255) taACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCAc                                                                                                                                      >  2:404137/1‑140 (MQ=255)  aaCCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACt                                                                                                                                     >  1:623861/1‑140 (MQ=255)      aGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGAt                                                                                                                                 >  2:845053/1‑140 (MQ=255)        ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATgc                                                                                                                               >  2:1247940/1‑140 (MQ=255)        ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATgc                                                                                                                               >  1:985855/1‑140 (MQ=255)        ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATgc                                                                                                                               >  2:579763/1‑140 (MQ=255)        ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATgc                                                                                                                               >  2:630959/1‑140 (MQ=255)        ccTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATgc                                                                                                                               >  2:692785/1‑140 (MQ=255)         cTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCG‑‑AAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGcgcg                                                                                                                            >  1:1035360/1‑140 (MQ=255)            ggCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCgcgg                                                                                                                           >  2:290817/1‑140 (MQ=255)            ggCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCgcgg                                                                                                                           >  2:834055/1‑140 (MQ=255)                   ttGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGc                                                                                                                    <  2:985855/140‑1 (MQ=255)                      aaCAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTAc                                                                                                                 >  2:116823/1‑140 (MQ=255)                           tGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTg                                                                                                            >  1:119375/1‑140 (MQ=255)                           tGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTg                                                                                                            >  1:894389/1‑140 (MQ=255)                              gACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCa                                                                                                         <  1:301852/140‑1 (MQ=255)                                 gCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAggc                                                                                                      >  2:413452/1‑140 (MQ=255)                                   ccAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAggcgg                                                                                                    >  1:768456/1‑140 (MQ=255)                                        gTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGaa                                                                                               >  1:1006095/1‑140 (MQ=255)                                        gTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGaa                                                                                               <  2:862476/140‑1 (MQ=255)                                         ttCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGaaa                                                                                              <  1:652373/140‑1 (MQ=255)                                           cGATCGGTAACCTGACAGCG‑‑AAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          <  2:1035360/140‑1 (MQ=255)                                            gATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAg                                                                                           <  2:320131/140‑1 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          <  1:116823/140‑1 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          >  1:991694/1‑140 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          <  1:210458/140‑1 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          <  1:290817/140‑1 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          <  1:845053/140‑1 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          >  1:1194646/1‑140 (MQ=255)                                             aTCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAgg                                                                                          <  1:404137/140‑1 (MQ=255)                                                 gTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAggcgg                                                                                                    <  1:646513/126‑1 (MQ=255)                                                 gTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAggcgg                                                                                                    >  2:646513/1‑126 (MQ=255)                                                  tAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCtata                                                                                     <  1:579763/140‑1 (MQ=255)                                                       tGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGc                                                                                <  2:848695/140‑1 (MQ=255)                                                          cAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCAtt                                                                             >  2:341658/1‑140 (MQ=255)                                                           aGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCtata                                                                                     <  1:426010/131‑1 (MQ=255)                                                           aGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCtata                                                                                     >  2:426010/1‑131 (MQ=255)                                                            gCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTg                                                                           >  1:78303/1‑140 (MQ=255)                                                             cGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc                                                                          <  2:1194646/140‑1 (MQ=255)                                                             cGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGc                                                                          <  2:991694/140‑1 (MQ=255)                                                               aaaaaGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGc                                                                                         <  2:870226/123‑1 (MQ=255)                                                               aaaaaGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGc                                                                                         >  1:870226/1‑123 (MQ=255)                                                                   aGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGAcc                                                                    <  1:891990/140‑1 (MQ=255)                                                                           gATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCAt                                                            <  1:1010434/140‑1 (MQ=255)                                                                              cgcgCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGc                                                         >  2:995953/1‑140 (MQ=255)                                                                               gcgcCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGcc                                                        >  1:516156/1‑140 (MQ=255)                                                                               gcgcCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGcc                                                        >  1:835726/1‑140 (MQ=255)                                                                                   cATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTaa                                                    <  1:834055/140‑1 (MQ=255)                                                                                             gCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATgccg                                               <  1:1158469/135‑1 (MQ=255)                                                                                             gCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATgccg                                               >  2:1158469/1‑135 (MQ=255)                                                                                               tGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTggg                                        >  2:516027/1‑140 (MQ=255)                                                                                                  cGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTggg                                        >  2:708657/1‑137 (MQ=255)                                                                                                  cGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTggg                                        <  1:708657/137‑1 (MQ=255)                                                                                                       aCGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCaaaa                                >  2:680598/1‑140 (MQ=255)                                                                                                       aCGGTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCaaaa                                <  1:260255/140‑1 (MQ=255)                                                                                                         ggTATTGGTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAaga                              <  2:894389/140‑1 (MQ=255)                                                                                                                gTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCaaaa                                >  1:1043128/1‑131 (MQ=255)                                                                                                                gTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCaaaa                                <  2:1043128/131‑1 (MQ=255)                                                                                                                gTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGc                       >  1:311919/1‑140 (MQ=255)                                                                                                                gTCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGc                       >  1:312689/1‑140 (MQ=255)                                                                                                                 tCTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCa                      <  2:311919/140‑1 (MQ=255)                                                                                                                  cTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCAt                     <  2:78303/140‑1 (MQ=255)                                                                                                                  cTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCAt                     <  2:768456/140‑1 (MQ=255)                                                                                                                  cTGGTGGCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCAt                     <  2:623861/140‑1 (MQ=255)                                                                                                                   tggtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATa                    >  2:118487/1‑140 (MQ=255)                                                                                                                     gtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAaga                              >  2:585646/1‑128 (MQ=255)                                                                                                                     gtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAaga                              <  1:585646/128‑1 (MQ=255)                                                                                                                     gtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAgg                  >  2:558215/1‑140 (MQ=255)                                                                                                                     gtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAgg                  >  1:812233/1‑140 (MQ=255)                                                                                                                     gtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAgg                  >  1:1234360/1‑140 (MQ=255)                                                                                                                     gtggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAgg                  >  2:1062290/1‑140 (MQ=255)                                                                                                                       ggCACTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCa                <  1:1247940/140‑1 (MQ=255)                                                                                                                          aCTACGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTa             >  2:106346/1‑140 (MQ=255)                                                                                                                              cGTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAAccc         >  1:128446/1‑140 (MQ=255)                                                                                                                               gTAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCg        >  1:93281/1‑140 (MQ=255)                                                                                                                                tAATGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGc       <  2:516156/140‑1 (MQ=255)                                                                                                                                   tGCCAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGa    >  2:33019/1‑140 (MQ=255)                                                                                                                                    gccagccaCTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGAt   <  1:413452/140‑1 (MQ=255)                                                                                                                                    gccagccaCTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGAt   <  2:1234360/140‑1 (MQ=255)                                                                                                                                     ccAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGAtc  <  2:812233/140‑1 (MQ=255)                                                                                                                                     ccAGCCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGAtc  <  1:995953/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         TAACCAGCCTCGGCGCAATTGAACAGTGGGACGCCCAGCGTTCGATCGGTAACCTGACAGCGAAAAAGATGATGGATCGCGCCATTGAACTGGCTGCCGATCACGGTATTGGTCTGGTGGCACTACGTAATGCCAACCACTGGATGCGCGGCGGCAGCTACGGCTGGCAGGCGGCGGAAAAAGGCTATATTGGCATTTGCTGGACCAACTCCATCGCCGTAATGCCGCCGTGGGGCGCAAAAGAGTGTCGCATAGGCACTAACCCGCTGATC  >  NC_000913/3742941‑3743212

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: ‑