Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,114,502	C→T	*465* (TAG→TAA)	wcaM ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,114,502	0	C	T	100.0%	248.2 / NA	76	*465* (TAG→TAA)	wcaM	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/38);  total (38/38)

TGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCA  >  NC_000913/2114367‑2114619                                                                                                                                              |                                                                                                                      tGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTAc                                                                                                                    >  2:1485897/1‑139 (MQ=255) tGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTAc                                                                                                                    >  2:1291929/1‑139 (MQ=255)   ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTa                                                                                                                     >  1:380925/1‑136 (MQ=255)   ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                  >  1:1127139/1‑139 (MQ=255)   ttCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTc                                                                                                                  >  1:401519/1‑139 (MQ=255)               cGCCTGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAAATTACTCCCTCAGTTCCgg                                                                                                      <  1:1114535/139‑1 (MQ=255)                          aTTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGt                                                                                           <  1:1291929/139‑1 (MQ=255)                                   aaTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                               >  2:1322446/1‑110 (MQ=255)                                   aaTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCt                                                                                                               <  1:1322446/110‑1 (MQ=255)                                   aaTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                         >  2:107721/1‑132 (MQ=255)                                   aaTTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                         <  1:107721/132‑1 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAtcttc                                                                                   >  1:1666553/1‑139 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAt                                                                                       >  1:58479/1‑129 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAt                                                                                       <  2:58479/129‑1 (MQ=255)                                        aTTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCcaca                                                                             >  2:1572587/1‑139 (MQ=255)                                         ttACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACAt                                                                            >  2:159637/1‑139 (MQ=255)                                            cTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATTGCAGTAGCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCca                                                                               <  2:1666553/139‑1 (MQ=255)                                                  tACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAatat                                                                   >  2:707344/1‑139 (MQ=255)                                                       gtgATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGa                                                              <  1:1533029/139‑1 (MQ=255)                                                         gATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTcc                                                                                                                 >  2:1387603/1‑86 (MQ=255)                                                         gATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTcc                                                                                                                 <  1:1387603/86‑1 (MQ=255)                                                         gATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                         <  2:1339787/126‑1 (MQ=255)                                                         gATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTaa                                                                         >  1:1339787/1‑126 (MQ=255)                                                              tACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                       >  2:833079/1‑91 (MQ=255)                                                              tACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCg                                                                                                       <  1:833079/91‑1 (MQ=255)                                                               aCGCTGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAt                                                                                       <  1:97719/106‑1 (MQ=255)                                                               aCGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAAt                                                                                       >  2:97719/1‑106 (MQ=255)                                                                cgccgcTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                  >  2:1296741/1‑126 (MQ=255)                                                                cgccgcTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                  <  1:1296741/126‑1 (MQ=255)                                                                      ttCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                      >  2:304052/1‑84 (MQ=255)                                                                      ttCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                      <  1:304052/84‑1 (MQ=255)                                                                           tgtgGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGAc                                          <  2:1690633/139‑1 (MQ=255)                                                                            gtgGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                         <  2:1127139/139‑1 (MQ=255)                                                                                         gCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCttt                                  >  2:1357894/1‑139 (MQ=255)                                                                                            gTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTc                               <  1:1426300/139‑1 (MQ=255)                                                                                            gTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTc                               <  1:1485897/139‑1 (MQ=255)                                                                                             tAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGTAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTcc                              >  2:1160266/1‑139 (MQ=255)                                                                                              aGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCt                             <  1:159637/139‑1 (MQ=255)                                                                                               gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                         >  1:813301/1‑78 (MQ=255)                                                                                               gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTa                                                                                         <  2:813301/78‑1 (MQ=255)                                                                                               gCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCtt                            >  2:1612876/1‑139 (MQ=255)                                                                                                     ttCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                      <  1:792387/59‑1 (MQ=255)                                                                                                     ttCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCgg                                                                                                      >  2:792387/1‑59 (MQ=255)                                                                                                          aaGTATGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCAcc                 >  1:54217/1‑139 (MQ=255)                                                                                                              aTGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCAt             <  2:401519/139‑1 (MQ=255)                                                                                                               tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCca                                                                               >  2:1552269/1‑72 (MQ=255)                                                                                                               tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCca                                                                               <  1:1552269/72‑1 (MQ=255)                                                                                                               tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAAc                                                                        >  2:925116/1‑79 (MQ=255)                                                                                                               tGACTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAAc                                                                        <  1:925116/79‑1 (MQ=255)                                                                                                                 aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                 >  2:34394/1‑52 (MQ=255)                                                                                                                 aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATc                                                                                                 <  1:34394/52‑1 (MQ=255)                                                                                                                 aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCt                                                                                                >  1:1213366/1‑53 (MQ=255)                                                                                                                 aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCt                                                                                                <  2:1213366/53‑1 (MQ=255)                                                                                                                 aCTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGaa          >  1:238888/1‑139 (MQ=255)                                                                                                                  cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGTAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAc         <  1:1160266/139‑1 (MQ=255)                                                                                                                  cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAc         <  1:1612876/139‑1 (MQ=255)                                                                                                                  cTCCATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAAc         <  2:54217/139‑1 (MQ=255)                                                                                                                     cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCAt             <  2:783889/132‑1 (MQ=255)                                                                                                                     cATTTTTCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCAt             >  1:783889/1‑132 (MQ=255)                                                                                                                       tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcatt                     >  2:163302/1‑122 (MQ=255)                                                                                                                       tttttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTttcatt                     <  1:163302/122‑1 (MQ=255)                                                                                                                         tttCCAGGAATGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCa  <  2:238888/139‑1 (MQ=255)                                                                                                                                  aTGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                  <  2:20270/66‑1 (MQ=255)                                                                                                                                  aTGGTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATg                                                                  >  1:20270/1‑66 (MQ=255)                                                                                                                                    ggTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAt                                                         <  1:61881/73‑1 (MQ=255)                                                                                                                                    ggTCGCAAATTTACTCCCTCAGTTCCGGCAATATAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAAt                                                         >  2:61881/1‑73 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACa                                                                       <  2:1414358/58‑1 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACa                                                                       >  1:1414358/1‑58 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAt                                                             <  2:1300190/68‑1 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGAt                                                             >  1:1300190/1‑68 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                                <  2:1337500/81‑1 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGat                                                >  1:1337500/1‑81 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                         >  1:968576/1‑88 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACg                                         <  2:968576/88‑1 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCa  <  2:1604969/127‑1 (MQ=255)                                                                                                                                     gTCGCAAATTTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCa  >  1:1604969/1‑127 (MQ=255)                                                                                                                                              |                                                                                                                      TGTTCGATTTTTTTCGCCAGCATCAATTACCCTGAATTGATTACTGAATTACTTGTGATGTTACGCCGCTTCGTTGTGGATT‑‑‑‑‑‑GCAGTAGCATTGTTCCTAAGTATGACTCCATTTTTCCAGGAATGGTCGCAAATCTACTCCCTCAGTTCCGGCAATCTAAAGTTAATCTTCTCCACATTAACAATATGGTGATTAATCCTGTCGATATCGACGGAGCTTTGTCCTTTTTCATTCACCGCATGAACATTTGCA  >  NC_000913/2114367‑2114619

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