BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,049	T→A	9.8%	intergenic (‑189/+10)	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,818,049	0	T	A	9.8%	666.5 / 32.1	208	intergenic (‑189/+10)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base T (108/77); new base A (7/15); total (116/92)
Fisher's exact test for biased strand distribution p-value = 2.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

CTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGA  >  NC_000913/2817901‑2818197

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

 taAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCACGCCCTTCGGGCCGTTGCTAAAGCAACGTTCTCATCCCGAGTGCCTGCGGTTAATTCCGCTACCACTTCGATACCTTCATTCCTGAAGATTTCg                                                                                                                                                    >  2:2618893/1‑150 (MQ=17)

        atatCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAata                                                                                                                                             >  1:2983237/1‑150 (MQ=255)

         tatCAATTAAGAGATGGTGCATCCGGGAGCAATCCTCGCCTAGCCGCCCGCCCCTCCAGCCTCTGAAAAAGAAACGTTATCATCCCTAGAGCCTGCGGTTAAATCCGCTACCTCTTCGATACCTGCATTGCTGAAGATTTCGgagactgt                                                                                                                                            >  1:2756065/1‑146 (MQ=18)

          atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                           >  2:4101299/1‑150 (MQ=255)

          atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                           <  1:1268960/150‑1 (MQ=255)

          atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                           >  2:4091019/1‑150 (MQ=255)

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

             aaTTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTg                                                                                                                                        <  1:2447761/150‑1 (MQ=255)

               ttAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTCTCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATATCGGAGAATATGGTGCa                                                                                                                                      >  1:2477006/1‑150 (MQ=255)

               ttAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCa                                                                                                                                      >  1:1161896/1‑150 (MQ=255)

                  agagATGGTGCATCCGGGAGGATTACTCCCCTTGCGGCTCGCCCTTCGGGCCGTTGCAAAAGCAACGATATCATACCTAGTGCCTGCGGTTAATTACGCTACCACTACGATACCTTCAATGCTGAAGATTTCGGAGAATATGGTGAATcc                                                                                                                                   >  1:2783814/1‑150 (MQ=18)

                      aTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGgga                                                                                                                               <  2:535498/150‑1 (MQ=255)

                         gTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggagga                                                                                                                            <  2:1875562/150‑1 (MQ=255)

                         gTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggagga                                                                                                                            <  1:1535672/150‑1 (MQ=255)

                            cATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTc                                                                                                                         <  1:2845122/150‑1 (MQ=255)

                             aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCg                                                                                                                        <  1:1505805/150‑1 (MQ=255)

                                cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAAc                                                                                                                     >  2:1594824/1‑150 (MQ=255)

                                cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGATACCTCTTCGATACCATCATTGCTGAACATTTCGGACAATATGGTCCATCCGGGACGATTAGAAc                                                                                                                     >  2:2444727/1‑150 (MQ=17)

                                  ggaggaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCt                                                                                                                   >  2:1895076/1‑150 (MQ=255)

                                   gaggaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTc                                                                                                                  >  1:1639319/1‑150 (MQ=255)

                                   gaggaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTc                                                                                                                  <  1:2639050/150‑1 (MQ=255)

                                      gatgccTCGCCTTGCGGCCCGCCCTTCGGGCAGTTGATAAAGCATCGTTACCAGCCCTAGTTCCTGCGGTTATTTACGCTACCTCTTCCATACGCTCTTCGCTGCAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCTAACCTCCGa                                                                                                               <  2:2578725/145‑1 (MQ=14)

                                      gaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                               <  2:3637213/150‑1 (MQ=255)

                                      gaTTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                               <  1:3651836/150‑1 (MQ=255)

                                            tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGc                                                                                                                                       <  2:552134/120‑1 (MQ=255)

                                            tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGc                                                                                                                                       >  1:552134/1‑120 (MQ=255)

                                            tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                         >  2:1744798/1‑150 (MQ=255)

                                              gCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCgg                                                                                                       >  2:3209302/1‑150 (MQ=255)

                                                 ttGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                                                    <  2:365890/150‑1 (MQ=255)

                                                 ttGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATAATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGGTc                                                                                                    >  1:1838981/1‑150 (MQ=255)

                                                       cTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                                                              >  1:800307/1‑150 (MQ=255)

                                                         cGCCCTTCGGGCCGGTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTTCTGAAGATTTCGGTGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                            <  1:470253/150‑1 (MQ=255)

                                                           cccTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGt                                                                                          >  1:2640691/1‑150 (MQ=255)

                                                            ccTTCGGGCCGTTGCTAAAGCAACGTTATCATACCTAGTGCCAGCGGTTAATTACGCTACCTCTTCGATAACATCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCACCGACCGCGCGGTTCGTAGCCGAGTa                                                                                         >  1:1923502/1‑150 (MQ=17)

                                                              ttCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATAAATTCATTGCTGAAGATTTCGGATAATATGGTGAATCCGGGATGATTAGAACATAAGACCGATCGGTTAGTAGCCGAGTAct                                                                                       >  2:1969920/1‑150 (MQ=18)

                                                              ttCGGGCCGTTGCTAAAGCAACGTTATCATCCCAAGTGCCAGCGGTTAATTCCGCAACCTCTTCGAAACCTTCATTTCTGAAGATTTCGGAGACAATGGAGCAACCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAcg                                                                                       >  1:3082416/1‑149 (MQ=17)

                                                                     cGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                                <  2:666618/150‑1 (MQ=255)

                                                                       ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGTTTCGTAGCCGAGTACTCTATCCAGc                                                                              >  1:1217005/1‑150 (MQ=255)

                                                                       ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                              >  1:2308834/1‑150 (MQ=255)

                                                                         gCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                                            >  1:2943387/1‑150 (MQ=255)

                                                                         gCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCAACCTCTTCGATACCTTCATTGCTGAAGAGTGCGGAGCGTCTGGTGCAGCCGGGAGGATTCGACCCTCCGACAGCTCGGTTCGGGGCCGAGAAATCTATCCCGCTg                                                                            >  2:1228856/1‑150 (MQ=17)

                                                                          cTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACATTCATTGCTGAAGATTTCGGAGACTATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCCGCTGa                                                                           >  2:2814284/1‑150 (MQ=255)

                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                                          >  1:1004708/1‑150 (MQ=255)

                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                                          >  1:2692248/1‑150 (MQ=255)

                                                                           tAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAg                                                                          >  1:3431482/1‑150 (MQ=255)

                                                                            aaaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                                                         >  2:1554828/1‑150 (MQ=255)

                                                                            aaaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGc                                                                         >  2:3916182/1‑150 (MQ=255)

                                                                               gCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGTATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAc                                                                      >  2:2300362/1‑150 (MQ=255)

                                                                                cAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACg                                                                     >  1:3318401/1‑150 (MQ=255)

                                                                                  aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                                   >  2:4002319/1‑150 (MQ=255)

                                                                                  aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGCACTCTATCCAGCTGAGCTACGGa                                                                   >  2:1041543/1‑150 (MQ=255)

                                                                                  aCGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCCTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGa                                                                   <  2:3377547/150‑1 (MQ=255)

                                                                                          atcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                           >  1:567071/1‑150 (MQ=255)

                                                                                          atcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCgg                                                           >  1:200019/1‑150 (MQ=255)

                                                                                              cTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGaaa                                                       <  2:3119851/150‑1 (MQ=255)

                                                                                                            tAATACCTCTATCTCTTCGATACCTTCATTGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTActgct                                         <  1:1406083/150‑1 (MQ=255)

                                                                                                              aTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTTTTTTACTGCTGa                                       <  1:3209302/150‑1 (MQ=255)

                                                                                                              aTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                       <  1:1611152/150‑1 (MQ=255)

                                                                                                                tCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                                               >  2:3607803/1‑108 (MQ=40)

                                                                                                                 ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                                    <  2:3505853/150‑1 (MQ=255)

                                                                                                                  cGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATAtt                                   <  2:2917123/150‑1 (MQ=255)

                                                                                                                   gCTATCTCTTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTg                                  <  2:1317912/150‑1 (MQ=9)

                                                                                                                    cTATCTCTTCGAGACCTTCATTGCTGAAGTTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGa                                 <  2:3326878/150‑1 (MQ=255)

                                                                                                                    cTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTATTTTACTGCTGATATTGa                                 >  2:3940514/1‑150 (MQ=255)

                                                                                                                     tACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGAt                                <  1:2157451/150‑1 (MQ=255)

                                                                                                                       tctctTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAc                              <  1:3770592/149‑1 (MQ=17)

                                                                                                                        ctctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTACTTTACTGCAGATATTtataaa                             >  2:483400/1‑148 (MQ=255)

                                                                                                                         tctTCGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGAt                                <  1:396767/146‑1 (MQ=18)

                                                                                                                             cGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGctac                        <  2:1161896/150‑1 (MQ=255)

                                                                                                                                tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaa                     <  2:1004708/150‑1 (MQ=255)

                                                                                                                                 aCCTTCATTGCTGAAGTTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                    <  1:3543304/150‑1 (MQ=25)

                                                                                                                                 aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                    <  2:1932254/150‑1 (MQ=18)

                                                                                                                                 aCCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                    <  2:1761626/150‑1 (MQ=18)

                                                                                                                                 aCCTTCATTGCTGAAGATTTCTGAGAATATGGTGCATCCGTGAGGATTCGTACCTCCGTCCGCTCTGTTCGTAGCCGAGTTCTCTATCAATCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaaa                    <  1:3282700/150‑1 (MQ=14)

                                                                                                                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc                 >  2:218373/1‑150 (MQ=255)

                                                                                                                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc                 >  1:2441247/1‑150 (MQ=255)

                                                                                                                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc                 <  1:3077467/150‑1 (MQ=255)

                                                                                                                                    ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACAGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGATACGGATGCATCGGGACATTTATTTTACTGCTGATATTGATACCGCTAGTAAAGcc                 >  1:2809671/1‑150 (MQ=255)

                                                                                                                                     tCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCg                                                                                                        <  1:1597949/62‑1 (MQ=40)

                                                                                                                                     tCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCg                                                                                                        >  2:1597949/1‑62 (MQ=40)

                                                                                                                                         tGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATc            <  2:2692248/150‑1 (MQ=255)

                                                                                                                                         tGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATc            <  1:4101299/150‑1 (MQ=255)

                                                                                                                                           cTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa          <  2:1118074/150‑1 (MQ=18)

                                                                                                                                           cTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCACCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATGTATTTTACTGCTGATATTGATACCGCTACTAAAGCCCTGTCaa          >  2:1325800/1‑150 (MQ=11)

                                                                                                                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa          >  2:820843/1‑150 (MQ=255)

                                                                                                                                           cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa          >  1:1183981/1‑150 (MQ=255)

                                                                                                                                           atGAAGTTTACGTTGAATTTGGTGCTTCGGGTCGTATTCTAACCTCGGACCGCAGGGTTGGTTGCGGTTTACTGTATCCAGCTGTGATACGGATGCATGGGGAAATTTATTCTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa          <  2:1674511/149‑1 (MQ=2)

                                                                                                                                            tGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg         <  1:6222/150‑1 (MQ=21)

                                                                                                                                            tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg         >  1:1734268/1‑150 (MQ=255)

                                                                                                                                            tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg         <  2:2308834/150‑1 (MQ=255)

                                                                                                                                             gAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGt        >  2:1856184/1‑150 (MQ=255)

                                                                                                                                              aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTCAAGCCATATCAAGTa       >  2:2394431/1‑150 (MQ=255)

                                                                                                                                              aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTAGTTTACTGCTGATATTTTTACAGCTACTATAGCCAAATCaataa       >  1:653015/1‑147 (MQ=11)

                                                                                                                                               aGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGTATTCATCTGGAAATTTATTGTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTaa      <  2:1247953/150‑1 (MQ=255)

                                                                                                                                               aGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTaa      <  2:91654/150‑1 (MQ=255)

                                                                                                                                               aGATTTCGGAGAATAGGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTaa      <  2:2820383/150‑1 (MQ=255)

                                                                                                                                                 gtttCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                                                         <  1:1042209/64‑1 (MQ=18)

                                                                                                                                                 aTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTa                                                                                         <  1:776734/65‑1 (MQ=37)

                                                                                                                                                 aTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAaga    <  2:2051987/150‑1 (MQ=255)

CTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGA  >  NC_000913/2817901‑2818197

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: `‑`