Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,897,439	G→A	*222* (TAG→TAA)	yieH →	phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,897,439	0	G	A	100.0%	252.3 / NA	78	*222* (TAG→TAA)	yieH	phosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/37);  total (41/37)

AGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGT  >  NC_000913/3897307‑3897573                                                                                                                                       |                                                                                                                                       aGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAAttcttc                                                                                                                                  >  2:1222847/1‑139 (MQ=255)  ggTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTca                                                                                                                                 <  2:1574639/139‑1 (MQ=255)   gTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                                <  2:410245/139‑1 (MQ=255)         aaGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTccc                                                                                                                          <  1:137201/139‑1 (MQ=255)             gtTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTTCCTTCa                                                                                                                      >  2:226379/1‑139 (MQ=255)             gtTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCa                                                                                                                      >  1:272851/1‑139 (MQ=255)              ttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcaca                                                                                                                               >  1:983816/1‑129 (MQ=255)              ttCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATGCTTcaca                                                                                                                               <  2:983816/129‑1 (MQ=255)                 tACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTTCCTTCACTTa                                                                                                                  <  1:226379/139‑1 (MQ=255)                      cTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCg                                                                                                             >  2:1107139/1‑139 (MQ=255)                        gcgcCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGTTATACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCg                                                                                                             >  1:965317/1‑139 (MQ=255)                             gACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAAt                                                                                                      <  1:250087/139‑1 (MQ=255)                              aCCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAACtt                                                                                                     >  2:521333/1‑139 (MQ=255)                                 cGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGgc                                                                                                  >  1:751500/1‑139 (MQ=255)                                   cacaATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGgcgc                                                                                                >  2:790961/1‑139 (MQ=255)                                      aaTAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             <  1:790961/139‑1 (MQ=255)                                          aGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         <  1:916179/139‑1 (MQ=255)                                          aGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         <  2:751500/139‑1 (MQ=255)                                          aGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         <  2:525331/139‑1 (MQ=255)                                                tCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgtaagta                                                                                   >  1:4916/1‑139 (MQ=255)                                                  gTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 <  1:1222847/139‑1 (MQ=255)                                                    tCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                                >  2:711187/1‑90 (MQ=255)                                                    tCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTcac                                                                                                                                <  1:711187/90‑1 (MQ=255)                                                         cGAAAGTCACCACCTTTACCCATCTTTCGCAGTTATACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAg                                                                            <  2:965317/139‑1 (MQ=255)                                                           aaaGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGCTTg                                                                        >  1:215089/1‑139 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAc                                                                                                                 <  1:518967/77‑1 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTAc                                                                                                                 >  2:518967/1‑77 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               >  1:315714/1‑95 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               <  2:315714/95‑1 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             <  2:1211707/97‑1 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             >  1:1211707/1‑97 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         <  2:1383173/101‑1 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         >  1:69412/1‑101 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         >  1:1383173/1‑101 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         <  2:69412/101‑1 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 <  2:1274457/109‑1 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTaaa                                                                                 >  1:1274457/1‑109 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                   >  2:1133222/1‑139 (MQ=255)                                                                                cTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCggg                                                   >  2:1363452/1‑139 (MQ=255)                                                                                 tttCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGt                                                  >  1:1188208/1‑139 (MQ=255)                                                                                 tttCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGt                                                  >  1:1233743/1‑139 (MQ=255)                                                                                 tttCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGt                                                  <  2:1188208/139‑1 (MQ=255)                                                                                    cGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTa                                                 <  1:1276864/137‑1 (MQ=255)                                                                                    cGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTa                                                 >  2:1276864/1‑137 (MQ=255)                                                                                         t‑‑tACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATc                                          <  2:272851/139‑1 (MQ=255)                                                                                             aCCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        >  2:299223/1‑123 (MQ=255)                                                                                             aCCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        <  1:299223/123‑1 (MQ=255)                                                                                                     tgtgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               >  2:1466811/1‑76 (MQ=255)                                                                                                     tgtgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               <  1:1466811/76‑1 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               >  2:1102466/1‑74 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               <  1:1102466/74‑1 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCa                                                          <  1:118407/111‑1 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCa                                                          >  2:118407/1‑111 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  2:1244430/1‑139 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  1:282413/1‑139 (MQ=255)                                                                                                       tgGAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAtt                              >  2:490764/1‑139 (MQ=255)                                                                                                         gAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                       <  1:1101158/80‑1 (MQ=255)                                                                                                         gAAAGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGgta                                                                                       >  2:1101158/1‑80 (MQ=255)                                                                                                           aaGCGCGTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGCTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGc                          <  2:215089/139‑1 (MQ=255)                                                                                                              cgcgTGGTTGGGATATTACGGCATAATTCTTCACAATCCCTTCACTTACCCCGCTTAAATTGACGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGt                       >  1:1558054/1‑139 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             >  1:1000106/1‑66 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCa                                                                                             <  2:1000106/66‑1 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        <  1:1284686/103‑1 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACg                                                        >  2:1284686/1‑103 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATc                                          <  2:947181/117‑1 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATc                                          >  1:947181/1‑117 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtttt                    >  2:949525/1‑139 (MQ=255)                                                                                                                 gTGGTTGGGATATTACGACATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTGAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGtttt                    >  1:921557/1‑139 (MQ=255)                                                                                                                   ggTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAg                  <  2:282413/139‑1 (MQ=255)                                                                                                                   ggTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAg                  <  1:1363452/139‑1 (MQ=255)                                                                                                                    gTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         >  2:1323056/1‑67 (MQ=255)                                                                                                                    gTTGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAgg                                                                                         <  1:1323056/67‑1 (MQ=255)                                                                                                                     ttGGGATATTACGGCATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa                <  2:291718/139‑1 (MQ=255)                                                                                                                                   cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               <  1:1321458/46‑1 (MQ=255)                                                                                                                                   cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCt                                                                                               >  2:1321458/1‑46 (MQ=255)                                                                                                                                   cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAATGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGa                <  2:829939/125‑1 (MQ=255)                                                                                                                                   cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGgtgt  >  2:265578/1‑139 (MQ=255)                                                                                                                                   cATAATTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGGTTTAGGa                >  1:829939/1‑125 (MQ=255)                                                                                                                                       |                                                                                                                                       AGGTATGGAAGTGTTCTACTTCTGCGCCGACCCGCACAATAAGCCGATCGTTCACCCGAAAGTCACCACCTTTACCCATCTTTCGCAGT‑‑TACCTGAACTGTGGAAAGCGCGTGGTTGGGATATTACGGCATAGTTCTTCACACTCCCTTCACTTACCCCGCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGT  >  NC_000913/3897307‑3897573

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