BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,250,165	C→T	R74R (CGC→CGT)	yhaH →	putative cytochrome

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,250,165	0	C	T	98.8%	879.6 / ‑2.9	248	R74R (CGC→CGT)	yhaH	putative cytochrome
Reads supporting (aligned to +/- strand): ref base C (1/2); new base T (77/167); total (78/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 271 total aligned reads displayed.

AGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCA  >  NC_000913/3250075‑3250254

aGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCt                                                                                         <  1:8587008/93‑1 (MQ=255)

aGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCt                                                                                         <  2:5378307/93‑1 (MQ=255)

aGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCt                                                                                         >  2:1825851/1‑93 (MQ=255)

aGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCt                                                                                         >  2:2278977/1‑93 (MQ=255)

 ggCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTg                                                                                        <  2:3791631/93‑1 (MQ=255)

   cTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCa                                                                                      >  2:7955925/1‑93 (MQ=255)

     ggCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACg                                                                                    >  1:4923297/1‑93 (MQ=255)

        aaCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGaca                                                                                 >  2:363363/1‑93 (MQ=255)

        aaCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGaca                                                                                 <  2:4672713/93‑1 (MQ=255)

        aaCGCGCTGGAGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGGGGGCGGGTCAGTTCCGCCGTCTGCACGaca                                                                                 >  1:4456201/1‑93 (MQ=255)

         aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGacac                                                                                <  2:2873141/93‑1 (MQ=255)

         aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGacac                                                                                >  1:7465941/1‑93 (MQ=255)

         aCGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGacac                                                                                <  2:4040129/93‑1 (MQ=255)

          cgcgCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACacc                                                                               >  2:2077497/1‑93 (MQ=255)

           gcgcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccg                                                                              <  2:4233604/93‑1 (MQ=255)

            cgcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccga                                                                             >  2:576962/1‑93 (MQ=255)

             gcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccgac                                                                            <  2:1767612/93‑1 (MQ=255)

             gcTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccgac                                                                            <  1:175464/93‑1 (MQ=255)

              cTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccgacc                                                                           >  1:687899/1‑93 (MQ=255)

              cTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACaccgacc                                                                           >  1:3242028/1‑93 (MQ=255)

                ggtggCGAAGGCATCCTGACGACCATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGc                                                                         <  2:5938493/93‑1 (MQ=255)

                ggtggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGc                                                                         >  1:7650177/1‑93 (MQ=255)

                ggtggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGc                                                                         >  1:6248788/1‑93 (MQ=255)

                 gtggCGAAGGCATCCTGACGACACTCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCt                                                                        <  2:4757130/93‑1 (MQ=255)

                 gtggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCt                                                                        <  1:5565122/93‑1 (MQ=255)

                  tggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTc                                                                       <  2:7690589/93‑1 (MQ=255)

                  tggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTc                                                                       <  2:5138757/93‑1 (MQ=255)

                  tggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTc                                                                       <  1:1908000/93‑1 (MQ=255)

                   ggCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCg                                                                      >  1:3651587/1‑93 (MQ=255)

                    gCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCgg                                                                     >  1:945174/1‑93 (MQ=255)

                    gCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCgg                                                                     <  2:3124991/93‑1 (MQ=255)

                    gCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGCGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCgg                                                                     >  1:1268714/1‑93 (MQ=255)

                     cGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGc                                                                    <  2:8393216/93‑1 (MQ=255)

                       aaGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgt                                                                  <  1:4815312/93‑1 (MQ=255)

                       aaGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCACTCGGCgt                                                                  >  2:3440263/1‑93 (MQ=255)

                         ggCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtgg                                                                >  2:6304699/1‑93 (MQ=255)

                          gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                               >  1:538693/1‑93 (MQ=255)

                          gCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                               <  1:5979780/93‑1 (MQ=255)

                            aTCCTGACGACAATCTACGGTCTTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGtggtgg                                                             >  1:8235110/1‑93 (MQ=255)

                              ccTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGc                                                           <  2:1394868/93‑1 (MQ=255)

                               cTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCa                                                          <  1:3141836/93‑1 (MQ=255)

                               cTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCa                                                          <  1:1570237/93‑1 (MQ=255)

                                 gacgacAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACt                                                        <  2:79251/93‑1 (MQ=255)

                                    gacAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCt                                                     >  1:4736159/1‑93 (MQ=255)

                                     acAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTa                                                    <  2:1349807/93‑1 (MQ=255)

                                      cAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAt                                                   >  2:758056/1‑93 (MQ=255)

                                      cAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAt                                                   <  2:2820997/93‑1 (MQ=255)

                                       aaTCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTAtt                                                  <  1:703313/93‑1 (MQ=255)

                                         tCTACGTTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCt                                                <  1:3555863/93‑1 (MQ=255)

                                           tACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCAGACGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTa                                              >  2:7062508/1‑93 (MQ=255)

                                            aCGGTAGTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTaa                                             <  1:7221220/93‑1 (MQ=255)

                                              ggTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATc                                           >  1:5540906/1‑93 (MQ=255)

                                              cgTATTCTGGTGTTTTTTCCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATc                                           <  2:1150704/92‑1 (MQ=255)

                                               gTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATcc                                          <  2:4983220/93‑1 (MQ=255)

                                               gTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATcc                                          >  1:4175463/1‑93 (MQ=255)

                                                  ttCTGGTGTTTTTTCCGTGGTGGGCGGTTCAGTTCCTCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGt                                       <  2:7449950/93‑1 (MQ=255)

                                                  ttCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGt                                       <  1:1661026/93‑1 (MQ=255)

                                                    cTGGTTTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGttt                                     <  2:5991786/93‑1 (MQ=255)

                                                    cTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGttt                                     <  1:6935463/93‑1 (MQ=255)

                                                       gtgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATc                                  >  1:1882562/1‑93 (MQ=255)

                                                        tgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCg                                 <  2:262509/93‑1 (MQ=255)

                                                        tgtTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCg                                 <  1:629305/93‑1 (MQ=255)

                                                          tttttACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                               <  2:2708325/61‑1 (MQ=255)

                                                          tttttACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCgtggt                                                               >  1:2708325/1‑61 (MQ=255)

                                                          tttttACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggc                               <  1:7638787/93‑1 (MQ=255)

                                                          tttttACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggc                               <  2:7032785/93‑1 (MQ=255)

                                                           ttttACCGTGGTGGGCGGTTCAGTTCCGCCGGCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggct                              <  2:6603442/93‑1 (MQ=255)

                                                            tttACCGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggctg                             <  2:8350247/93‑1 (MQ=255)

                                                             ttACCGCGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCggctgg                            <  1:6742403/93‑1 (MQ=255)

                                                                ccGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtt                         >  2:900865/1‑93 (MQ=255)

                                                                 cGTGGTGGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCtta                        <  2:1119583/93‑1 (MQ=255)

                                                                  gtggtgGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCTCTGCTATTCTTAATCCCGTTTATCGGCTGGCttat                       <  2:8719026/93‑1 (MQ=255)

                                                                  gtggtgGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttat                       <  1:3486177/93‑1 (MQ=255)

                                                                  gtggtgGGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttat                       <  2:6030694/93‑1 (MQ=255)

                                                                   tggtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatt                      <  1:1825851/93‑1 (MQ=255)

                                                                   tggtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatt                      <  1:2278977/93‑1 (MQ=255)

                                                                   aggtggGCGGTTCAGTTCCGCCGTATGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatt                      <  1:3813822/92‑1 (MQ=255)

                                                                    ggtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCttatta                     <  1:4370961/93‑1 (MQ=255)

                                                                     gtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTtattat                    <  2:6971558/93‑1 (MQ=255)

                                                                     gtggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTtattat                    <  1:8776056/93‑1 (MQ=255)

                                                                      tggGCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatc                   <  1:1348436/93‑1 (MQ=255)

                                                                      ggggCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatc                   <  1:7661677/92‑1 (MQ=255)

                                                                       gggCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatca                  <  1:2470984/93‑1 (MQ=255)

                                                                       gggCGGTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTatca                  <  1:7918334/93‑1 (MQ=255)

                                                                            gTTCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGtt             >  2:1711812/1‑92 (MQ=255)

                                                                            gTTCAGTGCCGCCGTCTGCCTGCCACCGCCCGCTCGGCTTGTTTGGCACTTCTATTCTTAATCTCGTTTATCGATTGGCTTATTATCATCGTc             <  2:5503556/93‑1 (MQ=255)

                                                                             ttCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCt            >  1:7365331/1‑93 (MQ=255)

                                                                              tCAGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCtt           <  1:5863938/93‑1 (MQ=255)

                                                                                aGTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTa         <  2:6205734/93‑1 (MQ=255)

                                                                                 gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTaa        <  2:3297754/93‑1 (MQ=255)

                                                                                 gTTCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTaa        >  1:2972966/1‑93 (MQ=255)

                                                                                  ttCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAAc       <  1:5143969/93‑1 (MQ=255)

                                                                                   tCCGCCGTCTTCACGACACCGACCGGTCGGCGGGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACt      <  2:5537509/93‑1 (MQ=255)

                                                                                   tCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACt      <  2:7415611/93‑1 (MQ=255)

                                                                                   tCCGCCGTCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACt      <  2:2470271/93‑1 (MQ=255)

                                                                                     cgccgtCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGc    <  2:1570190/93‑1 (MQ=255)

                                                                                      gccgtcTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGcc   <  1:4018905/93‑1 (MQ=255)

                                                                                      gccgtcTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGcc   <  1:6484373/93‑1 (MQ=255)

                                                                                      gccgtcTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGcc   <  2:6176878/93‑1 (MQ=255)

                                                                                       ccgtcTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCa  <  2:8072837/93‑1 (MQ=255)

AGGCTGGCAACGCGCTGGTGGCGAAGGCATCCTGACGACAATCTACGGTATTCTGGTGTTTTTACCGTGGTGGGCGGTTCAGTTCCGCCGCCTGCACGACACCGACCGCTCGGCGTGGTGGGCACTGCTATTCTTAATCCCGTTTATCGGCTGGCTTATTATCATCGTCTTTAACTGCCA  >  NC_000913/3250075‑3250254

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