BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	526,302	(G)_6→7	coding (3818/4281 nt)	rhsD →	RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	526,296	1	.	G	99.2%	404.5 / ‑3.6	128	coding (3812/4281 nt)	rhsD	RhsD
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (73/54); total (74/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

Only 100 of 148 total aligned reads displayed.

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑GGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGA  >  NC_000913/526204‑526383

ggaTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATTGGATTACTGCAGACTTGGGATGATGACAGGTCTGGAGCATGTAc                                                                                          >  2:3308267/3‑93 (MQ=255)

gAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAc                                                                                          >  1:76527/1‑93 (MQ=255)

 aaTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑g                                                                                        <  1:3546927/93‑1 (MQ=255)

  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGCTGATGCCAGCTCTGGAGCATGTAC‑gg                                                                                       >  2:4087816/1‑93 (MQ=255)

  aTTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑gg                                                                                       <  2:2140221/93‑1 (MQ=255)

   tttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑ggg                                                                                      <  2:2748300/93‑1 (MQ=255)

    ttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑gggg                                                                                     <  2:1947768/93‑1 (MQ=255)

    ttATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑gggg                                                                                     <  2:1169331/93‑1 (MQ=255)

     tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑ggggg                                                                                    >  1:4462941/1‑93 (MQ=255)

     tataTCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑ggggg                                                                                    >  2:4061757/1‑93 (MQ=255)

      atgtCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑gggggg                                                                                   <  2:386950/90‑1 (MQ=255)

       tatCAGTCTCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGACAGATCTGGAGCATGTAC‑gggggcg                                                                                  >  2:2928739/1‑91 (MQ=255)

       tatCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑ggggggg                                                                                  >  2:322380/1‑92 (MQ=255)

        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑ggggggga                                                                                 >  2:91858/1‑91 (MQ=255)

        atCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑ggggggga                                                                                 >  2:1675893/1‑91 (MQ=255)

            gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGttt                                                                              >  1:2134523/1‑93 (MQ=255)

            gTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGttt                                                                              >  1:618379/1‑93 (MQ=255)

              aTCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTtgt                                                                            >  2:4436917/1‑93 (MQ=255)

              aTCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTtgt                                                                            <  1:4567331/93‑1 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTtgtg                                                                           >  1:3467136/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTtgtg                                                                           >  2:2281146/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTtgtg                                                                           >  1:2050512/1‑93 (MQ=255)

               tCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGGTtgtg                                                                           >  1:3347248/1‑93 (MQ=255)

                ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTgtgg                                                                          >  2:1350197/1‑93 (MQ=255)

                ccTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTgtgg                                                                          >  2:3269099/1‑93 (MQ=255)

                 cTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTgtggt                                                                         <  1:1448994/93‑1 (MQ=255)

                  tttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTgtggtg                                                                        >  1:556614/1‑93 (MQ=255)

                   ttAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGgtgt                                                                       <  2:785377/93‑1 (MQ=255)

                    tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGtt                                                                      >  2:2388779/1‑93 (MQ=255)

                    tAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGtt                                                                      >  1:4641342/1‑93 (MQ=255)

                     aaaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGttt                                                                     >  1:1018215/1‑92 (MQ=255)

                      aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCt                                                                    >  2:3995482/1‑93 (MQ=255)

                      aaTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCt                                                                    >  2:1222502/1‑93 (MQ=255)

                       aTCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCtt                                                                   <  2:2715111/93‑1 (MQ=255)

                        tCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCttt                                                                  <  2:2772475/93‑1 (MQ=255)

                         ccACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGGGTTTGTGGTGTTCTTTc                                                                 >  2:1224039/1‑93 (MQ=255)

                         ccACGACAACAAATTGACCCTATGGGATTAGTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTc                                                                 <  2:1168928/93‑1 (MQ=255)

                          cACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCa                                                                >  1:1919145/1‑93 (MQ=255)

                             tACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGGGTTTGTGGTGTTCTTTCACGt                                                             >  2:3854815/1‑93 (MQ=255)

                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCCGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTataa                                                         >  1:848887/1‑93 (MQ=255)

                                 acaaATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTataa                                                         <  1:3059877/93‑1 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTataata                                                       >  1:1989164/1‑93 (MQ=255)

                                   aaaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTataata                                                       <  1:2355273/93‑1 (MQ=255)

                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGTATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAg                                                      <  1:2518764/93‑1 (MQ=255)

                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAg                                                      >  1:3320532/1‑93 (MQ=255)

                                    aaTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAg                                                      <  2:2157505/93‑1 (MQ=255)

                                     aTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAgg                                                     >  2:596882/1‑93 (MQ=255)

                                     aTTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAgg                                                     <  1:3611688/93‑1 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                    >  2:4155343/1‑93 (MQ=255)

                                      ttGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGa                                                    <  2:3463996/93‑1 (MQ=255)

                                         aCCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCa                                                 >  1:3776745/1‑93 (MQ=255)

                                          cccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCaa                                                <  1:4545166/93‑1 (MQ=255)

                                           ccTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGGGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAg                                               >  2:1768446/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             <  2:1449885/93‑1 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  2:3627430/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  2:3531408/1‑93 (MQ=255)

                                             tATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGCTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTa                                             >  2:253326/1‑93 (MQ=255)

                                              aTGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaa                                            >  1:1100248/1‑93 (MQ=255)

                                               tGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTaaa                                           >  1:945236/1‑93 (MQ=255)

                                                  gATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAAttt                                        <  2:3721042/93‑1 (MQ=255)

                                                   attattGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                       <  2:2078278/93‑1 (MQ=255)

                                                   attattGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                       <  1:1446179/93‑1 (MQ=255)

                                                   attattGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTg                                       >  1:2211250/1‑93 (MQ=255)

                                                      attGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑GGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  1:1197436/1‑93 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  1:995739/1‑93 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  1:520141/1‑93 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   >  1:1486509/1‑93 (MQ=255)

                                                       ttGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAg                                   <  1:2112637/93‑1 (MQ=255)

                                                        tGCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGt                                  >  2:486809/1‑93 (MQ=255)

                                                         gCAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTa                                 <  2:3957631/93‑1 (MQ=255)

                                                          cAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                <  2:84187/93‑1 (MQ=255)

                                                          cAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                <  2:451052/93‑1 (MQ=255)

                                                          cAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                <  2:3043656/93‑1 (MQ=255)

                                                          cAGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTAc                                <  2:4211283/93‑1 (MQ=255)

                                                           aGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                               <  2:684599/93‑1 (MQ=255)

                                                           aGACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACt                               >  1:4514985/1‑93 (MQ=255)

                                                            gACTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTg                              >  1:2525160/1‑93 (MQ=255)

                                                              ctggggATGATGCCATATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                            <  2:1626961/90‑1 (MQ=255)

                                                              cTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                            >  1:3394005/1‑93 (MQ=255)

                                                              cTTGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCa                            >  2:1434402/1‑93 (MQ=255)

                                                               ttGGGATGATGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAg                           >  2:1766300/1‑93 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       <  2:688030/93‑1 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       >  1:517809/1‑93 (MQ=255)

                                                                   gatgatGCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgc                       <  1:1394775/93‑1 (MQ=255)

                                                                    atgatgCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgct                      >  2:1353794/1‑93 (MQ=255)

                                                                     tgatgCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGAtgctg                     >  1:3006699/1‑93 (MQ=255)

                                                                       atgCCAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCg                   >  2:228338/1‑93 (MQ=255)

                                                                          ccAGATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTa                <  1:576174/93‑1 (MQ=255)

                                                                             gATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt             >  2:394523/1‑93 (MQ=255)

                                                                             gATCTGGAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGAt             >  1:1698599/1‑93 (MQ=255)

                                                                                   gAGCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTg       <  1:1480460/93‑1 (MQ=255)

                                                                                     gCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaa     <  2:3523495/93‑1 (MQ=255)

                                                                                     gCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaa     <  1:1364889/93‑1 (MQ=255)

                                                                                     gCATGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaa     <  2:894035/93‑1 (MQ=255)

                                                                                        tGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  <  2:3926350/93‑1 (MQ=255)

                                                                                        tGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  <  1:4091671/93‑1 (MQ=255)

                                                                                        tGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  >  2:3353839/1‑93 (MQ=255)

                                                                                        tGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  >  2:1413937/1‑93 (MQ=255)

                                                                                        tGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  >  2:1141230/1‑93 (MQ=255)

                                                                                        tGTACGGGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTgaaaga  <  2:2625125/93‑1 (MQ=255)

GAATTTATATCAGTATCCTTTAAATCCACTACAACAAATTGACCCTATGGGATTATTGCAGACTTGGGATGATGCCAGATCTGGAGCATGTAC‑GGGGGGAGTTTGTGGTGTTCTTTCACGTATAATAGGACCAAGTAAATTTGATAGTACTGCAGATGCTGCGTTAGATGCTTTGAAAGA  >  NC_000913/526204‑526383

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