BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,706,692	T→C	E482E (GAA→GAG)	yhjW ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,706,692	0	T	C	97.5%	852.8 / 6.7	282	E482E (GAA→GAG)	yhjW	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (3/4); new base C (126/149); total (129/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 292 total aligned reads displayed.

TCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTGATC  >  NC_000913/3706600‑3706781

   aGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCa                                                                                        <  2:3497345/93‑1 (MQ=255)

   aGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCa                                                                                        <  1:1920618/93‑1 (MQ=255)

   aGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCa                                                                                        <  2:3905242/93‑1 (MQ=255)

    gATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCAt                                                                                       >  1:3071941/1‑93 (MQ=255)

       atTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTaa                                                                                    >  1:224517/1‑93 (MQ=255)

       atTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTaa                                                                                    >  1:1871398/1‑93 (MQ=255)

        tttATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAAt                                                                                   <  1:31546/93‑1 (MQ=255)

        tttATCTGACATCCAGACCATCATCGGCACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCAGTAAt                                                                                   >  1:4681593/1‑93 (MQ=255)

         ttATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAAtt                                                                                  <  1:1179870/93‑1 (MQ=255)

          tATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTg                                                                                 >  1:4581841/1‑93 (MQ=255)

           aTCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGa                                                                                >  1:4518905/1‑93 (MQ=255)

           aTCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGa                                                                                <  2:1532465/93‑1 (MQ=255)

             cTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACt                                                                              <  1:3594150/93‑1 (MQ=255)

              tGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTc                                                                             <  1:3130236/93‑1 (MQ=255)

              tGACATCCAGACCATCATCGGTACGCGGAACTGCT‑CGGCGGTGCCAGTTCACGCGGCGTTCCGTGCAGGTGCTCGAGCTCATTAATTGACTCa                                                                            >  1:2905033/1‑93 (MQ=255)

               taCATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCa                                                                            <  2:3097369/92‑1 (MQ=255)

               gACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCa                                                                            <  1:2812978/93‑1 (MQ=255)

                aCATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCAc                                                                           <  2:2180000/93‑1 (MQ=255)

                   tCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACGGt                                                                        >  1:2692680/1‑93 (MQ=255)

                   tCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGt                                                                        >  1:2773162/1‑93 (MQ=255)

                   tCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGt                                                                        <  1:3769891/93‑1 (MQ=255)

                   tCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGt                                                                        <  1:808416/93‑1 (MQ=255)

                    ccAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTg                                                                       >  1:4488785/1‑93 (MQ=255)

                     cAGCCCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTgg                                                                      <  1:2525731/93‑1 (MQ=255)

                     cAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTgg                                                                      >  2:940648/1‑93 (MQ=255)

                      aGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGt                                                                     >  2:2067169/1‑93 (MQ=255)

                      aGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGt                                                                     >  2:2770869/1‑93 (MQ=255)

                      aGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGt                                                                     <  1:839058/93‑1 (MQ=255)

                       gACCATCATCGGTACGCGGGACTGCTCCGGCGGTGGCAGGTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTgggt                                                                    >  2:4344250/1‑91 (MQ=255)

                       gACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTc                                                                    >  2:4649562/1‑93 (MQ=255)

                       gACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTACACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCCTGGTc                                                                    >  2:3827920/1‑93 (MQ=255)

                        aCCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCa                                                                   >  1:4442533/1‑93 (MQ=255)

                         ctatcatcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAg                                                                  >  1:3520760/3‑93 (MQ=255)

                          catcatCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGc                                                                 <  2:895474/93‑1 (MQ=255)

                          catcatCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGc                                                                 >  2:2709138/1‑93 (MQ=255)

                            tcatcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTg                                                               >  2:234824/1‑93 (MQ=255)

                            tcatcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTg                                                               >  2:2635308/1‑93 (MQ=255)

                            tcatcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTg                                                               <  2:867064/93‑1 (MQ=255)

                            tcatcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTg                                                               >  1:1057781/1‑93 (MQ=255)

                             catcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGc                                                              <  1:4380198/93‑1 (MQ=255)

                              atcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCg                                                             >  2:3819108/1‑93 (MQ=255)

                               tcGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGt                                                            >  2:3284042/1‑93 (MQ=255)

                                   tACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGaa                                                        <  2:2310374/93‑1 (MQ=255)

                                   tACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGaa                                                        >  1:476674/1‑93 (MQ=255)

                                    aCGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAac                                                       <  1:323047/93‑1 (MQ=255)

                                     cgcgGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAaca                                                      >  2:3608814/1‑93 (MQ=255)

                                       cgGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGTTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAAcaca                                                    >  2:3246221/1‑93 (MQ=255)

                                        ggAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACaa                                                   >  2:2555948/1‑93 (MQ=255)

                                          aaCTGCTCCGGCGGTGCCAGTTCACGCGTCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATc                                                 <  2:3629875/93‑1 (MQ=255)

                                           aCTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCg                                                <  2:3145911/93‑1 (MQ=255)

                                           aCTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCg                                                >  2:1774891/1‑93 (MQ=255)

                                              gCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCtt                                             <  1:3305356/93‑1 (MQ=255)

                                                tCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTc                                           <  2:1142558/93‑1 (MQ=255)

                                                 ccGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCt                                          <  1:1418134/93‑1 (MQ=255)

                                                   ggcggTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTa                                        <  2:2979615/93‑1 (MQ=255)

                                                   ggcggTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTa                                        <  1:3375016/93‑1 (MQ=255)

                                                    gcggTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTAt                                       >  2:3713735/1‑93 (MQ=255)

                                                    gcggTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTAt                                       <  2:1940221/93‑1 (MQ=255)

                                                       gTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATcgc                                    <  2:1838862/93‑1 (MQ=255)

                                                       gTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATcgc                                    >  1:349209/1‑93 (MQ=255)

                                                         gCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGa                                  <  2:4658659/93‑1 (MQ=255)

                                                          ccAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGaa                                 >  1:1897799/1‑93 (MQ=255)

                                                          ccAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGaa                                 <  2:287180/93‑1 (MQ=255)

                                                             gTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCt                              <  1:1028033/93‑1 (MQ=255)

                                                             gTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCt                              >  2:3444515/1‑93 (MQ=255)

                                                             gTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCt                              <  2:4169995/93‑1 (MQ=255)

                                                             gTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCt                              <  2:566292/93‑1 (MQ=255)

                                                             gTTCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCt                              <  2:2600361/93‑1 (MQ=255)

                                                              ttCACGCCGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTg                             <  1:3664531/93‑1 (MQ=255)

                                                               tCACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGa                            <  1:247883/93‑1 (MQ=255)

                                                                cACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGat                           <  1:427671/93‑1 (MQ=255)

                                                                cACGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGat                           <  2:3327677/93‑1 (MQ=255)

                                                                 aCGCGGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatc                          <  2:3989046/93‑1 (MQ=255)

                                                                  cgcgGCGTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatca                         <  2:1124882/93‑1 (MQ=255)

                                                                   gcggcgTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatcaa                        >  1:2405708/1‑93 (MQ=255)

                                                                   gcggcgTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatcaa                        <  2:3358660/93‑1 (MQ=255)

                                                                   gcggcgTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatcaa                        <  1:4005479/93‑1 (MQ=255)

                                                                     ggcgTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatcaatc                      <  1:1829681/93‑1 (MQ=255)

                                                                      gcgTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatcaatca                     <  1:3780671/93‑1 (MQ=255)

                                                                      gcgTGCCGTGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGatcaatca                     >  1:3800260/1‑93 (MQ=255)

                                                                      gcgTGCCGGGCAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGTACCTGatcaatca                     >  2:4119502/1‑93 (MQ=255)

                                                                       cgtgccgtgcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCAc                    <  1:2379812/93‑1 (MQ=255)

                                                                        gtgccgtgcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACg                   >  2:2415524/1‑93 (MQ=255)

                                                                          gccgtgcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCt                 <  1:3916665/93‑1 (MQ=255)

                                                                             gtgcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGga              >  2:4476069/1‑93 (MQ=255)

                                                                             gtgcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGga              <  2:1028284/93‑1 (MQ=255)

                                                                               gcAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGgaga            <  2:4537372/93‑1 (MQ=255)

                                                                               gcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGgaga            <  2:4087122/93‑1 (MQ=255)

                                                                               gcAGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGgaga            <  1:3847982/93‑1 (MQ=255)

                                                                                 aGGTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATa          <  2:141184/93‑1 (MQ=255)

                                                                                  ggTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATaa         <  2:683186/93‑1 (MQ=255)

                                                                                  ggTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATaa         >  1:1389525/1‑93 (MQ=255)

                                                                                   gTGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATaaa        <  2:1007084/93‑1 (MQ=255)

                                                                                    tGCTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAg       <  1:2054935/93‑1 (MQ=255)

                                                                                      cTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTg     <  1:1550920/93‑1 (MQ=255)

                                                                                      cTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTg     >  1:662291/1‑93 (MQ=255)

                                                                                      cTCGCGCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTg     <  2:2788013/93‑1 (MQ=255)

                                                                                        cgcgCTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTGAt   <  1:3043669/93‑1 (MQ=255)

                                                                                         gcgcTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTGATc  >  2:4246951/1‑93 (MQ=255)

                                                                                         gcgcTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTGATc  >  2:1232881/1‑93 (MQ=255)

TCCAGATATTTATCTGACATCCAGACCATCATCGGTACGCGGAACTGCTCCGGCGGTGCCAGTTCACGCGGCGTGCCGTGCAGGTGCTCGCGTTCATTAATTGACTCACCGTGGTCAGCTGCGTAGAACACAATCGCTTTCTTATCGCGAACCTGATCAATCACGCTGGAGATAAAGTGATC  >  NC_000913/3706600‑3706781

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: `‑`