BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,817,916	T→G	6.7%	intergenic (‑56/+143)	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,916	0	T	G	6.7%	1116.1 / 60.3	346	intergenic (‑56/+143)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (134/189); new base G (11/12); total (145/201)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

Only 100 of 348 total aligned reads displayed.

AGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGC  >  NC_000913/2817769‑2818064

  gTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAag                                                                                                                                                  <  1:2836034/150‑1 (MQ=255)

  gTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGCTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAag                                                                                                                                                  >  2:1039928/1‑150 (MQ=255)

    tACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAagag                                                                                                                                                <  1:4191674/150‑1 (MQ=255)

       ggAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATG                                                                                                                                             <  1:934978‑M1/150‑1 (MQ=255)

         agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGT                                                                                                                                           <  2:6800337‑M1/150‑1 (MQ=255)

         agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGT                                                                                                                                           <  2:3028215‑M1/150‑1 (MQ=255)

             tatGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCAT                                                                                                                                       >  1:3675483‑M1/1‑150 (MQ=255)

                   gCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGA                                                                                                                                 <  2:2043096‑M1/150‑1 (MQ=255)

                   gCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGA                                                                                                                                 <  1:299798‑M1/150‑1 (MQ=255)

                      tCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGA                                                                                                                              <  1:4213971‑M1/150‑1 (MQ=255)

                      tCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGA                                                                                                                              <  1:2528826‑M1/150‑1 (MQ=255)

                          ggaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCAACTAAAGCCATATCAATTAAGAGATGGTGCAGCCGGGAGGATTcc                                                                                                                          >  2:6173750‑M1/1‑148 (MQ=255)

                           gaggaTTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACt                                                                                                                         >  2:5775757/1‑150 (MQ=255)

                                ttCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCt                                                                                                                    >  1:5798966/1‑150 (MQ=255)

                                 tCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGA                                                                                                                                 >  2:794432‑M1/1‑136 (MQ=255)

                                     aCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCgg                                                                                                               <  1:1396231/150‑1 (MQ=255)

                                     aCCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCgg                                                                                                               <  1:4433104/150‑1 (MQ=255)

                                      ccTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGc                                                                                                              >  1:7068067/1‑150 (MQ=255)

                                      ccTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGc                                                                                                              >  1:3496671/1‑150 (MQ=255)

                                        tCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTc                                                                                                            <  1:3370313/150‑1 (MQ=255)

                                        tCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTc                                                                                                            <  1:1601649/150‑1 (MQ=255)

                                         cagACCGCTCGGTTCGTATCCGAGTACTCTATCCCTCTGAGCTAAGGATGCATCGGGACAACTTATTATTCTGCTTATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCg                                                                                                           <  2:5923890/148‑1 (MQ=17)

                                           gACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGcc                                                                                                         <  1:2189721/150‑1 (MQ=255)

                                             ccGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCt                                                                                                       >  2:804103/1‑150 (MQ=255)

                                              cGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCtt                                                                                                      <  1:6489555/150‑1 (MQ=255)

                                                 tCGGTTCGTAGCCGAGTACTCTATCCAGCTGATCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   <  1:5388136/150‑1 (MQ=255)

                                                 tCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCgg                                                                                                   <  1:3912428/150‑1 (MQ=255)

                                                  cGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCggg                                                                                                  <  2:6366673/150‑1 (MQ=255)

                                                     ttCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAAATTATTCTAATGATGATATTGATACCGCTAATAATGCCATAACAATTAAGAGATGGAGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCg                                                                                               >  2:132504/1‑150 (MQ=17)

                                                      tCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGt                                                                                              <  1:5582916/150‑1 (MQ=255)

                                                       cGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGtt                                                                                             <  2:572568/150‑1 (MQ=255)

                                                            ccGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaa                                                                                        >  2:7069880/1‑150 (MQ=255)

                                                            ccGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCCGTGCGGGCCGTTGCTaa                                                                                        >  1:4030836/1‑150 (MQ=255)

                                                             cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       <  2:2266732/150‑1 (MQ=255)

                                                             cGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTaaa                                                                                       <  1:1145467/150‑1 (MQ=255)

                                                              gAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAg                                                                                      <  1:3113191/150‑1 (MQ=255)

                                                                   ctctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccag                                                                                 <  2:32201‑M1/150‑44 (MQ=255)

                                                                    tctATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTagaaccacagaccgctcggttcgtggccgagtactctatccagc                                                                                >  2:655753‑M1/1‑106 (MQ=255)

                                                                       aTCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTat                                                                             >  1:5320426/1‑150 (MQ=255)

                                                                            gCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATcc                                                                        >  1:2617148/1‑150 (MQ=255)

                                                                             cTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATccc                                                                       <  1:1823332/150‑1 (MQ=255)

                                                                               gaccttcGGATGCATCGGGAAAACTTATACTACTGCTTTAATTGATACATATACTAAAGACATATAATTTATGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTa                                                                     <  2:6995789/144‑1 (MQ=17)

                                                                                  cTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTg                                                                  >  1:1788801/1‑150 (MQ=255)

                                                                                  cTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTg                                                                  >  1:4942149/1‑150 (MQ=255)

                                                                                   tACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGc                                                                 >  1:1934394/1‑150 (MQ=255)

                                                                                    aCGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGcc                                                                <  2:2870211/150‑1 (MQ=255)

                                                                                     cGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCt                                                               >  1:5277887/1‑150 (MQ=255)

                                                                                     cGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCt                                                               >  1:777706/1‑150 (MQ=255)

                                                                                      ggATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCCCTGCGGGCCGTTGCTAAAGCAACGGTATCATCCCTAGTGCCTg                                                              >  2:5113005/1‑150 (MQ=255)

                                                                                      ggATGCATCGGGAAAACATATACTACTTCTGATATTGATACCTATACTAAAGCCATATAAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTg                                                              <  2:5256936/150‑1 (MQ=18)

                                                                                        aTGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCg                                                            >  1:1356553/1‑150 (MQ=255)

                                                                                           cATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCATTGCGGCTCGCCCTACGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGtt                                                         >  2:3494656/1‑150 (MQ=255)

                                                                                            aTCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATGTCAATTGAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGGTGATAATGCACCGTTATCATCCGTAGTGCCTGCGGTTa                                                        >  1:1231243/1‑150 (MQ=18)

                                                                                            aTCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTa                                                        <  1:3268810/150‑1 (MQ=255)

                                                                                            aTCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTa                                                        >  1:6868442/1‑150 (MQ=255)

                                                                                            aTCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTAGCGGCTCGCCCTTCGGGCCGTTGCGAAAGCACCGTTATCATCCCTAGTGCATGCGGtga                                                        >  2:4498605/1‑148 (MQ=255)

                                                                                             tctgtacaaCTTATTCTAATTCTGATATTGATACAGCTACTAAAGCCATATCATTTATTAGATGGTGCATACGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa                                                       <  1:429711/143‑1 (MQ=17)

                                                                                             tCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa                                                       <  1:6281813/150‑1 (MQ=255)

                                                                                             tCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTaa                                                       <  2:2115881/150‑1 (MQ=255)

                                                                                               gggAATACTTATTCTACTGCTTATATTGTTAACGATAATATAGCCATTTTAATTATTAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAAtt                                                     <  1:1756393/150‑1 (MQ=17)

                                                                                               gggAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCTCGTTGGGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAAtt                                                     >  1:5397907/1‑150 (MQ=255)

                                                                                                 ggaaattTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggc                                                   >  2:1515581‑M1/3‑77 (MQ=255)

                                                                                                 gAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTcc                                                   >  1:2323598/1‑150 (MQ=255)

                                                                                                  aaaaCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCg                                                  <  1:6032242/150‑1 (MQ=255)

                                                                                                   aaaCTTATTCTACTGCTGATATTGATACCTCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGc                                                 <  1:6810901/150‑1 (MQ=255)

                                                                                                   aaaCTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGc                                                 <  1:2692337/150‑1 (MQ=255)

                                                                                                      tgtatTTTACAGCTGATTCTGATACCGCTACTAACGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcgggga                                              <  2:4273957‑M1/148‑79 (MQ=255)

                                                                                                       ttATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaactggcggggagccggggat                                             >  2:315924‑M1/1‑71 (MQ=255)

                                                                                                       ttATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTAcc                                             >  1:2280848/1‑150 (MQ=255)

                                                                                                            cTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTc                                        <  2:608568/150‑1 (MQ=255)

                                                                                                            cTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTc                                        <  1:2067532/150‑1 (MQ=255)

                                                                                                              aCTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGa                                      >  2:1585792/1‑150 (MQ=255)

                                                                                                                tgctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATa                                    >  2:582692/1‑150 (MQ=255)

                                                                                                                  ctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATT                                                                                                                            <  1:146662‑M1/60‑1 (MQ=255)

                                                                                                                  ctgATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGA                                                                                                                                 <  1:6351822‑M1/55‑1 (MQ=255)

                                                                                                                       atTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCAt                             >  1:2447929/1‑150 (MQ=255)

                                                                                                                       atTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccc                                     <  1:6405154‑M1/142‑88 (MQ=255)

                                                                                                                          gATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGc                          >  1:4681764/1‑150 (MQ=255)

                                                                                                                            tACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTg                        <  1:820446/150‑1 (MQ=255)

                                                                                                                             aCCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGa                       <  1:5531905/150‑1 (MQ=255)

                                                                                                                              ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                      <  2:2403513/150‑1 (MQ=255)

                                                                                                                              ccGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGaa                      <  1:3894784/150‑1 (MQ=255)

                                                                                                                               cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGc                                                                                                              <  1:5719030/61‑1 (MQ=40)

                                                                                                                               cGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATT                                                                                                                            >  2:4410542‑M1/1‑47 (MQ=255)

                                                                                                                               cGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtgaggcggggattcgaaccccggatgcagctt                         >  1:3772419‑M1/1‑47 (MQ=255)

                                                                                                                                gCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTACTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGa                    >  2:7143857/1‑150 (MQ=255)

                                                                                                                                gCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCATAGAGGCTCGCCCTTCGTGCCGTTGCTAAAGAAACGTAATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGa                    >  2:3555180/1‑150 (MQ=17)

                                                                                                                                gCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTcgaacctccgaccgctcggttcgtagccgagtactctatccagctgagctacggatgcaaatggcggtga                                                      >  2:1812029‑M1/1‑46 (MQ=255)

                                                                                                                                 ctactaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAt                   >  2:2636931/1‑150 (MQ=255)

                                                                                                                                   actaCAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCGTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAttt                 <  2:5609911/150‑1 (MQ=255)

                                                                                                                                   actaAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGAttt                 <  2:5549515/150‑1 (MQ=255)

                                                                                                                                     taACGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg               <  1:5698707/150‑1 (MQ=255)

                                                                                                                                     taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCg               <  1:6104504/150‑1 (MQ=255)

                                                                                                                                       taGCCATATCATTTAGTTGATGGTGCTCCCGGGAGGTTTAGTCGCCTTGAGGCTCGCCCTTCGGGCCGTTGCTACAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga             <  2:4845896/149‑1 (MQ=18)

                                                                                                                                       acgCCATATCAATTCAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTGCGGGCCGTGGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga             <  1:4072095/148‑1 (MQ=255)

                                                                                                                                       aaGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga             <  2:384835/150‑1 (MQ=255)

                                                                                                                                        agaaatatAAATTAAGCGATTGTGCATCCGGGAGGATTACTAGCCTTGCGTCTCGCACTTCTGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCTGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGgag            <  2:3903968/146‑1 (MQ=18)

                                                                                                                                         gCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTAAAATCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGgaga           >  2:3751164/1‑150 (MQ=255)

                                                                                                                                          ccATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGaa          >  1:3202078/1‑150 (MQ=255)

                                                                                                                                               tCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATgg     <  1:1039613/150‑1 (MQ=255)

AGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAACTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGC  >  NC_000913/2817769‑2818064

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