BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,197,178	T→C	V35V (GTA→GTG)	glnE ←	GlnE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,197,178	0	T	C	98.7%	735.9 / ‑5.8	236	V35V (GTA→GTG)	glnE	GlnE
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (122/110); minor base A (3/0); total (126/110)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

TTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGA  >  NC_000913/3197088‑3197267

ttCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCAc                                                                                         >  1:5081300/1‑93 (MQ=255)

ttCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCAc                                                                                         >  1:4486169/1‑93 (MQ=255)

    tcgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGAc                                                                                     <  2:730839/93‑1 (MQ=255)

     cgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACt                                                                                    >  1:416186/1‑93 (MQ=255)

     cgGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCAC‑CGTTCCTCCAAAAAAATCACAAAAAATAAACACTAACt                                                                                    >  1:449635/1‑93 (MQ=255)

      ggCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACtt                                                                                   <  2:5559216/93‑1 (MQ=255)

      ggCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGAACTGACtt                                                                                   >  1:136283/1‑93 (MQ=255)

       gCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTc                                                                                  >  2:4265261/1‑93 (MQ=255)

       gCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTc                                                                                  >  2:1598828/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGcc                                                                               >  2:3558241/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGcc                                                                               >  2:4900980/1‑93 (MQ=255)

          tGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGcc                                                                               >  2:519615/1‑93 (MQ=255)

           gcggcgGTGGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCt                                                                              <  1:5356347/93‑1 (MQ=255)

            cggcggTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCtg                                                                             <  1:2436927/93‑1 (MQ=255)

               cggTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgc                                                                          <  1:5976795/93‑1 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGGACTGACTTCGCCTGTgcgc                                                                        >  1:4117057/1‑93 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  2:2728562/1‑93 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  1:2295072/1‑93 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  1:1304516/1‑93 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  2:1443038/1‑93 (MQ=255)

                 gTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTgcgc                                                                        >  2:4227955/1‑93 (MQ=255)

                  ttGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCt                                                                       >  1:981491/1‑93 (MQ=255)

                   tGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTa                                                                      <  1:4728413/93‑1 (MQ=255)

                      cTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGt                                                                   <  2:4056967/93‑1 (MQ=255)

                       tttGCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTg                                                                  <  1:3936793/93‑1 (MQ=255)

                          ccAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGAtt                                                               <  1:3270089/93‑1 (MQ=255)

                             gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                            >  2:1824935/1‑93 (MQ=255)

                             gTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCt                                                            >  2:5892646/1‑93 (MQ=255)

                              ttCCGTCAGCCACTCTTGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAATCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           >  2:2729916/1‑93 (MQ=255)

                              ttCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           <  1:4885868/93‑1 (MQ=255)

                              ttCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTc                                                           >  2:4149667/1‑93 (MQ=255)

                                ccGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCgg                                                         >  2:6046815/1‑93 (MQ=255)

                                  gTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCt                                                       >  1:1212590/1‑93 (MQ=255)

                                     aGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTaaa                                                    <  2:4298586/93‑1 (MQ=255)

                                     aGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTaaa                                                    <  2:2082946/93‑1 (MQ=255)

                                     aGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTaaa                                                    <  1:555841/93‑1 (MQ=255)

                                      gCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAg                                                   >  2:5910091/1‑93 (MQ=255)

                                       ccACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAgg                                                  <  1:3208075/93‑1 (MQ=255)

                                        cACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGc                                                 >  2:2853037/1‑93 (MQ=255)

                                        cACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGc                                                 <  1:5648817/93‑1 (MQ=255)

                                         aCTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGct                                                >  2:810893/1‑93 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctt                                               >  1:4221370/1‑92 (MQ=255)

                                          ctctGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACGAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctc                                               <  2:4005572/93‑1 (MQ=255)

                                           tctGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                              <  2:2290479/93‑1 (MQ=255)

                                           tctGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                              <  1:1501437/93‑1 (MQ=255)

                                           tctGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGctct                                              <  2:2691539/93‑1 (MQ=255)

                                            ctGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTg                                             <  1:5900426/93‑1 (MQ=255)

                                             tGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTgg                                            <  1:1389267/93‑1 (MQ=255)

                                             tGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTgg                                            >  2:1343800/1‑93 (MQ=255)

                                                aTGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAg                                         <  2:2895921/93‑1 (MQ=255)

                                                 tGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGc                                        <  2:3805771/93‑1 (MQ=255)

                                                  gcgcAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGcc                                       >  2:704606/1‑93 (MQ=255)

                                                    gcAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGc                                     >  1:1038321/1‑93 (MQ=255)

                                                     cAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCt                                    >  1:342309/1‑93 (MQ=255)

                                                      aaTCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTc                                   <  1:1002657/93‑1 (MQ=255)

                                                      aaTCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTc                                   <  2:661072/93‑1 (MQ=255)

                                                      aaTCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTc                                   <  1:6022553/93‑1 (MQ=255)

                                                       aTCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTca                                  >  1:2662912/1‑93 (MQ=255)

                                                         cACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaac                                <  1:726464/93‑1 (MQ=255)

                                                         cACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTcaac                                >  1:4917464/1‑93 (MQ=255)

                                                             cTGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACg                            >  2:4476553/1‑93 (MQ=255)

                                                              tGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACgg                           >  1:2254248/1‑93 (MQ=255)

                                                              tGTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACgg                           >  2:2486924/1‑93 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCACCAACGGt                          >  1:1125887/1‑93 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          <  1:1568744/93‑1 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          >  1:3858487/1‑93 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          >  1:4283227/1‑93 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          <  2:5427195/93‑1 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          <  2:2558281/93‑1 (MQ=255)

                                                               gTCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGt                          >  1:1792107/1‑93 (MQ=255)

                                                                tCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                         >  2:5596171/1‑93 (MQ=255)

                                                                tCCTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTc                         >  1:3604950/1‑93 (MQ=255)

                                                                 ccTGC‑ACAAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCt                        <  2:120263/93‑1 (MQ=255)

                                                                     c‑acaAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTTATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCa                    >  2:5749085/1‑93 (MQ=255)

                                                                       ccaAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg                   <  1:550712/92‑1 (MQ=255)

                                                                       acaAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAg                   <  2:5269973/93‑1 (MQ=255)

                                                                        caAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGt                  >  2:1140126/1‑93 (MQ=255)

                                                                        caAAATCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGt                  <  1:5692306/93‑1 (MQ=255)

                                                                         aaaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTa                 >  1:692497/1‑93 (MQ=255)

                                                                         aaaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTa                 <  1:2052558/93‑1 (MQ=255)

                                                                          aaaTCACTAAAAGTAAGCCCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc                <  1:323496/93‑1 (MQ=255)

                                                                          aaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc                >  1:2793790/1‑93 (MQ=255)

                                                                          aaaTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTAc                >  1:2968291/1‑93 (MQ=255)

                                                                            aTCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActg              >  1:4053696/1‑93 (MQ=255)

                                                                             tCACTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgc             >  1:2074413/1‑93 (MQ=255)

                                                                              cACTAAAAGTAAGCACTGACTTCTCCTGTGCGTTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgct            <  1:2840566/93‑1 (MQ=255)

                                                                               aCTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg           >  2:1713164/1‑93 (MQ=255)

                                                                               aCTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg           >  1:5794228/1‑93 (MQ=255)

                                                                               aCTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTActgctg           >  1:1614593/1‑93 (MQ=255)

                                                                                cTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGt          >  2:4698114/1‑93 (MQ=255)

                                                                                cTAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGt          <  2:5006104/93‑1 (MQ=255)

                                                                                 tAAAAGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTa         >  2:6128371/1‑93 (MQ=255)

                                                                                  aaaaGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTaa        >  1:952169/1‑93 (MQ=255)

                                                                                  aaaaGTAAGCACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTaa        <  1:4788472/93‑1 (MQ=255)

                                                                                     agtaagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACgg     <  1:4155792/93‑1 (MQ=255)

                                                                                     agtaagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACgg     <  2:5044059/93‑1 (MQ=255)

                                                                                       taagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTg   <  2:1896450/93‑1 (MQ=255)

                                                                                       taagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTg   >  2:1963066/1‑93 (MQ=255)

                                                                                       taagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTg   >  1:2259347/1‑93 (MQ=255)

                                                                                        aagcaCTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGa  <  2:3292600/93‑1 (MQ=255)

TTCGTCGGCCTGCGGCGGTTGGCTTTCCAGTTCCGTCAGCCACTCTGGATGCGCAATCACGCTGTCCTGC‑ACAAAATCACTAAAAGTAAGTACTGACTTCGCCTGTGCGCTAAGTGATTCCTCGGCTAAAGGCTCTGGCAGCCGCTCAACAACGGTCTGCCAGTACTGCTGTAACGGTGA  >  NC_000913/3197088‑3197267

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