BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,817,931	T→C	47.8%	intergenic (‑71/+128)	argQ ← / ← argZ	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,817,931	0	T	C	47.8%	0.5 / 215.8	156	intergenic (‑71/+128)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (42/38); new base C (35/39); total (78/78)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-02

Only 100 of 158 total aligned reads displayed.

GTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATA  >  NC_000913/2817833‑2818030

gTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATc                                                                                                    >  2:3058253/1‑100 (MQ=255)

 tACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCAccc                                                                                                   <  2:731422/100‑1 (MQ=255)

 tACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATcc                                                                                                   <  1:1798448/100‑1 (MQ=255)

   ctctATCCAGCTGAGCTACGGATGCATCGGGAAATTTTTTTCTACTGCTGATTTTGATACCCCTAGTAAAGCCATATAAATTAAAAGATGTTGCA‑CCggg                                                                                                >  1:991847/1‑100 (MQ=11)

    tctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCggg                                                                                                >  2:1139357/1‑100 (MQ=255)

      tATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGggag                                                                                              <  1:2399547/100‑1 (MQ=255)

       aTCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATTTTGATACCGCTACTAAAACCATATCAATTAAGAGATGTTGCACCCggg                                                                                                >  1:2500567/1‑97 (MQ=17)

            gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTAc                                                                                        >  2:4275580/1‑100 (MQ=255)

            gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTAc                                                                                        <  2:3947299/100‑1 (MQ=255)

            gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTAc                                                                                        <  1:917768/100‑1 (MQ=255)

              tGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTc                                                                                      <  1:1444079/100‑1 (MQ=255)

               gAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCg                                                                                     >  2:2527150/1‑100 (MQ=255)

                 gCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGcc                                                                                   >  2:2664231/1‑100 (MQ=255)

                   tACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGAttcgtagc                                                                                    >  1:1736404/1‑91 (MQ=255)

                   tACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCtt                                                                                 >  2:3562159/1‑100 (MQ=255)

                    aCGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTg                                                                                >  2:1619253/1‑100 (MQ=255)

                    aCGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGAATCCGGGAGGATTACTCGCCtt                                                                                 <  1:2474169/99‑1 (MQ=255)

                     cGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGc                                                                               <  2:2942468/100‑1 (MQ=255)

                      ggATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCg                                                                              <  2:4168694/100‑1 (MQ=255)

                       gATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCgg                                                                             >  1:3927759/1‑100 (MQ=255)

                          gCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTgag                                                                              >  1:3284476/1‑94 (MQ=255)

                          gCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTc                                                                          >  1:543687/1‑100 (MQ=255)

                           cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCg                                                                         >  1:4168847/1‑100 (MQ=255)

                           cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCg                                                                         <  1:100670/100‑1 (MQ=255)

                             tCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGcc                                                                       <  1:46225/100‑1 (MQ=255)

                             tCGGAAAAATTTATTCTTGTAGTGATTTTGATACCGCTACTAAAGCCGTATCAAGTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCg                                                                              <  1:1588190/93‑1 (MQ=11)

                               gggAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCTTTGCGGCTCGCCCt                                                                     >  2:2856036/1‑100 (MQ=255)

                               gggAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTAATCGCCTCGCGGCCCGCCCt                                                                     >  1:109123/1‑100 (MQ=255)

                                 gAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTc                                                                   >  1:2946961/1‑100 (MQ=255)

                                    aaCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCggg                                                                <  2:2366243/100‑1 (MQ=255)

                                     aCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGc                                                               >  1:2432192/1‑100 (MQ=255)

                                     aCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGc                                                               >  2:1806781/1‑100 (MQ=255)

                                     aCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGc                                                               >  2:1806749/1‑100 (MQ=255)

                                     aCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGc                                                               <  2:527764/100‑1 (MQ=255)

                                      cTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGcc                                                              >  2:1063189/1‑100 (MQ=255)

                                      cTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGc                                                               <  2:1248658/99‑1 (MQ=255)

                                       ttATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCg                                                             <  2:2726613/100‑1 (MQ=255)

                                        tATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGt                                                            >  2:749876/1‑100 (MQ=255)

                                          ttCTACTGTTGATATTGATACCGCTCCTAAAGCAATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTg                                                          <  1:59619/100‑1 (MQ=255)

                                           tCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGc                                                         <  1:2744576/100‑1 (MQ=255)

                                           tCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGc                                                         <  2:3881095/100‑1 (MQ=255)

                                            cTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCt                                                        >  2:228034/1‑100 (MQ=255)

                                            cTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCt                                                        >  2:2237396/1‑100 (MQ=255)

                                              aCTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaa                                                      >  2:1648757/1‑100 (MQ=255)

                                              aCTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCACGCCCTTCGGGCCGTTGCTa                                                       <  1:749874/99‑1 (MQ=255)

                                               ctgctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                     <  2:4106962/100‑1 (MQ=255)

                                                 gctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGc                                                   >  1:1217959/1‑100 (MQ=255)

                                                  ctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCa                                                  >  1:1717932/1‑100 (MQ=255)

                                                  ctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCa                                                  >  1:218586/1‑100 (MQ=255)

                                                  ctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATAGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCa                                                  >  2:1007941/1‑100 (MQ=255)

                                                   tgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCCCCTCGCGCCTCGGCCTTCGGGCCGTTTCTaaa                                                     >  1:1996224/1‑96 (MQ=14)

                                                    gATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAc                                                >  1:2300961/1‑100 (MQ=255)

                                                    gATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAAc                                                >  2:2983350/1‑100 (MQ=255)

                                                      tatTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGt                                              <  2:3549418/100‑1 (MQ=40)

                                                      tatTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGt                                              >  2:1697212/1‑100 (MQ=255)

                                                       atTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGtt                                             >  2:1615566/1‑100 (MQ=255)

                                                        ttGATACCGCTACTAAAGCCATATCAATTAAGAGGTGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTa                                            >  1:1759275/1‑100 (MQ=255)

                                                        ttGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCGTTCGGGCCGTTGCTTTAGCAACggcc                                            <  2:2592734/100‑4 (MQ=11)

                                                        ttGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTa                                            <  2:461003/100‑1 (MQ=255)

                                                         tGATACCGCTACTAGAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTTTAGCAACggccc                                           <  1:1238659/100‑5 (MQ=11)

                                                         tGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTTTAGCAACggccc                                           <  1:2163305/100‑5 (MQ=14)

                                                          gATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatc                                          >  1:253680/1‑100 (MQ=38)

                                                          gATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatc                                          >  1:3102247/1‑100 (MQ=38)

                                                           aTACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatca                                         <  1:5915/100‑1 (MQ=40)

                                                           aTACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTTTAGCAACggcccga                                         <  1:2823688/100‑7 (MQ=14)

                                                             aCCGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTatcatc                                       <  2:2219705/100‑1 (MQ=255)

                                                              ctgCTACTAGAGCCATATCAATTAGGAGATGGTGCACCCGGGAGGAGTACTCGCATTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATcc                                      <  1:3724908/98‑1 (MQ=14)

                                                              ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATcc                                      >  2:946784/1‑100 (MQ=255)

                                                              ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATcc                                      >  2:3307187/1‑100 (MQ=255)

                                                              ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTTTAGCAACggcccgaagg                                      <  2:3160846/100‑10 (MQ=14)

                                                               cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCCTccc                                     >  1:1055451/1‑100 (MQ=25)

                                                                  tactaAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAg                                  <  2:818522/100‑1 (MQ=255)

                                                                   actaAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCCTCCCTGGt                                 >  1:430393/1‑100 (MQ=255)

                                                                   aataAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCCAAAGCCACGTTATCATCCCTAGt                                 >  2:2890036/3‑100 (MQ=255)

                                                                     taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGc                               <  1:4098347/100‑1 (MQ=255)

                                                                     taAAGCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGc                               >  2:1442447/1‑100 (MQ=255)

                                                                         gCCATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGc                           <  2:2546614/100‑1 (MQ=255)

                                                                           cATATCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCgg                         >  2:2673211/1‑100 (MQ=255)

                                                                            atatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGTGGt                        <  2:1003684/100‑1 (MQ=255)

                                                                            atatCAATTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGt                        <  1:854773/100‑1 (MQ=255)

                                                                             tatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                       >  1:637659/1‑100 (MQ=255)

                                                                             tatCAATCAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                       >  2:1485539/1‑100 (MQ=255)

                                                                              atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTa                      <  2:1872902/100‑1 (MQ=255)

                                                                               tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa                     >  2:186712/1‑100 (MQ=40)

                                                                                cAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAAt                    <  2:4172115/100‑1 (MQ=255)

                                                                                  aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCTTTGCGGCTCGCGCCTCGGGCCGTGTCTAAAGCAAGGTTATCATCCCTAGTGCCTGCGGTTAGtt                   >  1:973096/1‑99 (MQ=14)

                                                                                  aTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTc                  <  2:594702/100‑1 (MQ=255)

                                                                                  aTTAAGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTc                  >  1:643435/1‑100 (MQ=255)

                                                                                     aaGAGATGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGc               >  1:3059404/1‑100 (MQ=255)

                                                                                      agagATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCt              >  2:2504576/1‑100 (MQ=255)

                                                                                       gagaTGGTGCATCCGGGAGGATTAGTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGc               <  1:1036782/98‑1 (MQ=255)

                                                                                       gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTagatcggaag             >  1:3896987/1‑91 (MQ=255)

                                                                                        agaTGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCCTCCCTAGTGCCTGCGGTTAATTCCGCCAc            >  1:2391813/1‑100 (MQ=255)

                                                                                        agaTGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAc            <  1:3520707/100‑1 (MQ=255)

                                                                                        agaTGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAc            >  2:1698721/1‑100 (MQ=255)

                                                                                         gaTGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAcc           <  1:1629956/100‑1 (MQ=255)

                                                                                         gaTGGTGCACCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAcc           <  2:2451702/100‑1 (MQ=255)

                                                                                              tGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTc      <  2:1904343/100‑1 (MQ=255)

                                                                                               gCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCg     <  2:3827361/100‑1 (MQ=255)

                                                                                                  tCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATa  >  2:4155732/1‑100 (MQ=255)

GTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATA  >  NC_000913/2817833‑2818030

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