BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,790,609	G→A	D239D (GAC→GAT)	htrL ←	HtrL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,790,609	0	G	A	99.5%	821.8 / ‑5.5	218	D239D (GAC→GAT)	htrL	HtrL
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (103/114); total (103/115)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 220 total aligned reads displayed.

GCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAAC  >  NC_000913/3790519‑3790698

gcGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATc                                                                                         <  1:5439383/93‑1 (MQ=255)

gcGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATc                                                                                         <  2:1690671/93‑1 (MQ=255)

  gAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCAt                                                                                       >  1:5363679/1‑93 (MQ=255)

  gAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCAt                                                                                       >  2:5810913/1‑93 (MQ=255)

  gAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCAt                                                                                       >  2:5044985/1‑93 (MQ=255)

  gAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCAt                                                                                       <  2:2351233/93‑1 (MQ=255)

    aaCAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCg                                                                                     >  2:3191793/1‑93 (MQ=255)

      cAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTc                                                                                   >  2:3878324/1‑93 (MQ=255)

          tCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCcaca                                                                               >  1:1960513/1‑93 (MQ=255)

           cGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACaa                                                                              <  2:3738157/93‑1 (MQ=255)

           cGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACaa                                                                              <  2:1573861/93‑1 (MQ=255)

            gAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAAt                                                                             <  2:4459995/93‑1 (MQ=255)

             aaCCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATg                                                                            >  2:4236231/1‑93 (MQ=255)

              aCCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGt                                                                           >  1:4388941/1‑93 (MQ=255)

              aCCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGt                                                                           <  2:106268/93‑1 (MQ=255)

              aCCATTTTCCTCCCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGt                                                                           <  2:5131200/93‑1 (MQ=255)

               ccATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGtt                                                                          <  1:2809376/93‑1 (MQ=255)

               ccATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGtt                                                                          <  1:1909559/93‑1 (MQ=255)

                cATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGtta                                                                         <  1:234335/93‑1 (MQ=255)

                 aTTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGttat                                                                        <  1:5976561/93‑1 (MQ=255)

                 aTTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGttat                                                                        <  2:2636409/93‑1 (MQ=255)

                 aTTTTCCGCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGttat                                                                        <  1:2158721/93‑1 (MQ=255)

                  ttttCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGttatt                                                                       <  1:5715354/93‑1 (MQ=255)

                   tttCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGttatta                                                                      <  2:5631054/93‑1 (MQ=255)

                   tttCCTCCCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGGTTGGAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGttatta                                                                      >  1:3334707/1‑93 (MQ=255)

                    ttCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTtattat                                                                     <  1:1839187/93‑1 (MQ=255)

                     tCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTattatt                                                                    >  1:6143422/1‑93 (MQ=255)

                     tCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTattatt                                                                    >  2:4927419/1‑93 (MQ=255)

                      ccTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTAttt                                                                   <  1:5429200/93‑1 (MQ=255)

                       ctctCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTa                                                                  <  2:2759154/93‑1 (MQ=255)

                        tctcCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAg                                                                 >  1:1141547/1‑93 (MQ=255)

                        tctcCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAg                                                                 <  2:2538325/93‑1 (MQ=255)

                        tctcCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAg                                                                 <  2:6134449/93‑1 (MQ=255)

                        tctcCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAg                                                                 >  2:5546779/1‑93 (MQ=255)

                          tcCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAgtg                                                               <  1:4861969/93‑1 (MQ=255)

                            cccAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGtt                                                             >  2:5655847/1‑93 (MQ=255)

                             ccAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTa                                                            >  2:929495/1‑93 (MQ=255)

                              cAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAt                                                           >  2:470247/1‑93 (MQ=255)

                              cAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAt                                                           >  1:4698903/1‑93 (MQ=255)

                               aGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAtt                                                          >  1:4966780/1‑93 (MQ=255)

                                gATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAttt                                                         >  1:151521/1‑93 (MQ=255)

                                gATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAttt                                                         >  2:4035162/1‑93 (MQ=255)

                                  tAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAttttt                                                       <  2:2420732/93‑1 (MQ=255)

                                   aGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAtttttt                                                      <  2:1185866/93‑1 (MQ=255)

                                   aGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTAtttttt                                                      <  1:5186602/93‑1 (MQ=255)

                                    gTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTg                                                     >  2:2197263/1‑93 (MQ=255)

                                     tttAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGa                                                    >  2:4798528/1‑93 (MQ=255)

                                     tttAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGa                                                    >  2:2786881/1‑93 (MQ=255)

                                       tAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACt                                                  >  1:3548774/1‑93 (MQ=255)

                                       tAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACt                                                  >  1:775953/1‑93 (MQ=255)

                                        aGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACtt                                                 <  2:3007097/93‑1 (MQ=255)

                                         ggTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACttt                                                >  2:5739263/1‑93 (MQ=255)

                                          gTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTc                                               <  2:3695694/93‑1 (MQ=255)

                                           ttGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTcc                                              <  1:310770/93‑1 (MQ=255)

                                           ttGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTcc                                              >  2:344316/1‑93 (MQ=255)

                                            tGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCa                                             <  2:4578420/93‑1 (MQ=255)

                                            tGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCa                                             <  2:3035023/93‑1 (MQ=255)

                                            tGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCa                                             <  1:5624544/93‑1 (MQ=255)

                                                gagaTCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAAc                                         >  1:3006368/1‑93 (MQ=255)

                                                  gaTCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTa                                       >  1:2469834/1‑93 (MQ=255)

                                                  gaTCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTa                                       >  2:3520110/1‑93 (MQ=255)

                                                   aTCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTaa                                      <  1:2254995/93‑1 (MQ=255)

                                                   aTCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTaa                                      <  1:1248079/93‑1 (MQ=255)

                                                    tCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAAt                                     <  1:2073341/93‑1 (MQ=255)

                                                    tCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAAt                                     >  2:4766782/1‑93 (MQ=255)

                                                      ggggCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAAttt                                   >  2:2888291/1‑93 (MQ=255)

                                                      ggggCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAAttt                                   >  1:3706896/1‑93 (MQ=255)

                                                      ggggCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAAttt                                   >  2:1524548/1‑93 (MQ=255)

                                                       gggCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTta                                  <  2:2066895/93‑1 (MQ=255)

                                                        ggCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTtat                                 <  2:2288830/93‑1 (MQ=255)

                                                         gCGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTtata                                <  2:3794302/93‑1 (MQ=255)

                                                          cGTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAg                               >  1:1690060/1‑93 (MQ=255)

                                                           gTTTGTAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGa                              >  2:3670979/1‑93 (MQ=255)

                                                                tAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCt                         <  2:6008546/93‑1 (MQ=255)

                                                                tAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCt                         <  1:5044350/93‑1 (MQ=255)

                                                                tAATAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCt                         <  2:765419/93‑1 (MQ=255)

                                                                 aaTAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCtt                        <  1:425318/93‑1 (MQ=255)

                                                                 aaTAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCtt                        <  1:3165357/93‑1 (MQ=255)

                                                                  aTAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCttt                       <  2:5052104/93‑1 (MQ=255)

                                                                   tAGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTc                      <  1:5728990/93‑1 (MQ=255)

                                                                    aGCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTcc                     <  1:2562242/93‑1 (MQ=255)

                                                                     gCACATGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCa                    <  1:1234652/93‑1 (MQ=255)

                                                                      cacaTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCAt                   <  2:1865300/93‑1 (MQ=255)

                                                                        caTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCAttt                 >  1:2686505/1‑93 (MQ=255)

                                                                         aTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTa                <  2:4694175/93‑1 (MQ=255)

                                                                         aTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTa                <  2:638964/93‑1 (MQ=255)

                                                                         aTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTa                <  1:2851725/93‑1 (MQ=255)

                                                                         aTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTa                <  1:376818/93‑1 (MQ=255)

                                                                         aTGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTa                >  2:5752576/1‑93 (MQ=255)

                                                                          tGACAAAAATGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTAt               >  2:4454042/1‑93 (MQ=255)

                                                                               aaaaTGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCtgt          <  1:942864/93‑1 (MQ=255)

                                                                               aaaaTGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCtgt          >  1:4190675/1‑93 (MQ=255)

                                                                               aaaaTGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCtgt          >  1:4041715/1‑93 (MQ=255)

                                                                                aaaTGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCtgtg         >  1:1723147/1‑93 (MQ=255)

                                                                                   tGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAAc      <  1:4178585/93‑1 (MQ=255)

                                                                                   tGCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAAc      >  2:1309330/1‑93 (MQ=255)

                                                                                    gCCCTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAAca     >  1:4457944/1‑92 (MQ=255)

                                                                                      ccTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACcaa   >  1:6016444/1‑93 (MQ=255)

                                                                                      ccTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACcaa   <  2:3547656/93‑1 (MQ=255)

                                                                                       cTGATCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACcaac  <  1:3845386/93‑1 (MQ=255)

GCGAAACAAATCGAACCATTTTCCTCTCCCCAGATAGTTTAGGTTGAAGAGATCGGGGCGTTTGTAATAGCACATGACAAAAATGCCCTGGTCATCGTCCACAATGTTATTATTTAGTGTTATTTTTTGACTTTCCAGAACTAATTTATAGAATTCTTTCCATTTATGCTGTGAACCAAC  >  NC_000913/3790519‑3790698

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