BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,904,576	G→A	P170L (CCC→CTC)	yebH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,904,576	0	G	A	94.5%	463.1 / 8.5	145	P170L (CCC→CTC)	yebH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/6); new base A (71/66); total (73/72)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-02

Only 100 of 147 total aligned reads displayed.

GGTGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAAC  >  NC_000913/1904486‑1904663

ggTGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGc                                                                                       <  2:1447421/93‑1 (MQ=255)

 gTGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCt                                                                                      <  1:3665355/93‑1 (MQ=255)

 gTGCATGAAAATGGACTTCATTGTAAATCAGCCCCTTCAGCTCCAGCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCt                                                                                      <  1:3665356/93‑1 (MQ=255)

     aTGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCac                                                                                  <  1:4520258/93‑1 (MQ=255)

     aTGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCac                                                                                  >  2:2295575/1‑93 (MQ=255)

     aTGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCac                                                                                  >  1:1177639/1‑93 (MQ=255)

     aTGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGCTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCac                                                                                  <  2:1894627/93‑1 (MQ=255)

       gAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCacta                                                                                <  1:3508603/93‑1 (MQ=255)

       gAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCacta                                                                                >  2:2922903/1‑93 (MQ=255)

        aaGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCactac                                                                               >  2:1170349/1‑93 (MQ=255)

         aGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCactact                                                                              >  2:1329181/1‑93 (MQ=255)

         aGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCactact                                                                              >  2:1907060/1‑93 (MQ=255)

           aTGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCg                                                                            >  2:3983144/1‑93 (MQ=255)

            tGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCGTCACTACTCGt                                                                           <  2:2748860/93‑1 (MQ=255)

            tGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGt                                                                           >  2:3567823/1‑93 (MQ=255)

                 ttCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTaa                                                                      >  1:3051874/1‑93 (MQ=255)

                  tCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAt                                                                     >  1:2526199/1‑93 (MQ=255)

                  tCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAt                                                                     >  2:931511/1‑93 (MQ=255)

                    aTTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAttt                                                                   >  2:1060427/1‑92 (MQ=255)

                    aTTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAttc                                                                   >  2:1060428/1‑93 (MQ=255)

                    aTTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAttc                                                                   >  2:1023220/1‑93 (MQ=255)

                      tGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAttctt                                                                 <  1:4712277/93‑1 (MQ=255)

                       gTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAAttcttc                                                                <  1:3937580/93‑1 (MQ=255)

                        tAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATtcttct                                                               <  1:1603381/93‑1 (MQ=255)

                            tCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCtt                                                           <  2:2463084/93‑1 (MQ=255)

                             cAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCttt                                                          <  2:4145185/93‑1 (MQ=255)

                              aGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCtttc                                                         >  2:1218183/1‑92 (MQ=255)

                              aGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTa                                                         <  2:2658838/93‑1 (MQ=255)

                               gCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTac                                                        <  2:2796872/93‑1 (MQ=255)

                                ccccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTata                                                       >  1:3513136/1‑91 (MQ=255)

                                ccccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTaca                                                       <  2:533937/93‑1 (MQ=255)

                                ccccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGATAATTCTTCTGCTTTaca                                                       <  1:3052285/93‑1 (MQ=255)

                                 cccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac                                                      >  2:4615782/1‑93 (MQ=255)

                                 cccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac                                                      >  2:4700618/1‑93 (MQ=255)

                                 cccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac                                                      <  1:3755167/93‑1 (MQ=255)

                                 cccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac                                                      >  1:1093172/1‑93 (MQ=255)

                                 cccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac                                                      <  2:3102988/93‑1 (MQ=255)

                                 cccTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTacac                                                      <  1:125/93‑1 (MQ=255)

                                   cTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACgg                                                    >  1:2918413/1‑93 (MQ=255)

                                   cTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACgg                                                    >  2:4254674/1‑93 (MQ=255)

                                    ttCAGCTCCATCAAATGTCGCGGTCCCGGCGTGTCAGTAACGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGgc                                                   <  1:3208404/93‑1 (MQ=255)

                                    ttCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGgc                                                   >  2:4434491/1‑93 (MQ=255)

                                    ttCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGgc                                                   >  2:1744645/1‑93 (MQ=255)

                                    ttCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGgc                                                   <  1:1102514/93‑1 (MQ=255)

                                     tCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGgcg                                                  >  2:3173967/1‑93 (MQ=255)

                                      cAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGgcgc                                                 <  1:2577046/93‑1 (MQ=255)

                                       cgCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGcgcg                                                <  1:174884/92‑1 (MQ=255)

                                       aGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGcgcg                                                >  1:2288351/1‑93 (MQ=255)

                                       aGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGcgcg                                                >  2:636534/1‑93 (MQ=255)

                                       aGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGcgcg                                                >  2:2968362/1‑93 (MQ=255)

                                        gCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGt                                               >  2:229619/1‑93 (MQ=255)

                                        gCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTCCTTTACACGGCGCGt                                               >  1:2558610/1‑93 (MQ=255)

                                          tCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTAGACGGCGCGTaa                                             <  2:2749428/93‑1 (MQ=255)

                                          tCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTaa                                             >  1:1966613/1‑93 (MQ=255)

                                            cATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGCGCTAATTCTTCTGCTTTACACGGCGCGTAAAt                                           <  2:3440495/93‑1 (MQ=255)

                                             aTCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATa                                          <  1:3393018/93‑1 (MQ=255)

                                             aTCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGACCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATa                                          >  2:1172092/1‑93 (MQ=255)

                                              tCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATAc                                         >  1:4292133/1‑93 (MQ=255)

                                               cAAATGTCGCGGTGCCGGCGTGTGAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACg                                        <  2:1648132/93‑1 (MQ=255)

                                                aaaTGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGt                                       >  1:4105997/1‑93 (MQ=255)

                                                 aaTGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGta                                      >  1:3570025/1‑93 (MQ=255)

                                                  aTGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGtat                                     <  1:1215138/93‑1 (MQ=255)

                                                   tGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGtatt                                    >  2:4273657/1‑93 (MQ=255)

                                                    gTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGtatta                                   <  2:209417/93‑1 (MQ=255)

                                                     gcgcgGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTTCTCGTGCTAATTCTTCTGCTTTACACGGCTCGTAAATACGtattat                                  <  1:2872862/92‑1 (MQ=255)

                                                       gcgGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATgg                                <  1:3565886/93‑1 (MQ=255)

                                                        cgGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGc                               >  1:311466/1‑93 (MQ=255)

                                                             ccGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCa                          >  2:854707/1‑93 (MQ=255)

                                                              cGGCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCAt                         >  1:3733034/1‑93 (MQ=255)

                                                               ggCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATa                        <  2:2831059/93‑1 (MQ=255)

                                                               ggCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATa                        >  1:1314003/1‑93 (MQ=255)

                                                                gCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATAc                       <  2:84552/93‑1 (MQ=255)

                                                                gCGTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATAc                       <  1:4517864/93‑1 (MQ=255)

                                                                  gTGGCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTg                     >  2:2138787/1‑93 (MQ=255)

                                                                    ggCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGt                   >  1:4307627/1‑93 (MQ=255)

                                                                    ggCAGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCGTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGt                   <  2:1291842/93‑1 (MQ=255)

                                                                       aGTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATc                <  2:3717512/93‑1 (MQ=255)

                                                                        gTAATGCCCCAGACGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGAAGTCCGTTCTAGTTTCg               <  2:4227062/93‑1 (MQ=255)

                                                                        gTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCg               <  1:2145564/93‑1 (MQ=255)

                                                                        gTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCg               <  2:3747076/93‑1 (MQ=255)

                                                                        gTAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCg               >  2:1465251/1‑93 (MQ=255)

                                                                         tAATGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCgg              <  2:4124460/93‑1 (MQ=255)

                                                                          aaTGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGa             >  1:3396246/1‑93 (MQ=255)

                                                                           aTGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGaa            >  1:1849268/1‑93 (MQ=255)

                                                                           aTGACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGaa            >  1:3639070/1‑93 (MQ=255)

                                                                             gACCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGaaaa          <  2:4514086/93‑1 (MQ=255)

                                                                              gcccAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAAc         <  1:2110079/92‑1 (MQ=255)

                                                                              ccccAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAAc         <  1:1732542/92‑1 (MQ=255)

                                                                              aCCCAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAAc         >  2:3515159/1‑93 (MQ=255)

                                                                               cccAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCTTACTGGTATCGGAAAACg        >  1:3011864/1‑93 (MQ=255)

                                                                               cccAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACg        <  2:1513834/93‑1 (MQ=255)

                                                                                ccAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACgg       <  2:3875100/93‑1 (MQ=255)

                                                                                 cAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGc      >  2:964745/1‑93 (MQ=255)

                                                                                 cAGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGc      <  1:2202411/93‑1 (MQ=255)

                                                                                  aGCCGCCGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCa     >  2:1733695/1‑93 (MQ=255)

                                                                                   gccgccGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCaa    >  1:2875324/1‑93 (MQ=255)

                                                                                   gccgccGAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCaa    >  2:2700850/1‑93 (MQ=255)

                                                                                    ccgccgAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCaaa   <  2:4524053/93‑1 (MQ=255)

                                                                                    ccgccgAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCaaa   <  2:4025033/93‑1 (MQ=255)

                                                                                     cgccgAGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAAc  <  2:4174020/93‑1 (MQ=255)

GGTGCATGAAGATGTACTTCATTGTAAATCAGCCCCTTCAGCTCCATCAAATGTCGCGGTCCCGGCGTGGCAGTAATGACCCAGCCGCCGGGCTTCACTACTCGTGCTAATTCTTCTGCTTTACACGGCGCGTAAATACGTATTATGGCGTCCATACTGGTATCGGAAAACGGCAAAC  >  NC_000913/1904486‑1904663

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