Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	276,670	C→T	intergenic (‑8/‑260)	unknown ← / → unknown	Peptidoglycan lipid II flippase MurJ/SSU ribosomal protein S20p

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	276,670	0	C	T	98.2%	141.3 / ‑4.6	54	intergenic (‑8/‑260)	unknown/unknown	Peptidoglycan lipid II flippase MurJ/SSU ribosomal protein S20p
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (25/28);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

ATAAATCGCCGGATTTCTGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATCACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTGTCACTATTAACTAACAGAAATTTT  >  1/276425‑276916                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       aTAAATCGCCGGATTTCTGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTActc                                                                                                                                                                                                                                                    <  1:656887/250‑1 (MQ=255) aTAAATCGCCGGATTTCTGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTActc                                                                                                                                                                                                                                                    <  1:253420/250‑1 (MQ=255)     aTCGCCGGATTTCTGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAAc                                                                                                                                                                                                                                               <  1:16657/251‑1 (MQ=255)                 cTGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAACTTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTg                                                                                                                                                                                                                                    >  1:427859/1‑250 (MQ=255)                  tGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGcc                                                                                                                                                                                                                                  <  1:485605/251‑1 (MQ=255)                                cacaGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGaaa                                                                                                                                                                                                                    <  1:113511/251‑1 (MQ=255)                                      cACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGaaa                                                                                                                                                                                                                    >  1:311171/1‑245 (MQ=255)                                      cACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGt                                                                                                                                                                                                              >  1:194720/1‑251 (MQ=255)                                            ggCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTa                                                                                                                                                                                                         >  1:686966/1‑250 (MQ=255)                                                    aaaaCGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTatat                                                                                                                                                                                                 <  1:157541/250‑1 (MQ=255)                                                      aaCGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCt                                                                                                                                                                                              <  1:280201/251‑1 (MQ=255)                                                          ccccTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAgg                                                                                                                                                                                          <  1:340708/251‑1 (MQ=255)                                                          ccccTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAgg                                                                                                                                                                                          <  1:369749/251‑1 (MQ=255)                                                           cccTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAggg                                                                                                                                                                                         <  1:311380/251‑1 (MQ=255)                                                              ttCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTaa                                                                                                                                                                                            >  1:677306/1‑245 (MQ=255)                                                                                  ggaAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTc                                                                                                                                                                                               >  1:135980/1‑222 (MQ=255)                                                                                            aaTGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTc                                                                                                                                                        <  1:339738/251‑1 (MQ=255)                                                                                               gCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAAtct                                                                                                                                                                       <  1:896228/233‑1 (MQ=255)                                                                                                cGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAAc                                                                                                                                                                                                                                               >  1:559926/1‑160 (MQ=255)                                                                                                           tAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTaa                                                                                                                                                                                            <  1:8998/200‑1 (MQ=255)                                                                                                           tAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGaa                                                                                                                                         >  1:20425/1‑251 (MQ=255)                                                                                                             aaaaTACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAg                                                                                                                                       <  1:31770/251‑1 (MQ=255)                                                                                                             aaaaTACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAg                                                                                                                                       <  1:37590/251‑1 (MQ=255)                                                                                                                            cAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAatat                                                                                                                         <  1:829990/250‑1 (MQ=255)                                                                                                                            cAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAatat                                                                                                                         <  1:148857/250‑1 (MQ=255)                                                                                                                               cgcCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTaa                                                                                                                                                                                            >  1:39508/1‑180 (MQ=255)                                                                                                                                 ccTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATCACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGa                                                                                                                    <  1:97673/250‑1 (MQ=255)                                                                                                                                                     aaTCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTa                                                                                                                              >  1:229602/1‑220 (MQ=255)                                                                                                                                                       tCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCa                                                                                             >  1:677089/1‑251 (MQ=255)                                                                                                                                                          ccacATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCg                                                                                          >  1:45444/1‑251 (MQ=255)                                                                                                                                                          ccacATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCg                                                                                          >  1:819765/1‑251 (MQ=255)                                                                                                                                                          ccacATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCg                                                                                          <  1:66894/251‑1 (MQ=255)                                                                                                                                                                         aGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAg                                                                                                                                                                                                                       >  1:247048/1‑111 (MQ=255)                                                                                                                                                                          gCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGtc                                                                                                >  1:201164/1‑229 (MQ=255)                                                                                                                                                                                      cgAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACaaa                                                              <  1:94132/251‑1 (MQ=255)                                                                                                                                                                                          aTTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTGAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGta                                                          >  1:630702/1‑251 (MQ=255)                                                                                                                                                                                          aTTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGta                                                          >  1:213893/1‑251 (MQ=255)                                                                                                                                                                                                 tCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAgg                                                                                                                                                                                          <  1:467016/116‑1 (MQ=255)                                                                                                                                                                                                       cgcTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGc                                              <  1:533190/250‑1 (MQ=255)                                                                                                                                                                                                                  aTGCCGGATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGcc                                   <  1:869037/250‑1 (MQ=255)                                                                                                                                                                                                                       ggATTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAAt                                                                                                                                                                                                                                       >  1:746892/1‑49 (MQ=255)                                                                                                                                                                                                                         aTTTTAATAATCTTTTGCTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTAttt                            <  1:515329/250‑1 (MQ=255)                                                                                                                                                                                                                         aTTTTAATAATCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTAttt                            <  1:519385/250‑1 (MQ=255)                                                                                                                                                                                                                             taataaTCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGaaa                                                                                                                                                                                                                    <  1:691695/62‑1 (MQ=255)                                                                                                                                                                                                                             taataaTCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCg                                                                                          >  1:458763/1‑184 (MQ=255)                                                                                                                                                                                                                                  atCTTTTACTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTGTCACTAtt                  <  1:707070/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTGTCACTATtaactaa            >  1:574606/1‑250 (MQ=255)                                                                                                                                                                                                                                         aCTCAAAACGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTGTCACTATtaact              >  1:504444/1‑248 (MQ=255)                                                                                                                                                                                                                                                aCGGATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGtt                                                                                                                               <  1:283333/128‑1 (MQ=255)                                                                                                                                                                                                                                                  ggATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAgg                                                                                                                                                                                          <  1:127018/67‑1 (MQ=255)                                                                                                                                                                                                                                                  ggATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTc                                                                                                                                                        >  1:843702/1‑101 (MQ=255)                                                                                                                                                                                                                                                  ggATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTa                                       >  1:569400/1‑214 (MQ=255)                                                                                                                                                                                                                                                  ggATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTg                          >  1:292879/1‑227 (MQ=255)                                                                                                                                                                                                                                                  ggATTACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTGTCACTATTAACTAACAGAAAtttt  >  1:484183/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ATAAATCGCCGGATTTCTGATATTCCGCCAACACAGGCACGAAGGCTTGTGAAAACGCCCCTTCGGCAAACAAGCGGCGGAGGAAATTCGGAATGCGATTGGCAAATAAAAATACGTCGGCGGCAGCGCCTGCGCCGATTAAATTGGCAATCACCACATCACGCACCAAGCCTAATACACGCGAAATTAAGGTCATGGCGCTGACGACGATGCCGGATTTTAATAATCTTTTACTCAAAACGGATCACTCTGAACGATAAAATTTGCCTAATTGTATAGAAATTTCGTTTTTACGCTATATTCTAAGGGAAAAAACTGTTAAAATCTCGCGCACATCGTTTCGTGAGCCGGTTTGAAAGTGCTTGTTAAATATAACGATGCTTTCTTTCGGTTGTGTCTCACCGAAATCAACACAAAATTCCCGTTGACAAAAGTATAAGAAATCGGCATATTTACGCCCCTATTTTGTCACTATTAACTAACAGAAATTTT  >  1/276425‑276916

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