BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,458,752	T→C	R86R (CGT→CGC)	b1396 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,458,752	0	T	C	99.6%	695.8 / ‑5.2	224	R86R (CGT→CGC)	b1396	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (98/125); total (98/126)
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.62e-01

Only 100 of 234 total aligned reads displayed.

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

aTTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGccc                                                                                         <  2:1537100/93‑1 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                                                                       <  2:4515687/93‑1 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                                                                       <  1:2308087/93‑1 (MQ=255)

    tCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGc                                                                                     <  1:4199323/93‑1 (MQ=255)

      aCTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCtt                                                                                   >  2:6552605/1‑93 (MQ=255)

        tGGCTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTg                                                                                 >  1:3277636/1‑93 (MQ=255)

           cTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCg                                                                              >  2:6351190/1‑93 (MQ=255)

           cTGATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCg                                                                              >  2:3214248/1‑93 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCgca                                                                            >  1:7216017/1‑91 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            >  2:2058466/1‑93 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            >  2:1965840/1‑93 (MQ=255)

             gATACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGga                                                                            <  2:4807135/93‑1 (MQ=255)

              aTACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgag                                                                           >  2:5425525/1‑93 (MQ=255)

               tACTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgaga                                                                          <  2:6471955/93‑1 (MQ=255)

                 cTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGaca                                                                        >  2:2455120/1‑93 (MQ=255)

                 cTGCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGaca                                                                        >  2:5808058/1‑93 (MQ=255)

                   gCCTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                                                                      <  1:1815978/93‑1 (MQ=255)

                   gCCTTTGCCTACGCCTGCAATAGTCAGGGACTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                                                                      <  2:906912/93‑1 (MQ=255)

                    ccTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     <  2:6115932/93‑1 (MQ=255)

                     cTTTGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCtt                                                                    >  2:6263047/1‑93 (MQ=255)

                      tttGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGGCGGGGACACCTTa                                                                   >  2:5910009/1‑93 (MQ=255)

                       ttGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTaa                                                                  <  1:4035870/93‑1 (MQ=255)

                       ttGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTaa                                                                  <  2:3816306/93‑1 (MQ=255)

                        tGCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAAc                                                                 <  1:750241/93‑1 (MQ=255)

                         gCCTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACt                                                                <  1:207039/93‑1 (MQ=255)

                         gCCTACGCCTGCAATAGTCAGGGGCAGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACt                                                                >  1:3364735/1‑93 (MQ=255)

                          ccTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTg                                                               >  1:4225941/1‑93 (MQ=255)

                          ccTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTg                                                               <  1:6932584/93‑1 (MQ=255)

                          ccTACGCCTGCAATAGTCAGGCGCTGGCAGCCGTCGCTTCTGCCTGCAAGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTg                                                               <  2:675216/93‑1 (MQ=255)

                           cTACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGc                                                              <  2:4367306/93‑1 (MQ=255)

                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  2:5300991/1‑93 (MQ=255)

                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  1:60688/1‑93 (MQ=255)

                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  1:5552275/1‑93 (MQ=255)

                            tACGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCt                                                             >  2:4964547/1‑93 (MQ=255)

                             aCGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTa                                                            >  2:84276/1‑93 (MQ=255)

                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  2:2675178/93‑1 (MQ=255)

                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  1:7853307/93‑1 (MQ=255)

                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  2:4233561/93‑1 (MQ=255)

                              cGCCTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  1:1217122/93‑1 (MQ=255)

                                 cTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        <  1:3876418/93‑1 (MQ=255)

                                 cTGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        >  2:3906807/1‑93 (MQ=255)

                                  tGCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcg                                                       >  1:4996547/1‑93 (MQ=255)

                                   gCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcgc                                                      <  1:7467235/93‑1 (MQ=255)

                                   gCAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcgc                                                      >  1:282870/1‑93 (MQ=255)

                                    cAATAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCa                                                     <  2:3602808/93‑1 (MQ=255)

                                     aaTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTAAACTGCTACTGCGCAg                                                    <  1:461710/93‑1 (MQ=255)

                                      aTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg                                                   <  1:1316804/93‑1 (MQ=255)

                                      aTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg                                                   <  1:3053434/93‑1 (MQ=255)

                                      aTAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg                                                   >  2:3330235/1‑93 (MQ=255)

                                       tAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGt                                                  >  1:2671667/1‑93 (MQ=255)

                                       tAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGt                                                  <  1:5322872/93‑1 (MQ=255)

                                       tAGTCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGt                                                  >  1:5272937/1‑93 (MQ=255)

                                        agtcAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTa                                                 <  1:2716701/90‑1 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtca                                            >  2:1526542/1‑93 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtca                                            <  2:6996497/93‑1 (MQ=255)

                                               ggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatc                                          >  2:4674004/1‑93 (MQ=255)

                                                gCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTGCTGCGCAGGTACGtcatca                                         <  2:2191434/93‑1 (MQ=255)

                                                gCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtcatca                                         <  2:259178/93‑1 (MQ=255)

                                                    gCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgc                                     >  1:3539625/1‑93 (MQ=255)

                                                    gCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgc                                     <  2:4749677/93‑1 (MQ=255)

                                                      aGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaa                                   >  1:6127284/1‑93 (MQ=255)

                                                       gCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaag                                  >  2:7170581/1‑93 (MQ=255)

                                                        ccGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagc                                 <  2:6018220/93‑1 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              >  1:5308408/1‑93 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              >  2:7782673/1‑93 (MQ=255)

                                                           tCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCaaa                              <  2:2781872/93‑1 (MQ=255)

                                                            cGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAc                             <  1:6354399/93‑1 (MQ=255)

                                                             gCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAcc                            >  2:5260853/1‑93 (MQ=255)

                                                             gCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAAcc                            >  1:6762583/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  2:1050331/1‑93 (MQ=255)

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          >  1:6470033/1‑93 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgt                         >  1:1459921/1‑93 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgt                         <  2:5603836/93‑1 (MQ=255)

                                                                tCTGCCTGCACGATTGATTTTTTGCGCCCAGGCGTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgt                         <  1:7650365/93‑1 (MQ=255)

                                                                 ctgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtg                        <  1:543663/93‑1 (MQ=255)

                                                                 ctgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtg                        >  2:1167816/1‑93 (MQ=255)

                                                                 ctgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtg                        >  2:7845808/1‑93 (MQ=255)

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  1:1060159/93‑1 (MQ=255)

                                                                  ctccttcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       <  1:1942759/91‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  1:7566464/93‑1 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    >  1:489715/1‑93 (MQ=255)

                                                                     ctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTa                    <  2:3420980/93‑1 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   >  1:4497491/1‑93 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   >  2:4223171/1‑93 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGGACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   >  1:5390256/1‑93 (MQ=255)

                                                                       gcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACg                  <  1:5719053/93‑1 (MQ=255)

                                                                         aCGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGAc                >  2:7653282/1‑93 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               <  2:2977234/93‑1 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               <  2:2161713/93‑1 (MQ=255)

                                                                           gattgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACAt              <  1:6965387/93‑1 (MQ=255)

                                                                            attgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATc             <  1:5072392/93‑1 (MQ=255)

                                                                            attgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATc             <  2:4459071/93‑1 (MQ=255)

                                                                               gattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaa          >  1:4793394/1‑93 (MQ=255)

                                                                                attTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGgaa         <  2:335022/93‑1 (MQ=255)

                                                                                attTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGaaa         <  1:6580543/93‑1 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        <  2:3436958/93‑1 (MQ=255)

                                                                                 ttttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAt        <  2:2137323/93‑1 (MQ=255)

                                                                                  tttttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAAtt       <  1:3114541/93‑1 (MQ=255)

                                                                                     ttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGtt    <  1:170631/93‑1 (MQ=255)

                                                                                     ttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGtt    >  2:4299308/1‑93 (MQ=255)

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

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