Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,049	T→A	8.0%	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	8.0%	318.3 / 12.1	97	intergenic (‑189/+10)	argQ/argZ	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (55/34);  new base A (4/4);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTA  >  NC_000913/2817908‑2818192                                                                                                                                              |                                                                                                                                                cATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAat                                                                                                                                         <  2:901350/150‑1 (MQ=255)    atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                      <  2:2161415/150‑1 (MQ=255)    atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTCGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                      <  2:383397/150‑1 (MQ=255)           aaGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATc                                                                                                                               <  2:2135767/150‑1 (MQ=255)           aaGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATc                                                                                                                               >  1:607660/1‑150 (MQ=255)           aaGAGATGGTGCATCCGGGAGGATTACTCGCCGTGCGGCTCGCCCTACGGGCCGTTGCTATAGCAACGTTATCATCGCTATTGCGTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATc                                                                                                                               <  2:1616862/150‑1 (MQ=18)             gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCg                                                                                                                             >  1:1507880/1‑150 (MQ=255)             gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCg                                                                                                                             >  1:1507836/1‑150 (MQ=255)             gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCg                                                                                                                             >  1:1479392/1‑150 (MQ=255)             gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGATGCTAAAGCAACGATATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTAGATACCTTCATTGCTGAAGAATTCGGAGAATATTGTGCATCCg                                                                                                                             >  1:2522195/1‑150 (MQ=17)              agaTGGTGCATCCGGGAGGATTACTCGCGTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCTTACCTAGTGCCTGCGGTTAATTCCGCTAGCTCTTCGATACCTTCATTGCTGTAGTTTTCGGAGATTATGGTGCTTCCgg                                                                                                                            >  2:1249366/1‑150 (MQ=17)              agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:1034743/1‑150 (MQ=255)              agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  2:2473761/1‑150 (MQ=255)              agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:2314650/1‑150 (MQ=255)              agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:1877708/1‑150 (MQ=255)               gaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCACCGGTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTGCATAGCTGAAGATTTCGGAGAATATGGTGCATCCggg                                                                                                                           >  1:1790773/1‑150 (MQ=255)               gaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCggg                                                                                                                           >  1:1922603/1‑150 (MQ=255)               gaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCggg                                                                                                                           >  1:125828/1‑150 (MQ=255)                aTGGTGCACCCGGGAGGCTAACTCGCCCTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGgga                                                                                                                          <  2:1441510/150‑1 (MQ=255)                   gTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggagga                                                                                                                       >  2:1560254/1‑150 (MQ=255)                       aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCg                                                                                                                   >  1:297344/1‑150 (MQ=255)                         ccGGGAGGTGTTCTCGCCCTGCGGCTCGCCCTTCGGTCCGTTGCTTTAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGaa                                                                                                                 <  2:560233/150‑1 (MQ=18)                         ccGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTCCAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGaa                                                                                                                 <  2:656769/150‑1 (MQ=255)                          cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGACTATGGTGCATCCGGGAGGATTCGAAc                                                                                                                >  1:1851056/1‑150 (MQ=255)                           gggAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAAcc                                                                                                               >  1:1558479/1‑150 (MQ=255)                                    aCTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGc                                                                                                      <  2:1690662/150‑1 (MQ=255)                                      tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:776651/1‑150 (MQ=255)                                      tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:776745/1‑150 (MQ=255)                                      tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATCCGAACCTCCGACCGCTc                                                                                                    >  1:1939214/1‑150 (MQ=255)                                      tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATCCGAACCTCCGACCGCTc                                                                                                    >  1:1939247/1‑150 (MQ=255)                                      tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCAGCGGTTAATTCCGCTACCTCTTAGATACCTTCATTGCTGAAGATTTCGGAGAATATTGTGCATCCGGGAGGATCCGAACCTCCGACCGCTc                                                                                                    >  1:659252/1‑150 (MQ=255)                                           ttGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGGTc                                                                                               >  1:2530891/1‑150 (MQ=255)                                           ttGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATACCTAGTGCCTGCGGTTAATTACGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                                               >  1:2531024/1‑150 (MQ=255)                                                 ttCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga                                                                                                                                             >  2:2091546/2‑98 (MQ=255)                                                 ttCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga                                                                                                                                             <  1:2091546/97‑1 (MQ=255)                                                   cGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                       <  1:737596/150‑1 (MQ=255)                                                       cTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAc                                                                                   >  1:1018292/1‑150 (MQ=255)                                                           gggCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                               >  2:537518/1‑150 (MQ=255)                                                             gCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                             <  1:2021512/150‑1 (MQ=255)                                                                 ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                         >  2:1632641/1‑150 (MQ=255)                                                                 ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                         >  2:1203503/1‑150 (MQ=255)                                                                 ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                         >  2:1203482/1‑150 (MQ=255)                                                                       aaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                                   >  1:2012312/1‑150 (MQ=255)                                                                       aaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                                   >  1:138226/1‑150 (MQ=255)                                                                              gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                                                            >  1:2167543/1‑150 (MQ=255)                                                                                  tcatcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATc                                                        <  1:409864/150‑1 (MQ=255)                                                                                           gTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAattt                                               >  1:1091984/1‑150 (MQ=255)                                                                                           gTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAagtt                                               >  1:2258606/1‑150 (MQ=255)                                                                                                 gCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                         >  1:781070/1‑150 (MQ=255)                                                                                                  cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                        >  2:221285/1‑150 (MQ=255)                                                                                                        aTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                  <  2:2547131/150‑1 (MQ=255)                                                                                                           ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               <  2:1507836/150‑1 (MQ=255)                                                                                                           ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               <  2:1507880/150‑1 (MQ=255)                                                                                                           ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               <  2:1479392/150‑1 (MQ=255)                                                                                                            cGCTACCTCTTCGATACCTTCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGa                                                                                                                  <  1:322617/66‑1 (MQ=25)                                                                                                             gCTAGCTCTTCGATACCATCATTTCTGAAGATTTAGGATAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAttttt                             <  2:2531024/150‑2 (MQ=11)                                                                                                                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGACACTTACCTTACTGCTGATATTGAGAcac                       >  2:1171853/1‑148 (MQ=11)                                                                                                                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                       >  2:2395011/1‑150 (MQ=255)                                                                                                                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                       >  2:346433/1‑150 (MQ=255)                                                                                                                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACATCCGACCGCTCGGTTCGTAGCCGAGTACACAATACAGATGAGCTACGGATGCATCGGGAAATTTATGTTACTGCTGATATTGTTACCg                       >  1:2256093/1‑150 (MQ=14)                                                                                                                   tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACAGCCGACCGCTCGGTTCTTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGAATCGGGAAAGTTATTTTACTGCTGATATTGATAGCg                       >  1:1613155/1‑150 (MQ=14)                                                                                                                     ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGt                                                                                                                                    <  1:2020280/39‑1 (MQ=255)                                                                                                                     ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTATTTTACTGCTGATATTGATACCGct                     >  2:2159346/1‑150 (MQ=255)                                                                                                                      tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:811619/150‑1 (MQ=255)                                                                                                                      tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:171045/150‑1 (MQ=255)                                                                                                                      tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:579508/150‑1 (MQ=255)                                                                                                                       cGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGctac                   <  1:1973602/150‑1 (MQ=18)                                                                                                                          tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCtaataa                >  1:1038987/1‑150 (MQ=255)                                                                                                                          tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaa                >  1:431941/1‑150 (MQ=255)                                                                                                                          tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACAGCTACTaa                >  1:1059327/1‑150 (MQ=255)                                                                                                                          tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTAATTTACTGCTGATATTGTTAACGCtaataa                >  1:1059343/1‑150 (MQ=14)                                                                                                                            cttcgtttgttgAAGTTTTCGGAGAATATGGTGCATCGGGGAGGATTCGAACCTCCGACCGCTGGGTTCGTAGCCGTGTACTCTATCCAGCTGAGCTACGGATGCATGGGTAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:1143430/144‑1 (MQ=11)                                                                                                                            ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:1143429/150‑1 (MQ=255)                                                                                                                            ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:473396/150‑1 (MQ=255)                                                                                                                              ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTCCGGATGCATCGGGAAATATATTTGACGGCTGATATTGATACCGCTACTAAAGcc            >  2:2526085/1‑150 (MQ=11)                                                                                                                              ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc            >  2:1025157/1‑150 (MQ=18)                                                                                                                              ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc            >  2:191668/1‑150 (MQ=18)                                                                                                                              ttCATTGCTGAAGATTTCGGAGTATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                       <  2:1370478/107‑1 (MQ=25)                                                                                                                              ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                       >  1:1370478/1‑107 (MQ=35)                                                                                                                              ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc            >  1:1600896/1‑150 (MQ=255)                                                                                                                               tCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCa           >  1:1402424/1‑150 (MQ=255)                                                                                                                                  ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCatat        <  1:1203482/150‑1 (MQ=255)                                                                                                                                  ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCatat        <  1:1203503/150‑1 (MQ=255)                                                                                                                                  ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCatat        <  1:1632641/150‑1 (MQ=255)                                                                                                                                    ccTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCa      <  1:717246/149‑1 (MQ=18)                                                                                                                                     cTGAAGATTTCTGTGTATTTGGTGCTTGGGGGATGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa     <  1:1463915/150‑1 (MQ=14)                                                                                                                                     cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTTATAACGCTAATAAAGCCATATCaa     >  1:933721/1‑150 (MQ=255)                                                                                                                                     cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa     >  2:138882/1‑150 (MQ=255)                                                                                                                                     cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCACAGACCCCACCCGTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAGCGCGCAATTTATGTTCCGGCGGATAATGATACCGATACTCAAGCAATATcta     >  1:933781/1‑148 (MQ=9)                                                                                                                                      tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                          <  1:2398504/48‑1 (MQ=255)                                                                                                                                      tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg    <  2:1558479/150‑1 (MQ=255)                                                                                                                                       gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              >  2:1143464/1‑139 (MQ=18)                                                                                                                                       gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:1143464/139‑1 (MQ=18)                                                                                                                                       gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGt   <  1:2504094/150‑1 (MQ=25)                                                                                                                                        aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCt                                                                                                     >  1:1987118/1‑51 (MQ=255)                                                                                                                                        aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCt                                                                                                     <  2:1987118/51‑1 (MQ=255)                                                                                                                                        aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAACGCTACAAAAGCCATATCAAGTa  >  2:1777102/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                CATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTA  >  NC_000913/2817908‑2818192

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