Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,133,868	A→C	A595A (GCT→GCG)	wzc ←	colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,133,868	0	A	C	100.0%	275.9 / NA	70	A595A (GCT→GCG)	wzc	colanic acid production tyrosine‑protein kinase; autokinase; Ugd phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (37/33);  total (37/33)

ATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGG  >  NC_000913/2133727‑2134015                                                                                                                                              |                                                                                                                                                    atcaatcaACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGgtagta                                                                                                                                            <  2:87325/151‑1 (MQ=255)  tcaatcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGgtagta                                                                                                                                            >  1:579980/1‑150 (MQ=255)     atcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAatat                                                                                                                                        <  2:767675/151‑1 (MQ=255)     atcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAatat                                                                                                                                        <  1:768709/151‑1 (MQ=255)      tcaaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGCCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATc                                                                                                                                       >  1:765173/1‑151 (MQ=255)       caaCACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCg                                                                                                                                      >  2:74775/1‑151 (MQ=255)         acacCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGcc                                                                                                                                    >  2:1118267/1‑151 (MQ=255)          cacCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGccc                                                                                                                                   <  1:394169/151‑1 (MQ=255)             cAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGa                                                                                                                                <  1:786908/151‑1 (MQ=255)              aGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGAc                                                                                                                               >  1:452080/1‑151 (MQ=255)                gTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAACGCGCTCGATCATCAACAGTTAAGAAGGATATGGCGGCACCTGACCGCGCGGGATCAGGCCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGAAGTAATATCGCCCTGAcc                                                                                                                              >  2:1100267/1‑150 (MQ=255)                gTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCa                                                                                                                             <  2:1097944/151‑1 (MQ=255)                  cATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCaa                                                                                                                            <  2:23166/150‑1 (MQ=255)                   aTAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATc                                                                                                                          <  1:330143/151‑1 (MQ=255)                     aGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAAt                                                                                                                           >  2:989561/1‑148 (MQ=255)                     aGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCa                                                                                                                         <  1:535478/150‑1 (MQ=255)                             cTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTcc                                                                                                                 >  1:1003101/1‑150 (MQ=255)                               cGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCg                                                                                                                <  2:645036/149‑1 (MQ=255)                                gCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGAc                                                                                                              >  1:420603/1‑150 (MQ=255)                                gCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACa                                                                                                             <  2:631860/151‑1 (MQ=255)                                       tCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCAt                                                                                                       <  1:142630/150‑1 (MQ=255)                                             gTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGt                                                                                                 <  1:210421/150‑1 (MQ=255)                                              ttCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGtt                                                                                                >  2:960840/1‑150 (MQ=255)                                                 ggCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGttatt                                                                                             >  1:1080208/1‑150 (MQ=255)                                                   cAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGt                                                                                          >  1:87452/1‑151 (MQ=255)                                                    aaaGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTg                                                                                         >  1:616117/1‑151 (MQ=255)                                                     aaGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGc                                                                                        >  2:129309/1‑151 (MQ=255)                                                      aGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGcc                                                                                       >  2:322254/1‑151 (MQ=255)                                                            cGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAAc                                                                                   >  1:465374/1‑149 (MQ=255)                                                                    aaCAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTggg                                                                         <  2:179558/151‑1 (MQ=255)                                                                     aCAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGgt                                                                        >  2:675696/1‑151 (MQ=255)                                                                        gTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGgtgt                                                                      >  2:598991/1‑150 (MQ=255)                                                                              aaGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCttt                                                               >  2:1063181/1‑151 (MQ=255)                                                                               aGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTg                                                              <  1:446707/151‑1 (MQ=255)                                                                                   tttGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTgcgc                                                           >  2:829045/1‑150 (MQ=255)                                                                                   tttGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTgcgc                                                           >  2:516308/1‑150 (MQ=255)                                                                                     tGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCat                                                         >  2:237230/1‑150 (MQ=255)                                                                                       gcggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATc                                                      <  2:18800/151‑1 (MQ=255)                                                                                       gcggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATc                                                      >  2:162611/1‑151 (MQ=255)                                                                                         ggcACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGc                                                    >  2:148444/1‑151 (MQ=255)                                                                                            aCCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGt                                                 <  2:1065237/151‑1 (MQ=255)                                                                                             ccTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGt                                                 <  1:898752/150‑1 (MQ=255)                                                                                             ccTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTc                                                <  2:299183/151‑1 (MQ=255)                                                                                                gACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGAt                                             >  1:943400/1‑151 (MQ=255)                                                                                                gACCGCGCGGGATCAGGGCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACCGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCg                                               <  1:576452/149‑1 (MQ=255)                                                                                                 aCCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATc                                            <  2:272868/151‑1 (MQ=255)                                                                                                 aCCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATc                                            <  2:577448/151‑1 (MQ=255)                                                                                                   cgcgcgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATca                                           >  1:1034517/1‑150 (MQ=255)                                                                                                    gcgcgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATcaac                                         >  2:1008000/1‑151 (MQ=255)                                                                                                    gcgcgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATcaac                                         <  2:654534/151‑1 (MQ=255)                                                                                                      gcgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATcaaca                                        >  2:1072379/1‑150 (MQ=255)                                                                                                      gcgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATcaaca                                        <  1:452791/150‑1 (MQ=255)                                                                                                       cgGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCaacaac                                      <  1:1118393/151‑1 (MQ=255)                                                                                                           aTCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACAcgc                                  <  1:769427/151‑1 (MQ=255)                                                                                                             cAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACAcgcg                                 <  1:887891/150‑1 (MQ=255)                                                                                                              aGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGtt                               >  1:1002529/1‑151 (MQ=255)                                                                                                              aGGTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGtt                               >  1:730871/1‑151 (MQ=255)                                                                                                                gTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTa                             <  2:579896/151‑1 (MQ=255)                                                                                                                gTCAAATTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTa                             <  1:970341/151‑1 (MQ=255)                                                                                                                   aaaTTTGGCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTg                          <  1:973697/151‑1 (MQ=255)                                                                                                                          gCAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTgg                    >  1:406998/1‑150 (MQ=255)                                                                                                                           cAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCt                  <  2:50561/151‑1 (MQ=255)                                                                                                                           cAATAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCt                  <  1:545789/151‑1 (MQ=255)                                                                                                                             aTAGAGGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTg                 >  1:608801/1‑150 (MQ=255)                                                                                                                               agagGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGAt               <  1:409925/150‑1 (MQ=255)                                                                                                                                gagGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATc              >  1:863226/1‑150 (MQ=255)                                                                                                                                gagGTCGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCAATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATc              >  2:882963/1‑150 (MQ=255)                                                                                                                                     cGGTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCAccgcc        >  2:540069/1‑151 (MQ=255)                                                                                                                                      ggTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCAccgccg       >  2:218876/1‑151 (MQ=255)                                                                                                                                       gTTTAGCCGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACcgccgc      >  2:227206/1‑151 (MQ=255)                                                                                                                                           agccgcGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAgg  >  1:815543/1‑151 (MQ=255)                                                                                                                                           agccgcGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAgg  >  2:251153/1‑151 (MQ=255)                                                                                                                                              |                                                                                                                                                    ATCAATCAACACCAGGTCATAGTTTTTACTCGCCCAGTTCACCAGTTCGGCAAAGCGTTCGCTCATCAACAGTTCAGAAGGATTTGGCGGCACCTGACCGCGCGGGATCAGGTCAAATTTGGCAATAGAGGTCGGTTTAGCAGCGGTAGTAATATCGCCCTGACCAATCAGAATTTCCGACAGGCCATTAACGTTATTAGTGCCCAACAGCTCGTGGGTGTAGCCTTTGCGCATATCGCAGTCGATCAACAACACGCGTTTATTGGTCTGGCTGATCACCGCCGCCAGG  >  NC_000913/2133727‑2134015

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