BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	329,174	A→G	Y163C (TAC→TGC)	betT →	high‑affinity choline transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	329,174	0	A	G	100.0%	745.5 / NA	230	Y163C (TAC→TGC)	betT	high‑affinity choline transport
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (104/126); total (104/126)

Only 100 of 239 total aligned reads displayed.

AGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329084‑329263

aGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCt                                                                                         <  1:1349374/93‑1 (MQ=255)

 ggCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCtt                                                                                        <  2:788271/93‑1 (MQ=255)

  gCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCttt                                                                                       >  2:4614672/1‑93 (MQ=255)

      cgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCt                                                                                   <  1:4489346/93‑1 (MQ=255)

      cgcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCt                                                                                   <  1:3951467/93‑1 (MQ=255)

       gcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTa                                                                                  <  2:1969521/93‑1 (MQ=255)

       gcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTa                                                                                  <  1:5044402/93‑1 (MQ=255)

       gcgTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTa                                                                                  <  1:1319749/93‑1 (MQ=255)

         gTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATc                                                                                >  1:4047123/1‑93 (MQ=255)

         gTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATc                                                                                <  1:1646246/93‑1 (MQ=255)

            aGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGtt                                                                             >  2:4200503/1‑93 (MQ=255)

             ggCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGCTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTta                                                                            >  2:3575256/1‑93 (MQ=255)

              gCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTtat                                                                           <  1:1531706/93‑1 (MQ=255)

                gATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATaa                                                                         <  2:2953006/93‑1 (MQ=255)

                 aTGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAt                                                                        >  1:2760124/1‑93 (MQ=255)

                 aTGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAt                                                                        <  1:5835202/93‑1 (MQ=255)

                  tGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAtt                                                                       >  2:3894900/1‑93 (MQ=255)

                  tGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAAtt                                                                       <  1:6035453/93‑1 (MQ=255)

                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 >  2:5191321/1‑93 (MQ=255)

                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 <  2:1931125/93‑1 (MQ=255)

                        ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGc                                                                 <  2:5137064/93‑1 (MQ=255)

                          aCGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTc                                                               <  2:1655891/93‑1 (MQ=255)

                           cGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCa                                                              <  1:4357950/93‑1 (MQ=255)

                             cTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATAGTTATAATTTGCCGCTCAcc                                                            <  2:3202627/93‑1 (MQ=255)

                              tGTTTCGCTACTGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCa                                                           <  1:2631487/93‑1 (MQ=255)

                              tGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCa                                                           <  1:6693105/93‑1 (MQ=255)

                               gTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCAt                                                          >  1:6194367/1‑93 (MQ=255)

                                gttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATc                                                         <  1:2469509/92‑1 (MQ=255)

                                  tCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCg                                                       <  2:4369655/93‑1 (MQ=255)

                                   cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCCCCATCCGc                                                      <  2:4842587/93‑1 (MQ=255)

                                   cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGc                                                      <  1:3347541/93‑1 (MQ=255)

                                   cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGc                                                      >  1:886155/1‑93 (MQ=255)

                                   cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGc                                                      <  1:56982/93‑1 (MQ=255)

                                   cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGc                                                      >  2:5317634/1‑93 (MQ=255)

                                    aCTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCt                                                     >  2:2176309/1‑93 (MQ=255)

                                     cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  2:5885249/93‑1 (MQ=255)

                                     cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  2:5817314/93‑1 (MQ=255)

                                     cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  1:3010049/93‑1 (MQ=255)

                                     cTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTc                                                    <  1:1196247/93‑1 (MQ=255)

                                       aCGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCgg                                                  >  1:4858614/1‑93 (MQ=255)

                                       aCGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCgg                                                  <  2:52199/93‑1 (MQ=255)

                                       aCGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCgg                                                  >  1:3006698/1‑93 (MQ=255)

                                         ggCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcg                                                >  2:654928/1‑93 (MQ=255)

                                         ggCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcg                                                <  1:6115254/93‑1 (MQ=255)

                                          gCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcgc                                               <  1:4153611/93‑1 (MQ=255)

                                          gCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcgc                                               >  1:3199515/1‑93 (MQ=255)

                                          gCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcgc                                               <  2:1689336/93‑1 (MQ=255)

                                          gCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGgcgc                                               <  2:6724783/93‑1 (MQ=255)

                                           cTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCt                                              <  2:4940586/93‑1 (MQ=255)

                                           cTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCt                                              >  1:406320/1‑93 (MQ=255)

                                            ttAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTg                                             <  2:1045438/93‑1 (MQ=255)

                                             tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCCCCATCCGCTCGGCGCTGt                                            <  1:2956920/93‑1 (MQ=255)

                                             tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGt                                            >  2:1196583/1‑93 (MQ=255)

                                             tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGt                                            >  2:3336626/1‑93 (MQ=255)

                                               aCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAc                                          >  2:2427088/1‑93 (MQ=255)

                                                ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  2:4061593/93‑1 (MQ=255)

                                                ccGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAcc                                         <  1:6412454/93‑1 (MQ=255)

                                                 cGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAccc                                        <  1:1057689/93‑1 (MQ=255)

                                                  ggCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCg                                       >  1:1419410/1‑93 (MQ=255)

                                                   gCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGa                                      <  2:2374205/93‑1 (MQ=255)

                                                      ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   <  1:2138541/93‑1 (MQ=255)

                                                      ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   <  2:142870/93‑1 (MQ=255)

                                                      ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCt                                   >  1:2306437/1‑93 (MQ=255)

                                                        tCGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTc                                 >  2:398245/1‑93 (MQ=255)

                                                         cGATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCg                                >  2:5529391/1‑93 (MQ=255)

                                                          gATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCgg                               >  1:5207047/1‑93 (MQ=255)

                                                          gATGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCgg                               >  2:281573/1‑93 (MQ=255)

                                                           aTGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGt                              >  2:3755979/1‑93 (MQ=255)

                                                           aTGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGt                              >  1:5439678/1‑93 (MQ=255)

                                                            tGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTa                             >  2:2580243/1‑93 (MQ=255)

                                                            tGTATGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTa                             <  2:674536/93‑1 (MQ=255)

                                                               aTGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          <  1:752896/93‑1 (MQ=255)

                                                               aTGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          <  2:3665637/93‑1 (MQ=255)

                                                               aTGCGCTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAAc                          <  1:3508580/93‑1 (MQ=255)

                                                                  cgcgGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCTCTCGGCGCTGTACCCGATCTTCGGTAAACGGa                       <  1:6658535/93‑1 (MQ=255)

                                                                   gcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAt                      <  1:2741435/93‑1 (MQ=255)

                                                                   gcTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTAACCGATCTTCGGTAAACGGAt                      >  1:1268363/1‑93 (MQ=255)

                                                                    cTGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGAtt                     <  2:4768869/93‑1 (MQ=255)

                                                                     tGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTa                    >  2:3928845/1‑93 (MQ=255)

                                                                     tGATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTa                    >  1:1458689/1‑93 (MQ=255)

                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   >  1:3471323/1‑93 (MQ=255)

                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   >  2:2991225/1‑93 (MQ=255)

                                                                      gATGGGCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTaa                   <  2:4548821/93‑1 (MQ=255)

                                                                         gggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACgg                >  1:4718947/1‑93 (MQ=255)

                                                                          ggCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACggg               >  2:1500208/1‑93 (MQ=255)

                                                                           gCATGGCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGc              <  2:6269679/93‑1 (MQ=255)

                                                                               ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          >  2:4119281/1‑93 (MQ=255)

                                                                               ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          >  2:1813514/1‑93 (MQ=255)

                                                                               ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          <  1:3833242/93‑1 (MQ=255)

                                                                               ggCGCTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGAt          <  1:2562942/93‑1 (MQ=255)

                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  2:1616424/93‑1 (MQ=255)

                                                                                gcgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATa         <  2:1578419/93‑1 (MQ=255)

                                                                                 cgcTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAg        >  2:3419558/1‑93 (MQ=255)

                                                                                   cTCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGt      <  2:6017107/93‑1 (MQ=255)

                                                                                    tCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     >  1:1721856/1‑93 (MQ=255)

                                                                                    tCGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTc     <  2:883097/93‑1 (MQ=255)

                                                                                     cGGATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCa    >  2:3873763/1‑93 (MQ=255)

                                                                                      ggATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCAc   >  1:191360/1‑93 (MQ=255)

                                                                                      ggATGCTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGCCAc   >  2:207444/1‑93 (MQ=255)

                                                                                       gATGCTTTAGCTATCTTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACt  <  2:6201833/93‑1 (MQ=255)

AGGCCGCGCGTCAGGCGATGGTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACT  >  NC_000913/329084‑329263

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