BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	294,990	T→C	H444R (CAC→CGC)	intF ←	putative phage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	294,990	0	T	C	98.9%	1107.6 / 0.0	359	H444R (CAC→CGC)	intF	putative phage integrase
Reads supporting (aligned to +/- strand): ref base T (3/1); new base C (192/163); total (195/164)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

Only 100 of 383 total aligned reads displayed.

TATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATA  >  NC_000913/294900‑295079

tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         >  1:5038571/1‑93 (MQ=255)

tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         >  1:3562101/1‑93 (MQ=255)

tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         >  1:4061281/1‑93 (MQ=255)

tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         >  1:3067412/1‑93 (MQ=255)

    ttGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAtt                                                                                     >  2:965291/1‑93 (MQ=255)

       aTTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTc                                                                                  <  2:7073861/93‑1 (MQ=255)

       aTTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTc                                                                                  <  1:3708691/93‑1 (MQ=255)

        tttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCg                                                                                 >  2:1210190/1‑93 (MQ=255)

        tttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCg                                                                                 >  1:4610088/1‑93 (MQ=255)

          tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTa                                                                               >  2:4308849/1‑93 (MQ=255)

          tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTa                                                                               >  1:3401188/1‑93 (MQ=255)

          tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTa                                                                               >  1:4828515/1‑93 (MQ=255)

               tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGa                                                                          <  1:794714/93‑1 (MQ=255)

               tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGa                                                                          <  1:453862/93‑1 (MQ=255)

               aggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGa                                                                          <  1:2492303/92‑1 (MQ=255)

                 gAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATc                                                                        >  1:6068531/1‑93 (MQ=255)

                    tCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTaa                                                                     >  2:5054494/1‑93 (MQ=255)

                     cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAAc                                                                    <  1:6246685/93‑1 (MQ=255)

                     cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAAc                                                                    <  1:5935108/93‑1 (MQ=255)

                     cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAAc                                                                    <  2:4400439/93‑1 (MQ=255)

                        cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGga                                                                 <  2:257866/93‑1 (MQ=255)

                        cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGga                                                                 <  1:3922192/93‑1 (MQ=255)

                        cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGga                                                                 <  2:1098994/93‑1 (MQ=255)

                         gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgac                                                                >  2:5656088/1‑93 (MQ=255)

                         gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgac                                                                >  1:3239841/1‑93 (MQ=255)

                         gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgac                                                                <  2:4578890/93‑1 (MQ=255)

                         gCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgac                                                                >  2:2145736/1‑93 (MQ=255)

                           atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacga                                                              <  2:3386834/93‑1 (MQ=255)

                            tatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacgac                                                             >  1:2995828/1‑93 (MQ=255)

                            tatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacgac                                                             >  2:7137124/1‑93 (MQ=255)

                            tatCAACGTTGTTTTTTATGGTAATACCTGATTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTCACGgacgac                                                             >  2:5401358/1‑93 (MQ=255)

                              tCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACgg                                                           <  2:4693362/93‑1 (MQ=255)

                                aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCg                                                         >  2:191531/1‑93 (MQ=255)

                                aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCg                                                         >  2:6721660/1‑93 (MQ=255)

                                 aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGa                                                        >  2:5295753/1‑93 (MQ=255)

                                   gttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                      >  2:6942389/1‑93 (MQ=255)

                                   gttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                      <  1:2846604/93‑1 (MQ=255)

                                    ttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAg                                                     >  1:1441333/1‑93 (MQ=255)

                                    ttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAg                                                     >  1:7718307/1‑93 (MQ=255)

                                     tgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgt                                                    <  2:759275/93‑1 (MQ=255)

                                      ggtttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtg                                                   <  2:1369886/91‑1 (MQ=255)

                                       ttttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtgt                                                  >  2:5752320/1‑93 (MQ=255)

                                          tttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                               >  1:7558314/1‑93 (MQ=255)

                                          tttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                               <  1:3130393/93‑1 (MQ=255)

                                           ttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGttttt                                              >  1:239003/1‑93 (MQ=255)

                                            tagggTAATACCTGTTTCATTTAAGATCCATGTCTCAGTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGGCGGCGATAGTGTTTTTc                                             <  2:283468/90‑1 (MQ=255)

                                            tATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                             >  2:1694175/1‑93 (MQ=255)

                                            tATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTc                                             <  2:1346364/93‑1 (MQ=255)

                                               ggTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGc                                          <  1:1994439/93‑1 (MQ=255)

                                                gTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCa                                         >  1:7278130/1‑93 (MQ=255)

                                                 tAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCaa                                        >  1:7093883/1‑93 (MQ=255)

                                                 tAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCaa                                        >  2:7211226/1‑93 (MQ=255)

                                                  aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAag                                       >  2:6844805/1‑93 (MQ=255)

                                                   aTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                      >  1:4226843/1‑93 (MQ=255)

                                                    tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                     <  1:1832609/93‑1 (MQ=255)

                                                    tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                     >  1:6663877/1‑93 (MQ=255)

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

                                                     aCCTTCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgc                                    <  1:5715323/93‑1 (MQ=255)

                                                      ccTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcg                                   >  2:7619832/1‑93 (MQ=255)

                                                        tGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                 <  2:2937727/93‑1 (MQ=255)

                                                        tGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                 >  1:4008039/1‑93 (MQ=255)

                                                        tGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCt                                 >  1:4449627/1‑93 (MQ=255)

                                                           ttCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTggg                              >  2:6025685/1‑93 (MQ=255)

                                                            tCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGt                             >  1:5124467/1‑93 (MQ=255)

                                                             cATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGtt                            <  1:1008731/93‑1 (MQ=255)

                                                               tttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACCGATCTAACGGACGACAGCGATAGTGTTTTTCTGCACGAGCGCTGGGTTtg                          >  2:747473/1‑93 (MQ=255)

                                                               tttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtg                          <  2:5110438/93‑1 (MQ=255)

                                                               tttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtg                          <  1:7446531/93‑1 (MQ=255)

                                                               tttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtg                          >  1:6924307/1‑93 (MQ=255)

                                                                ttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                         <  1:2533870/93‑1 (MQ=255)

                                                                   aGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGGTTGTGTc                      >  1:2813453/1‑93 (MQ=255)

                                                                    gATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTcc                     >  1:1326724/1‑93 (MQ=255)

                                                                    gATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTcc                     >  1:7161417/1‑93 (MQ=255)

                                                                     aTCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTccc                    >  1:6352094/1‑93 (MQ=255)

                                                                     aTCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTccc                    >  1:1557093/1‑93 (MQ=255)

                                                                      tCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCa                   <  2:2399027/93‑1 (MQ=255)

                                                                      tCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCa                   >  1:5639558/1‑93 (MQ=255)

                                                                       ccATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGGCGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                  <  1:1683476/93‑1 (MQ=255)

                                                                       ccATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                  >  2:3307400/1‑93 (MQ=255)

                                                                       ccATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                  >  2:3698257/1‑93 (MQ=255)

                                                                        cATGTCTCAATTTTCTCGCGCCATTTTCGTAACTGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                 >  1:5569259/1‑93 (MQ=255)

                                                                        cATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATa                 >  1:4191867/1‑93 (MQ=255)

                                                                         aTGTTTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                >  1:5250910/1‑93 (MQ=255)

                                                                          tGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCTTAAt               >  1:1112234/1‑93 (MQ=255)

                                                                          tGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAt               >  1:5470187/1‑93 (MQ=255)

                                                                          tGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAt               <  2:6233025/93‑1 (MQ=255)

                                                                          tGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCACATAAt               >  1:6779399/1‑93 (MQ=255)

                                                                           gTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAtt              >  1:6533637/1‑93 (MQ=255)

                                                                            tctcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGGTTGTGTCCCATAAttt             >  1:847674/1‑93 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           <  2:614888/93‑1 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           <  1:2478857/93‑1 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           >  2:7898485/1‑93 (MQ=255)

                                                                                aaTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGc         <  1:7781196/93‑1 (MQ=255)

                                                                                aaTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGc         <  2:4120938/93‑1 (MQ=255)

                                                                                 aTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGca        >  1:2169376/1‑93 (MQ=255)

                                                                                  ttttCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaa       >  2:3969033/1‑93 (MQ=255)

                                                                                     tctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaacaa    <  2:6419041/93‑1 (MQ=255)

                                                                                     tctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaacaa    <  2:6638847/93‑1 (MQ=255)

                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  <  2:3592704/93‑1 (MQ=255)

                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  >  2:2543525/1‑93 (MQ=255)

TATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATA  >  NC_000913/294900‑295079

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