BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,160,251	A→G	T145T (ACT→ACC)	sufI ←	suppressor of ftsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,160,251	0	A	G	99.6%	827.5 / NA	255	T145T (ACT→ACC)	sufI	suppressor of ftsI
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (113/141); total (113/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 266 total aligned reads displayed.

GTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGG  >  NC_000913/3160161‑3160340

gTTGGGGATAGGCACCGACTTCCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGt                                                                                         <  1:823776/93‑1 (MQ=255)

  tGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTAt                                                                                       <  1:8552794/93‑1 (MQ=255)

  tGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTAt                                                                                       <  1:749997/93‑1 (MQ=255)

   ggggATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTAtt                                                                                      >  2:5237809/1‑93 (MQ=255)

      gATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGc                                                                                   >  2:10156/1‑93 (MQ=255)

      gATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGc                                                                                   <  2:1079576/93‑1 (MQ=255)

        tAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGt                                                                                 >  1:8336411/1‑93 (MQ=255)

         aGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTg                                                                                <  1:4809241/93‑1 (MQ=255)

          ggCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGa                                                                               <  1:8076990/93‑1 (MQ=255)

          ggCAGCGACTTGCTGACTTCATCTCCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGa                                                                               <  1:295595/93‑1 (MQ=255)

            cAGCGACTTGCTGACTTCATCTTCCCCCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATa                                                                             <  1:3286528/93‑1 (MQ=255)

            cAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATa                                                                             <  2:5454143/93‑1 (MQ=255)

             aGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATAc                                                                            <  2:877284/93‑1 (MQ=255)

             aGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATAc                                                                            <  1:516245/93‑1 (MQ=255)

              gCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATAcc                                                                           >  2:8657612/1‑93 (MQ=255)

              gCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATAcc                                                                           >  1:4555706/1‑93 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                         <  2:6918929/93‑1 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                         >  2:3225818/1‑93 (MQ=255)

                gACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcac                                                                         >  1:1509150/1‑93 (MQ=255)

                 aCTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACcaca                                                                        >  2:3842436/1‑93 (MQ=255)

                    tGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACagag                                                                     >  2:361414/1‑93 (MQ=255)

                     gCTGACTTCATCTTCCACCATCCACATTCCGGCAAGGCCGTTATAGCCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATCCCACAGAGt                                                                    <  2:6571272/93‑1 (MQ=255)

                     gCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGt                                                                    >  1:3230859/1‑93 (MQ=255)

                     gCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGt                                                                    <  2:3698155/93‑1 (MQ=255)

                     gCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGt                                                                    <  1:5629851/93‑1 (MQ=255)

                      cTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTa                                                                   >  2:3653170/1‑93 (MQ=255)

                       tGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAg                                                                  >  2:9770590/1‑93 (MQ=255)

                        gACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGc                                                                 <  2:3982769/93‑1 (MQ=255)

                           ttCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGc                                                              <  1:345870/93‑1 (MQ=255)

                           ttCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGc                                                              >  1:2335035/1‑93 (MQ=255)

                            tCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCg                                                             >  1:3698569/1‑93 (MQ=255)

                            tCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCg                                                             >  2:2245938/1‑93 (MQ=255)

                            tCATCTGCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCg                                                             <  1:9490496/93‑1 (MQ=255)

                             cATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGt                                                            >  2:6835688/1‑93 (MQ=255)

                             cATCTGCCGCCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGt                                                            <  2:7241036/93‑1 (MQ=255)

                             cATCTCCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGTGGTATTGGCGTGATACCACAGAGTAGCTGCGt                                                            <  1:2391110/93‑1 (MQ=255)

                             cATCTCCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGt                                                            <  1:6405034/93‑1 (MQ=255)

                              aTCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGtt                                                           <  1:1765826/93‑1 (MQ=255)

                              aTCTCCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGGGtt                                                           <  2:9732473/93‑1 (MQ=255)

                               tCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTc                                                          >  1:485301/1‑93 (MQ=255)

                                cTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCt                                                         <  2:4326304/93‑1 (MQ=255)

                                cTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCt                                                         <  1:1311707/93‑1 (MQ=255)

                                 ttCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTg                                                        <  2:2785981/93‑1 (MQ=255)

                                     cccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGa                                                    <  1:493193/92‑1 (MQ=255)

                                     accaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGa                                                    >  2:5114204/1‑93 (MQ=255)

                                     accaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGa                                                    >  2:10023550/1‑93 (MQ=255)

                                     accaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGa                                                    >  2:4751406/1‑93 (MQ=255)

                                      ccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                   <  1:7947143/93‑1 (MQ=255)

                                      ccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                   <  2:6010504/93‑1 (MQ=255)

                                      ccaGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGaa                                                   <  1:7146526/93‑1 (MQ=255)

                                       caGCCACATTCCGGCAAGGCCGTTATAGGCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAAt                                                  <  1:8394913/93‑1 (MQ=255)

                                       caGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAAt                                                  <  1:6626831/93‑1 (MQ=255)

                                        aGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATg                                                 >  2:6365914/1‑93 (MQ=255)

                                        aGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATg                                                 >  2:2390395/1‑93 (MQ=255)

                                        aGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATg                                                 <  2:256095/93‑1 (MQ=255)

                                         gCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGANTAGCTGCGTTCTGGCGAATgg                                                >  1:5853543/1‑93 (MQ=255)

                                         gCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATgg                                                >  1:1390417/1‑93 (MQ=255)

                                         gCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATgg                                                <  2:7541347/93‑1 (MQ=255)

                                         gCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGAGAATgg                                                >  2:3173964/1‑93 (MQ=255)

                                          ccACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATggg                                               <  2:343553/93‑1 (MQ=255)

                                            acaTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCa                                             <  2:1922539/93‑1 (MQ=255)

                                             caTTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAg                                            <  2:1263795/93‑1 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          <  1:557068/93‑1 (MQ=255)

                                               ttCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTa                                          >  2:9287121/1‑93 (MQ=255)

                                                tCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTAc                                         <  1:8025997/93‑1 (MQ=255)

                                                 cccgCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACg                                        <  1:771448/90‑1 (MQ=255)

                                                  cGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACgg                                       <  2:5828755/93‑1 (MQ=255)

                                                     cAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGGGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTg                                    <  1:6048683/93‑1 (MQ=255)

                                                      aaGGCCGTTATAGACCTGCTGAGCCTTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGc                                   >  1:5041742/1‑93 (MQ=255)

                                                        ggCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGccc                                 >  1:6463718/1‑93 (MQ=255)

                                                          ccGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAg                               <  1:8937310/93‑1 (MQ=255)

                                                           cGTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGCGCCCAGt                              >  1:4427923/1‑93 (MQ=255)

                                                            gTTATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTc                             >  2:6804418/1‑93 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                            >  1:2817875/1‑93 (MQ=255)

                                                             ttATAGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCa                            >  2:7191991/1‑93 (MQ=255)

                                                               ataGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGc                          >  2:1145760/1‑93 (MQ=255)

                                                                 aGACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGt                        <  1:3353928/93‑1 (MQ=255)

                                                                  gACCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGtt                       <  1:3103517/93‑1 (MQ=255)

                                                                   aCCTGCTGAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGttt                      >  2:2245156/1‑93 (MQ=255)

                                                                      tgctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGc                   >  1:7891912/1‑93 (MQ=255)

                                                                      tgctgAGCCGAGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGc                   >  2:1403478/1‑93 (MQ=255)

                                                                      tgctgAACCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCACTACGGGTGCCCAGTCAGCGTCTggg                   >  2:2730283/1‑92 (MQ=255)

                                                                       gctgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCg                  >  2:7738126/1‑93 (MQ=255)

                                                                         tgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAc                <  1:8746755/93‑1 (MQ=255)

                                                                         tgAGCCGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAc                >  1:6961395/1‑93 (MQ=255)

                                                                           aGCCGTGCGGTTGGGGGTATTGGGGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTGTGGCGAcat              >  1:292920/1‑93 (MQ=255)

                                                                            gCCGTGCGGTTGGGGGTATTGGCGTGATCCCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatc             <  2:5038927/93‑1 (MQ=255)

                                                                             ccGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatca            <  2:9091135/93‑1 (MQ=255)

                                                                              cGTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGAcatcat           >  1:1625075/1‑93 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTTTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc          >  2:4992510/1‑93 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc          <  2:7139715/93‑1 (MQ=255)

                                                                               gTGCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACatcatc          >  1:2637998/1‑93 (MQ=255)

                                                                                 gctgTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCg        <  1:5239314/90‑1 (MQ=255)

                                                                                 gCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCg        <  1:3942583/93‑1 (MQ=255)

                                                                                 gCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCg        <  1:6326919/93‑1 (MQ=255)

                                                                                 gCGGTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCg        <  1:9580538/93‑1 (MQ=255)

                                                                                    gtgggggGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGc     <  2:455752/90‑1 (MQ=255)

                                                                                    gTTGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGc     <  2:8267529/93‑1 (MQ=255)

                                                                                     ttGGGGGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGcc    >  2:1714703/1‑93 (MQ=255)

                                                                                       ggggGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCgg  >  2:7893456/1‑93 (MQ=255)

GTTGGGGATAGGCAGCGACTTGCTGACTTCATCTTCCACCAGCCACATTCCGGCAAGGCCGTTATAGACCTGCTGAGCCGTGCGGTTGGGAGTATTGGCGTGATACCACAGAGTAGCTGCGTTCTGGCGAATGGGCAGTACGGGTGCCCAGTCAGCGTTTGGCGACATCATCCGTGCCGG  >  NC_000913/3160161‑3160340

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