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	98.6%	445.1 / ‑2.7	143	R86R (CGT→CGC)	b1396	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (77/63); total (78/65)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 150 total aligned reads displayed.

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        >  1:3684301/1‑93 (MQ=255)

 ttttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCa                                                                                        <  2:2012709/93‑1 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAg                                                                                       <  1:1285714/93‑1 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                                                                       <  2:1366696/93‑1 (MQ=255)

  tttCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAg                                                                                       <  1:2678950/93‑1 (MQ=255)

    tCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGc                                                                                     <  1:3998978/93‑1 (MQ=255)

      aCTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCtt                                                                                   <  1:2868840/93‑1 (MQ=255)

      aCTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCtt                                                                                   <  1:894201/93‑1 (MQ=255)

       cTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCttt                                                                                  <  1:3360556/93‑1 (MQ=255)

        tGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTg                                                                                 <  2:3537676/93‑1 (MQ=255)

        tGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTg                                                                                 >  1:3570500/1‑93 (MQ=255)

        tGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTg                                                                                 >  1:2886390/1‑93 (MQ=255)

            tGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCgg                                                                             >  2:16385/1‑93 (MQ=255)

               tACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgaga                                                                          <  1:10285/93‑1 (MQ=255)

               tACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGgaga                                                                          >  1:3233918/1‑93 (MQ=255)

                  tGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGacac                                                                       <  1:2058426/93‑1 (MQ=255)

                  tGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTGTGCGCCCAGGGTTTGCCGGAGacac                                                                       >  1:1489439/1‑93 (MQ=255)

                   gCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACAcc                                                                      <  1:4344382/93‑1 (MQ=255)

                    ccTTTGCCTACGCCTGCGATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCt                                                                     <  1:4215761/93‑1 (MQ=255)

                      tttGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTa                                                                   >  2:4121602/1‑93 (MQ=255)

                        tGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAAc                                                                 <  2:2868618/93‑1 (MQ=255)

                         gCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACt                                                                >  2:435691/1‑93 (MQ=255)

                          ccTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACCCCTTAACTg                                                               >  2:3486491/1‑93 (MQ=255)

                           cTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGc                                                              >  2:1001070/1‑93 (MQ=255)

                           cTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGc                                                              >  2:3523485/1‑93 (MQ=255)

                           cTACGCCGGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGc                                                              <  2:4112832/93‑1 (MQ=255)

                              cGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  2:1204432/93‑1 (MQ=255)

                              cGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           <  2:1641629/93‑1 (MQ=255)

                              cGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTAc                                                           >  1:3619454/1‑93 (MQ=255)

                               gCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACt                                                          >  2:3503653/1‑93 (MQ=255)

                               gCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACt                                                          <  2:325185/93‑1 (MQ=255)

                                ccTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTg                                                         <  1:2369100/93‑1 (MQ=255)

                                 cTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        >  1:1732263/1‑93 (MQ=255)

                                 cTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgc                                                        <  2:901962/93‑1 (MQ=255)

                                   gCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTgcgc                                                      >  1:3598666/1‑93 (MQ=255)

                                    cAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCa                                                     <  2:4051552/93‑1 (MQ=255)

                                      aTAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAgg                                                   <  2:2979382/93‑1 (MQ=255)

                                       tAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGt                                                  <  2:224585/93‑1 (MQ=255)

                                        aGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTa                                                 >  1:4256007/1‑93 (MQ=255)

                                          ccAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACg                                               <  1:61344/93‑1 (MQ=255)

                                            aGGGGCTGGCAGCCGTCGCTTCTGCCTGCGCGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtc                                             >  2:1972663/1‑93 (MQ=255)

                                            aGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGtc                                             <  1:2833409/93‑1 (MQ=255)

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTNANTGCTACTGCGCAGNTACGncn                                            >  1:3435257/1‑92 (MQ=255)

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

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

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

                                             ggggCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGGACGtca                                            >  1:2559041/1‑93 (MQ=255)

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

                                                gCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCCCAGGTACGtcatca                                         <  1:1938451/93‑1 (MQ=255)

                                                   ggCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAggg                                      >  1:240174/1‑93 (MQ=255)

                                                   ggCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAggg                                      >  1:2798801/1‑93 (MQ=255)

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

                                                      aGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaa                                   <  1:4350540/93‑1 (MQ=255)

                                                         cGTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagca                                >  2:3691019/1‑93 (MQ=255)

                                                          gTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaa                               >  2:4519685/1‑93 (MQ=255)

                                                          gTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaa                               >  2:616607/1‑93 (MQ=255)

                                                          gTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaa                               >  1:1108476/1‑93 (MQ=255)

                                                          gTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaa                               >  1:408888/1‑93 (MQ=255)

                                                          gTCGCTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGgcaagcaa                               >  1:651183/1‑93 (MQ=255)

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

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

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

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

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

                                                              cTTCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCg                           <  1:2537274/93‑1 (MQ=255)

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

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

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

                                                               ttCTGCCTGCACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCgg                          <  1:3622174/93‑1 (MQ=255)

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

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

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

                                                                  tgcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGgtgt                       >  1:2081409/1‑93 (MQ=255)

                                                                   tccggcACGATTCATTTTTTGCGCCCAGGCTTTTCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGtt                      <  2:1837177/89‑1 (MQ=255)

                                                                   gcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGtt                      >  2:4363719/1‑93 (MQ=255)

                                                                   gcctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGtt                      <  1:4710972/93‑1 (MQ=255)

                                                                    cctgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGttt                     <  2:3158237/93‑1 (MQ=255)

                                                                      tgcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTAc                   <  2:1440406/93‑1 (MQ=255)

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

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

                                                                       gcACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGCCACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACg                  >  2:3169558/1‑93 (MQ=255)

                                                                        cACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGa                 >  2:103305/1‑93 (MQ=255)

                                                                        cACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGa                 <  1:3690995/93‑1 (MQ=255)

                                                                        cACGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGa                 >  2:4288659/1‑93 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               >  1:767276/1‑93 (MQ=255)

                                                                          cGATTGATTTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACGGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACa               >  2:4210940/1‑93 (MQ=255)

                                                                           gattgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACAt              >  2:29299/1‑93 (MQ=255)

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

                                                                            attgattTTTTGCG‑CCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGGCAACAGGGCAAGGCAACCGGTGTTTACGACATCg            >  2:4133653/1‑93 (MQ=255)

                                                                             ttgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCg            >  1:523443/1‑93 (MQ=255)

                                                                             ttgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCg            >  1:4087494/1‑93 (MQ=255)

                                                                             ttgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCg            <  2:2708674/93‑1 (MQ=255)

                                                                             ttgattTTTTGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCg            >  2:170423/1‑93 (MQ=255)

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

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

                                                                                    tttGCGCCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGt     <  2:591063/93‑1 (MQ=255)

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

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

                                                                                     ttGCGCCCACGCGTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGtt    <  1:163794/93‑1 (MQ=255)

                                                                                       gcgcCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTaa  <  2:3763214/93‑1 (MQ=255)

ATTTTCACTGGCTGATACTGCCTTTGCCTACGCCTGCAATAGCCAGGGGCTGGCAGCCGTCGCTTCTGCCTGCACGATTGATTTTTTGCGTCCAGGCTTTGCCGGAGACACCTTAACTGCTACTGCGCAGGTACGTCATCAGGGCAAGCAAACCGGTGTTTACGACATCGAAATTGTTAA  >  NC_000913/1458662‑1458841

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