BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,233,455	T→C	D21G (GAC→GGC)	b3084 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,233,455	0	T	C	99.4%	544.1 / ‑4.9	171	D21G (GAC→GGC)	b3084	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (91/78); total (92/79)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 179 total aligned reads displayed.

AGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGAC  >  NC_000913/3233365‑3233544

aGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTGCAACGTCTGTAGCGGGTTAACGCCa                                                                                         >  1:2361418/1‑93 (MQ=255)

  gATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATc                                                                                       <  2:1176762/93‑1 (MQ=255)

   aTCACCACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATcc                                                                                      <  2:3548101/93‑1 (MQ=255)

   aTCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATcc                                                                                      >  2:3770994/1‑93 (MQ=255)

     cAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGt                                                                                    >  2:2140759/1‑93 (MQ=255)

      aaCACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTc                                                                                   <  2:4502402/93‑1 (MQ=255)

      aaCACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTc                                                                                   >  1:180933/1‑93 (MQ=255)

       acacCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTcc                                                                                  >  1:2016170/1‑92 (MQ=255)

       acacCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCg                                                                                  >  2:1267809/1‑93 (MQ=255)

       acacCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCg                                                                                  >  1:2016171/1‑93 (MQ=255)

        cacCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGc                                                                                 >  1:4716822/1‑93 (MQ=255)

         acCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGcc                                                                                >  1:4397183/1‑93 (MQ=255)

            gggCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCgca                                                                             >  2:1519840/1‑91 (MQ=255)

            gggCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGa                                                                             <  1:315622/93‑1 (MQ=255)

               ccGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGaaaa                                                                          <  2:3194091/93‑1 (MQ=255)

                cgcgGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAAc                                                                         >  1:4354509/1‑93 (MQ=255)

                  cgGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGt                                                                       <  2:346880/93‑1 (MQ=255)

                   ggATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGtt                                                                      <  2:4041058/93‑1 (MQ=255)

                      tttCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTa                                                                   >  1:2184227/1‑93 (MQ=255)

                      tttCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTa                                                                   >  2:2049510/1‑93 (MQ=255)

                      tttCTGTGTCGGCCAACTGTTGCAGCAAATATTCAGCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTa                                                                   <  1:950801/93‑1 (MQ=255)

                       tgtggtgtCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAg                                                                  <  1:750459/89‑1 (MQ=255)

                        tCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAAGGCCATCCGTCGCCGGAAAACGCTGTAgt                                                                 >  2:784586/1‑93 (MQ=255)

                           gtgtCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTc                                                              >  1:3285029/1‑93 (MQ=255)

                             gtcgtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAg                                                            >  2:3827073/1‑93 (MQ=255)

                             gtcgtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAg                                                            <  2:2429820/93‑1 (MQ=255)

                             gtcgtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAg                                                            <  1:3490568/93‑1 (MQ=255)

                             gtcgtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCACACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGATGTGTccg                                                            >  1:1376093/1‑91 (MQ=255)

                               cgtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTg                                                          >  2:1231562/1‑93 (MQ=255)

                                gtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGa                                                         >  1:4296424/1‑93 (MQ=255)

                                gtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGa                                                         >  1:4099494/1‑93 (MQ=255)

                                gtcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGa                                                         >  2:3898656/1‑93 (MQ=255)

                                 tcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGaa                                                        >  2:3196346/1‑93 (MQ=255)

                                 tcCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGaa                                                        >  2:2457100/1‑93 (MQ=255)

                                     aCTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGa                                                    <  2:27958/93‑1 (MQ=255)

                                       tGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGaaa                                                  >  1:1569450/1‑93 (MQ=255)

                                       tGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGaaa                                                  >  1:4463571/1‑93 (MQ=255)

                                        tttGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAAc                                                 <  2:48845/92‑1 (MQ=255)

                                         ttGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAAcc                                                <  1:3231169/93‑1 (MQ=255)

                                         ttGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAAcc                                                <  2:4404891/93‑1 (MQ=255)

                                           gcagcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGt                                              <  1:1226812/93‑1 (MQ=255)

                                           gcagcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGt                                              >  2:1682196/1‑93 (MQ=255)

                                             agcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTg                                            <  2:667368/93‑1 (MQ=255)

                                             agcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTg                                            <  2:1676067/93‑1 (MQ=255)

                                             agcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTg                                            <  1:2358555/93‑1 (MQ=255)

                                              gcaAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTTAACGGAAACCGTTGt                                           >  2:3360213/1‑93 (MQ=255)

                                                aaaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCt                                         >  2:240907/1‑93 (MQ=255)

                                                 aaTATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTa                                        <  1:2815519/93‑1 (MQ=255)

                                                  atatTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTaa                                       >  1:3454894/1‑93 (MQ=255)

                                                   tatTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAg                                      >  1:192188/1‑93 (MQ=255)

                                                   tatTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAg                                      >  2:1746824/1‑93 (MQ=255)

                                                   tatTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAg                                      <  2:2776359/93‑1 (MQ=255)

                                                   aatTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAg                                      <  1:4134323/92‑1 (MQ=255)

                                                    atTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGt                                     <  2:930345/93‑1 (MQ=255)

                                                     ttCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTg                                    <  1:1020890/93‑1 (MQ=255)

                                                     ttCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTg                                    >  2:1051401/1‑93 (MQ=255)

                                                     ttCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTg                                    <  1:3748562/93‑1 (MQ=255)

                                                      tCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTgg                                   >  1:4205973/1‑93 (MQ=255)

                                                       cATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGc                                  <  1:3321142/93‑1 (MQ=255)

                                                        aTCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCt                                 >  1:3039228/1‑93 (MQ=255)

                                                        aTCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCt                                 >  1:1682777/1‑93 (MQ=255)

                                                        aTCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCt                                 >  1:1738500/1‑93 (MQ=255)

                                                          ccgccgCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTACGTGGCTCa                               >  1:3387749/1‑93 (MQ=255)

                                                          ccgccgCTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCa                               <  2:4682714/93‑1 (MQ=255)

                                                           cgccgcTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCAt                              <  1:1087895/93‑1 (MQ=255)

                                                            gccgcTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATc                             <  1:3904109/93‑1 (MQ=255)

                                                            gccgcTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATc                             <  2:2497360/93‑1 (MQ=255)

                                                            gccgcTTCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGCCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATc                             <  1:2609316/93‑1 (MQ=255)

                                                                 ttCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGccc                        >  2:630588/1‑93 (MQ=255)

                                                                 ttCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGccc                        >  2:284448/1‑93 (MQ=255)

                                                                 ttCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGccc                        <  2:1402317/93‑1 (MQ=255)

                                                                 ttCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGccc                        <  1:3861004/93‑1 (MQ=255)

                                                                 ttCCCACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGAAAACCGTTGTCTAAGTGGCTCATCGGccc                        >  1:539191/1‑93 (MQ=255)

                                                                   cccACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctc                      >  2:2613722/1‑93 (MQ=255)

                                                                   cccACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctc                      >  1:866935/1‑93 (MQ=255)

                                                                    ccACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctcc                     >  2:4180493/1‑93 (MQ=255)

                                                                     cACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctcct                    <  1:3794681/93‑1 (MQ=255)

                                                                     cACGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCcctcct                    <  1:1615099/93‑1 (MQ=255)

                                                                      aCGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTg                   >  1:368037/1‑93 (MQ=255)

                                                                      aCGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTg                   >  2:3693534/1‑93 (MQ=255)

                                                                       cGCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGa                  >  2:4211745/1‑93 (MQ=255)

                                                                        gCCTGTAGCGGGTTAACGCCATCCGTCGCCGGGAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGaa                 >  1:2106618/1‑93 (MQ=255)

                                                                        gCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGaa                 >  1:2532706/1‑93 (MQ=255)

                                                                        gCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGaa                 <  2:3281193/93‑1 (MQ=255)

                                                                        gCCTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGaa                 >  2:3156743/1‑93 (MQ=255)

                                                                         ccTGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAAt                <  2:4375551/93‑1 (MQ=255)

                                                                           tGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGa              <  1:4665298/93‑1 (MQ=255)

                                                                           tGTAGCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGa              >  2:2989787/1‑93 (MQ=255)

                                                                           tGTAGCGGGTTAACGCCATCCGTCGCAGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGa              >  2:614333/1‑93 (MQ=255)

                                                                               gCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATaaa          <  1:992056/93‑1 (MQ=255)

                                                                               gCGGGTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATaaa          <  1:2393347/93‑1 (MQ=255)

                                                                                 gggTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAAt        <  1:2958962/93‑1 (MQ=255)

                                                                                 gggTTAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTCGCTCATCGGCCCTCCTGAATGATAAAAt        >  1:2563885/1‑93 (MQ=255)

                                                                                    ttAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTc     <  1:1401044/93‑1 (MQ=255)

                                                                                    ttAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTc     <  1:1093857/93‑1 (MQ=255)

                                                                                    ttAACGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTc     >  1:1014783/1‑93 (MQ=255)

                                                                                      aaCGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGa   >  2:3456045/1‑93 (MQ=255)

                                                                                      aaCGCCATCCGGCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGa   <  1:4637816/93‑1 (MQ=255)

                                                                                       aCGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGAc  <  1:2029743/93‑1 (MQ=255)

                                                                                       aCGCCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGAc  <  1:1257322/93‑1 (MQ=255)

AGGATCAACACCGGGCCGCGGATTTCTGTGTCGTCCAACTGTTGCAGCAAATATTCATCCGCCGCTTCCCACGCCTGTAGCGGGTTAACGTCATCCGTCGCCGGAAAACGTTGTAGTGTCAGTGAACGGAAACCGTTGTCTAAGTGGCTCATCGGCCCTCCTGAATGATAAAATTTCGAC  >  NC_000913/3233365‑3233544

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