Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,532	T→C	L119L (TTA→CTA)	creD →	putative inner membrane protein CreD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,532	0	T	C	95.7%	132.8 / ‑2.1	48	L119L (TTA→CTA)	creD	putative inner membrane protein CreD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (23/22);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

ACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCT  >  NC_000913/4638287‑4638767                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aCCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTTGTTACCTGAGTCATTGATGGTTGATGGCAATAAGAACGTGGAAGAACGAAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTaaa                                                                                                                                                                                                                                         >  1:357375/1‑250 (MQ=255)  ccGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAg                                                                                                                                                                                                                                       <  1:184093/251‑1 (MQ=255)              ggAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTGATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGGTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATg                                                                                                                                                                                                                           <  1:152121/251‑1 (MQ=255)                  gATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTc                                                                                                                                                                                                                       >  1:342239/1‑251 (MQ=255)                            gCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGc                                                                                                                                                                                                             <  1:373069/251‑1 (MQ=255)                                       aGCGCGCCGCAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTCTCCGGTGCAGGAAGAGGATAAAACAGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATAGTTGATGGCAATCAGAACGTGGCAGAACGCAAGATAGGGGTTTATACCGGTCAGGTCTTGCACAGTGATTTAACTCAAAAAGCCGATTTCGCTGTTGCGCGTCCTAGCgcaccccacgc                                                                                                                                                                                                  >  2:378572/1‑241 (MQ=255)                                               gCAAAAACTAGTTAGGCCGCTTATCGCTAGTCCTGTGAACGCGCTTTATCCGGTGCAGTAAGAGGATAACACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGCTGATGGCAATCAGAACGTGGAGGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTTACGCTCAAAGCCGTGTACGATGTTTCGCGTCTTAGCGAGATCAACGCGCCacat                                                                                                                                                                                           >  2:242792/1‑247 (MQ=255)                                                cAAAAACTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAatat                                                                                                                                                                                         >  1:56709/1‑251 (MQ=255)                                                 aaaaaCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATc                                                                                                                                                                                        >  1:152849/1‑251 (MQ=255)                                                  aaaaCTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATc                                                                                                                                                                                        >  1:19010/1‑250 (MQ=255)                                                    aaCTAGTTGGGCCGCTTATCGCTATTCCTGTGACAGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCACTCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCAACGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCGCCGCGCCAAATATCatc                                                                                                                                                                                     >  2:266253/1‑249 (MQ=255)                                                      cTCGTTGGGCCGCTTATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCtt                                                                                                                                                                                   >  1:233656/1‑251 (MQ=255)                                                        cGTTGGGCAGCTTATCGCTATTCCTGTGTCCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAg                                                                                                                                                                                 <  1:307783/251‑1 (MQ=255)                                                            gggCCGCTTATCGCTCTTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGCAGTTTTATCCATTTTTGGTTACCTGAGTCATTGCTGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCa                                                                                                                                                                              >  1:284805/1‑250 (MQ=255)                                                                           aTTCCTGTGACCTAGCCTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATtgtg                                                                                                                                                                <  1:26692/249‑1 (MQ=255)                                                                           aTTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATtgtg                                                                                                                                                                <  2:152849/249‑1 (MQ=255)                                                                               cTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTCAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGc                                                                                                                                                          >  1:16417/1‑251 (MQ=255)                                                                                 gtgACCGAGCTTTATTCGGTGCAGGAAGAGGATAAAACCTTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTGGATGGCAATCAGAACGTGGAAGAACGCAAGATATGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCTCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCg                                                                                                                                                         <  2:486487/250‑1 (MQ=255)                                                                                         gCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCAAAATATCATCTTAGGCAAGCCCTTTATTGTGCTTAGCGTCGGGGATg                                                                                                                                                >  2:418483/1‑251 (MQ=255)                                                                                              ataCGGTGCAGGAAGAGGATAAAACCGGGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGACGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGcgcgg                                                                                                                                           >  1:149178/1‑250 (MQ=255)                                                                                                       aGGAAGAGGATAAAACCGTGGAGCGGATACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATtggt                                                                                                                                  <  2:348797/251‑1 (MQ=255)                                                                                                               gATAAAACCGGGGAGCGGAAACGAAGTTTTATCCATCTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACATTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCTGTTATTGTGATGAGCGTCGGGGATGCGCGTTGTATTGGTGTGGTGaa                                                                                                                          <  2:312091/251‑1 (MQ=255)                                                                                                                               ggAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaac                                                                                                          >  1:396787/1‑251 (MQ=255)                                                                                                                                 aaaCGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAAGCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacgg                                                                                                        <  2:19010/251‑1 (MQ=255)                                                                                                                                     gAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTaacggaacg                                                                                                    >  1:177072/1‑251 (MQ=255)                                                                                                                                       aGTTTTATCCAGTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGGGATTAGCGTCGGGGATGCGCTTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGgc                                                                                                  <  2:173303/251‑1 (MQ=255)                                                                                                                                             aTCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGTTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGGGTGTTGAAAGCGCCTGAAGTTAACGGAACGTCGCTGAc                                                                                            <  1:331233/251‑1 (MQ=255)                                                                                                                                             aTCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGAc                                                                                            >  1:222853/1‑251 (MQ=255)                                                                                                                                                catgtttGGTTACCTGAGTCACTGATGGTTGATGGCAATCAGAACGTGGAAGAAAGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCAt                                                                                         <  2:222853/247‑1 (MQ=255)                                                                                                                                                       ggTTACCTGAGTCATTGATGGTTGATGGCAAGCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTCAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAAccc                                                                                  <  2:16417/251‑1 (MQ=255)                                                                                                                                                          tACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGcc                                                                                                                                                                          >  2:193986/1‑160 (MQ=255)                                                                                                                                                          tACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGcc                                                                                                                                                                          <  1:193986/160‑1 (MQ=255)                                                                                                                                                                   cATTGATGGTTGAGGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTGTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGtt                                                                       <  2:56709/250‑1 (MQ=255)                                                                                                                                                                   cATTGATGGGTGCTGGCAAACAGAACGTGGAAGAAAGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTCACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCACCTATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGACGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGGAGTTAACGGAACGGCGCTGCCCATTGAACCCGGATCCGGGtt                                                                       >  2:251357/1‑250 (MQ=255)                                                                                                                                                                     ttgatggttgatggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGtt                                                                       >  1:105091/1‑248 (MQ=255)                                                                                                                                                                          ggttgatgtCAATAAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGAGTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAg                                                               <  2:316273/251‑1 (MQ=255)                                                                                                                                                                             tgatggCAATCAGAACGTGGAAGAACGCAAGATATGGGTTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGCGGGGAAAGCGCCTGAAGTTAACGGACCGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCCAGGCg                                                            >  2:24101/1‑251 (MQ=255)                                                                                                                                                                               atggCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCggg                                                          >  1:45351/1‑251 (MQ=255)                                                                                                                                                                                       taGAACGTGGCCGAACGTAAGATAGGGATTTAACCCGGTCTGGTCTGTCACAGTGATTTAACCATAAAAGCCGATTTCGATGTTTCGCGTCTTAGCTAGTTCAACGCGCCAAATATCCTCTTAGGCAAGCCATTTATTGTGAATAGCGCCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGTTGCTGACCATTGAACCCGGCCCCGGGTTAGAGCACGGCGGGccggtct                                                   <  2:130830/249‑7 (MQ=255)                                                                                                                                                                                               ggAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATcc                                          >  1:445126/1‑251 (MQ=255)                                                                                                                                                                                                 aagaaCGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGGGGGCATATcccc                                        >  2:324496/1‑250 (MQ=255)                                                                                                                                                                                                    aaCGCAAGATAGGGATTTATACCGGGCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTGCGCGTCGTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTAGTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTa                                     <  2:488736/251‑1 (MQ=255)                                                                                                                                                                                                     aCGCAAGATAGGGATTTATTCCGGTCAGGTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCTCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTGTAc                                    <  2:300284/251‑1 (MQ=255)                                                                                                                                                                                                                     tataCCGGTTGGGTCTGGCACAGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGTGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGc                     <  2:345690/250‑1 (MQ=255)                                                                                                                                                                                                                                ggTCTGTCACAGTGTTTTAACCCTAAAAGCCGATTTCGTTGTTTCTCGTCTTATCGAGCTCAACGCGCCAAATATCATCTTTGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCACTTTACCCGAAGGAGACAGGCGGACGCAGaccc         <  2:283673/251‑4 (MQ=255)                                                                                                                                                                                                                                ggTCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAAcc         <  2:165483/251‑1 (MQ=255)                                                                                                                                                                                                                                  tCTGGCACAGTGATTTAACGCTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGAAGGGGGAGCATATCCCTTTACCCGAAAGGGACTTGCGGagagagaacc         >  2:83155/1‑240 (MQ=255)                                                                                                                                                                                                                                       cacaGTGATTTAACGCTAAAATCCGATTTCGATGTTTCGCGTCTTAGCGAGCTCAACGCGCCAAATATCATCTTAGGCAAGCCCTTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCCGAAGGGGACTGGCGGAAGCAGAACCTGAAGCt  <  1:187730/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            ACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGGGCCGCAAAAACTCGTTGGGCCGCTCATCGCTATTCCTGTGACCGAGCTTTATACGGTGCAGGAAGAGGATAAAACCGTGGAGCGGAAACGAAGTTTTATCCATTTTTGGTTACCTGAGTCATTGATGGTTGATGGCAATCAGAACGTGGAAGAACGCAAGATAGGGATTTATACCGGTCAGGTCTGGCACAGTGATTTAACGTTAAAAGCCGATTTCGATGTTTCGCGTCTTAGCGAACTCAACGCGCCAAATATCACCTTAGGCAAGCCATTTATTGTGATTAGCGTCGGGGATGCGCGTGGTATTGGTGTGGTGAAAGCGCCTGAAGTTAACGGAACGGCGCTGACCATTGAACCCGGCACCGGGTTAGAGCAAGGCGGGCAGGGCGTGCATATCCCTTTACCTGAAGGGGACTGGCGGAAGCAGAACCTGAAGCT  >  NC_000913/4638287‑4638767

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