BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,991,301	C→T	P50L (CCG→CTG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,991,301	0	C	T	99.2%	459.3 / ‑5.1	133	P50L (CCG→CTG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (61/71); total (61/72)
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 133 total aligned reads displayed.

TTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGT  >  NC_000913/3991204‑3991395

ttGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTgg                                                                                              <  1:810318/100‑1 (MQ=255)

 tGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTggg                                                                                             <  1:2059881/100‑1 (MQ=255)

 tGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTggg                                                                                             <  2:2022010/100‑1 (MQ=255)

 tGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTggg                                                                                             >  2:2845190/1‑100 (MQ=255)

  gCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGt                                                                                            >  2:2776945/1‑100 (MQ=255)

     tgtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTAc                                                                                         >  1:874473/1‑100 (MQ=255)

      gtgGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTAcc                                                                                        <  1:2168557/100‑1 (MQ=255)

         gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTg                                                                                     >  1:821084/1‑100 (MQ=255)

         gATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTg                                                                                     >  1:1387410/1‑100 (MQ=255)

          aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGTCATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGa                                                                                    >  2:2083133/1‑100 (MQ=255)

          aTCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGa                                                                                    >  2:3159245/1‑100 (MQ=255)

           tCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGAt                                                                                   <  1:3168046/100‑1 (MQ=255)

            cgcgcgCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATg                                                                                  <  1:1758809/100‑1 (MQ=255)

            cgcgcgCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATg                                                                                  <  2:2306978/100‑1 (MQ=255)

             gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATgg                                                                                 <  1:152029/100‑1 (MQ=255)

             gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATgg                                                                                 <  1:493560/100‑1 (MQ=255)

             gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATgg                                                                                 <  1:1131620/100‑1 (MQ=255)

             gcgcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATgg                                                                                 <  1:2870331/100‑1 (MQ=255)

             gcgcgcTTGCTGCTATGGGGCCTGCATTACAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATgg                                                                                 <  1:144807/100‑1 (MQ=255)

               gcgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGCTGCACTATCACCATCCGCTAATGCTGGGCTACCTTGATggca                                                                               <  1:2207566/100‑3 (MQ=255)

               gcgcTTGCTGCTATGGGGCCTGCATACCAACAGGTCTACAGTCGACTGCCGACATTGTTGCAATATCACCATCCGCTAATGCTGGGTTACCTTGATgg                                                                                 <  2:220248/98‑1 (MQ=255)

                cgcTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTaa                                                                              >  1:985841/1‑100 (MQ=255)

                  cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAGCg                                                                            <  2:2974120/100‑1 (MQ=255)

                  cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                            >  2:1301755/1‑100 (MQ=255)

                  cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                            >  1:1083883/1‑100 (MQ=255)

                  cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                            >  1:215017/1‑100 (MQ=255)

                  cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                            >  1:127243/1‑100 (MQ=255)

                  cTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACg                                                                            >  2:2648709/1‑100 (MQ=255)

                    tgctgcTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATTCTGGGTTACCTTGATGGTAACGtt                                                                          >  2:2588728/1‑100 (MQ=255)

                      ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTcc                                                                        >  2:81204/1‑100 (MQ=255)

                      ctgctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTcc                                                                        >  2:1773254/1‑100 (MQ=255)

                        gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCa                                                                      <  1:2011238/100‑1 (MQ=255)

                        gctATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCa                                                                      >  2:2389985/1‑100 (MQ=255)

                          tATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCaaa                                                                    >  2:2556982/1‑100 (MQ=255)

                           aTGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAg                                                                   <  1:2792448/100‑1 (MQ=255)

                            tGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAgg                                                                  >  2:863109/1‑100 (MQ=255)

                                 ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAttt                                                             <  1:1789540/100‑1 (MQ=255)

                                 ccTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAttt                                                             >  2:701514/1‑100 (MQ=255)

                                  cTGCATTCCAACCGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                            <  1:2988127/100‑1 (MQ=255)

                                  cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                            <  1:2267839/100‑1 (MQ=255)

                                  cTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTg                                                            <  2:1297749/100‑1 (MQ=255)

                                    gCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGcc                                                          <  2:634597/100‑1 (MQ=255)

                                      aTTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCtt                                                        <  2:875487/100‑1 (MQ=255)

                                       ttCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCttt                                                       <  2:1881030/100‑1 (MQ=255)

                                        tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTa                                                      >  1:533560/1‑100 (MQ=255)

                                        tCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttg                                                      <  1:73024/100‑2 (MQ=255)

                                         ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTCTACCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTac                                                     >  2:1041304/1‑100 (MQ=255)

                                         ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGCTGGTAACGTTCCCAAAGGCATTTGCCTTTac                                                     >  1:203518/1‑100 (MQ=255)

                                         ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTac                                                     >  1:3150638/1‑100 (MQ=255)

                                         ccAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgg                                                     <  1:1643989/100‑3 (MQ=255)

                                          cAACAGGTCTACAGTCTACTGCCGACATTTTTGCACTATCATCATCCTCTAATTCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttggg                                                    <  1:1085040/100‑4 (MQ=255)

                                           aaCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTacac                                                   <  1:1509240/100‑1 (MQ=255)

                                            aCAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAATGCCTTTgggaacgttac                                                  <  1:912820/100‑11 (MQ=255)

                                              aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCTAAAGGCAAATGCCtttgggaacg                                                <  1:684149/100‑8 (MQ=255)

                                              aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgggaacg                                                <  1:6712/100‑8 (MQ=255)

                                              aGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCAAATGCCtttgggaacg                                                <  1:6864/100‑8 (MQ=255)

                                               ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCt                                               <  2:705055/100‑1 (MQ=255)

                                               ggTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCt                                               >  1:1224205/1‑100 (MQ=255)

                                                  cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                            >  2:1409127/1‑100 (MQ=255)

                                                  cTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgat                                            <  1:705948/100‑1 (MQ=255)

                                                    aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgatga                                          >  2:2166473/1‑100 (MQ=255)

                                                    aCAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCtgatga                                          >  2:1712116/1‑100 (MQ=255)

                                                     cAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaa                                         <  1:1310096/100‑1 (MQ=255)

                                                      aGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGaaa                                        <  1:768169/100‑1 (MQ=255)

                                                       gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                       >  1:1173250/1‑100 (MQ=255)

                                                       gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                       <  1:1036713/100‑1 (MQ=255)

                                                       gTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAAc                                       <  1:2552546/100‑1 (MQ=255)

                                                          tACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCa                                    >  1:345894/1‑100 (MQ=255)

                                                           aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCaa                                   >  1:1344635/1‑100 (MQ=255)

                                                           aCTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCaa                                   >  1:1976622/1‑100 (MQ=255)

                                                            cTGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAAc                                  <  1:627860/100‑1 (MQ=255)

                                                             tGCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACg                                 <  1:1849142/100‑1 (MQ=255)

                                                              gCCGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGc                                >  2:865414/1‑100 (MQ=255)

                                                                cGACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCa                              >  1:2826108/1‑100 (MQ=255)

                                                                 gACATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCAc                             >  1:970944/1‑100 (MQ=255)

                                                                   cATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTa                           >  2:820268/1‑100 (MQ=255)

                                                                   cATTGTTGCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTa                           >  2:821302/1‑100 (MQ=255)

                                                                         tgCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGaa                     >  2:2713586/1‑100 (MQ=255)

                                                                          gCACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAAc                    >  1:2233322/1‑100 (MQ=255)

                                                                           cACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                   <  2:1317919/100‑1 (MQ=255)

                                                                           cACTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACg                   <  2:3013418/100‑1 (MQ=255)

                                                                            aCTATCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGa                  >  2:2885947/1‑100 (MQ=255)

                                                                               aTCACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCt               >  1:2232505/1‑100 (MQ=255)

                                                                                ttaCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCtt              >  1:2208705/3‑100 (MQ=255)

                                                                                 cACCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg             <  1:2412788/100‑1 (MQ=255)

                                                                                 cACCATCCGCTAATGATGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTg             <  1:3208742/100‑1 (MQ=255)

                                                                                  aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa            >  1:2400068/1‑100 (MQ=255)

                                                                                  aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa            <  1:2417918/100‑1 (MQ=255)

                                                                                  aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa            >  2:1387459/1‑100 (MQ=255)

                                                                                  aCCATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGa            <  1:1252003/100‑1 (MQ=255)

                                                                                    cATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc          <  2:1374309/100‑1 (MQ=255)

                                                                                    cATCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAAc          <  1:987891/100‑1 (MQ=255)

                                                                                     aTCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt         <  1:2776980/100‑1 (MQ=255)

                                                                                     aTCCGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACt         <  1:925452/100‑1 (MQ=255)

                                                                                       ccGCTAATGCTAGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCCTTACACGCCTGATCAAACCCAACGCCACTACCTGAACGAGCTTCAACt         <  2:2655608/98‑1 (MQ=255)

                                                                                        cGCTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTa      >  1:1442907/1‑100 (MQ=255)

                                                                                          cTAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATc    >  1:1506531/1‑100 (MQ=255)

                                                                                           tAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGCAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCg   >  1:2588025/1‑100 (MQ=255)

                                                                                           tAATGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCg   >  1:2497714/1‑100 (MQ=255)

                                                                                            aaTGCTGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGt  <  2:1638336/100‑1 (MQ=255)

TTGCGTGTGGATCGCGCGCTTGCTGCTATGGGGCCTGCATTCCAACAGGTCTACAGTCTACTGCCGACATTGTTGCACTATCACCATCCGCTAATGCCGGGTTACCTTGATGGTAACGTTCCCAAAGGCATTTGCCTTTACACGCCTGATGAAACTCAACGCCACTACCTGAACGAGCTTGAACTGTATCGT  >  NC_000913/3991204‑3991395

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