Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,738,631	T→C	L263P (CTG→CCG)	purR →	DNA‑binding transcriptional repressor PurR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,738,631	0	T	C	95.7%	135.2 / ‑3.7	47	L263P (CTG→CCG)	purR	DNA‑binding transcriptional repressor PurR
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (22/23);  total (23/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 = 8.86e-01

CATCCCCGGCCCGCTGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTGGCTGACGGCCCGTTCCGCGAC‑TATCGTCGTTAATCACCCG  >  NC_000913/1738392‑1738876                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       cATCCCCGGCCCGCTGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGc                                                                                                                                                                                                                                             >  1:107955/1‑251 (MQ=255)                  aaCGTAACACCGGAGCAGGCCGCCTTGCCGGGTTTAGGAAGGCGATGGAAGAAGCGATGACCAAGGGGCCGGAAAGCTGGA‑TTGGGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGGC‑GCAGCCGCAGCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTGGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTATCGCTGATc                                                                                                                                                                                                                            <  2:227794/251‑1 (MQ=255)                          acCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGCTCAAGGTCCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGATGGCGATATCATGGCAATGAGCGCACTTTGTGCTGCTGATGAACTGGGCCCGCGCATCCCCCCGGATGTTTCgaggttcagttctgg                                                                                                                                                                                                                    >  2:283133/1‑237 (MQ=255)                                gcAGGCCGCCTTGCCGGTTTTAGGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGt                                                                                                                                                                                                              <  2:235164/251‑1 (MQ=255)                                   ggCCTCCTTGCCGGTTTTATGAAGTCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTgcg                                                                                                                                                                                                           <  1:420083/251‑1 (MQ=255)                                   ggCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGAGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTgcg                                                                                                                                                                                                           >  1:461227/1‑251 (MQ=255)                                       cgctcttccgaacgaggAAGGCGAGGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGACCCTGAACCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACGGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAc                                                                                                                                                                                                       <  2:207826/235‑1 (MQ=255)                                                gTTTTATTAAGGCGATGGAAGAAGCGATGAGCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGTATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTAt                                                                                                                                                                                              <  2:305646/251‑1 (MQ=255)                                                     aTGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAAc                                                                                                                                                                                                       >  1:207826/1‑237 (MQ=255)                                                      tGAAGGCGATGGGAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTa                                                                                                                                                                                          >  1:142706/1‑249 (MQ=255)                                                            cgctGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGCTATCATGGCAATGTGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCACGCAGGATGTTTCGCTGATCGGTTATGATAACGTACGCAACGCGCGCTATTTTACGCCGgcg                                                                                                                                                                                  >  2:470101/4‑251 (MQ=255)                                                            cGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGGC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGgcg                                                                                                                                                                                  <  2:348179/251‑1 (MQ=255)                                                                          atgatCAAGGTGCCGGAAAGCTGGACTTG‑GCAGGGTGACTTTGAACCTGAATACGGTTATCGCGCCATGCTGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTTCTGCTGATGCAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGGTCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATcc                                                                                                                                                                     <  2:458356/250‑1 (MQ=255)                                                                                 aGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGcc                                                                                                                                                              >  1:507026/1‑250 (MQ=255)                                                                                       cGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAgct                                                                                                                                                       >  1:20008/1‑249 (MQ=255)                                                                                                     ggtgCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGaaa                                                                                                                                          <  1:371939/250‑1 (MQ=255)                                                                                                       tgCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGCCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCCAAAGATTCGCTGGGTGAAAc                                                                                                                                         >  2:3814/1‑250 (MQ=255)                                                                                                          cgtgttaCTTTTAACCTGAATCCGGTTTTCTCGCCATGCAGCAAATCCTGTCTGC‑GGCGCATCTCCCTACTGCCGTTTTCTGTGGTGGCGATATCATGACAATGGGAGCACTTTGTGTTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCg                                                                                                                                     <  2:367365/245‑1 (MQ=255)                                                                                                           gggTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCTCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGt                                                                                                                                    >  1:110469/1‑251 (MQ=255)                                                                                                                     gAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGGCGCATCGCCCTACTGCCGTCTTCTGGGGTGGCGATATCATGGCAATGGGCGCGCTTTGTGCTGCTGAAGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAATATTCGCTGGGTGAAACGGCGTTCAACATGCt                                                                                                                          <  2:514289/251‑1 (MQ=255)                                                                                                                     gAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCTGTCTTCGGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGc                                                                                                                           <  1:201625/250‑1 (MQ=255)                                                                                                                       aCCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGt                                                                                                                        <  1:458531/251‑1 (MQ=255)                                                                                                                                  gTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgtatc                                                                                                             <  1:168241/251‑1 (MQ=255)                                                                                                                                    tATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGc                                                                                                                                      >  2:74167/1‑224 (MQ=255)                                                                                                                                    tATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGc                                                                                                                                      <  1:74167/224‑1 (MQ=255)                                                                                                                                     aTCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCCTGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGCTAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAACAATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTc                                                                                                          >  2:210040/1‑251 (MQ=255)                                                                                                                                      tCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTc                                                                                                          >  1:538122/1‑250 (MQ=255)                                                                                                                                       cgcgCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTca                                                                                                         >  2:35148/1‑250 (MQ=255)                                                                                                                                          gcCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGCAACAGCGTTCAACATGCTGTTGGatcgta                                                                                                               >  2:137151/1‑241 (MQ=255)                                                                                                                                          gcCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGatcgta                                                                                                               <  1:137151/241‑1 (MQ=255)                                                                                                                                                                  gccgccTCGCCCTACAGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGCTGAAATGGGCCCGCGCGTCCCGCAGGATGCTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTTCGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGTTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAATCGTGAAGAACCGCAGTCTATTg                                                                              <  2:142706/251‑1 (MQ=255)                                                                                                                                                                       aTCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACAATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTg                                                                         >  1:196086/1‑251 (MQ=255)                                                                                                                                                                             cTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGACCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCcg                                                                   >  1:559778/1‑251 (MQ=255)                                                                                                                                                                                  gCCGTCTTCTGGGGGGGCGATATCATGGCAATGGGCGCCCTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTGACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCtt                                                              <  2:538122/251‑1 (MQ=255)                                                                                                                                                                                    cGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAACGATTCGCTGGGTGAAACAGAGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCttg                                                             >  2:372660/1‑250 (MQ=255)                                                                                                                                                                                    cGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCttga                                                            <  1:3814/251‑1 (MQ=255)                                                                                                                                                                                           tgtgGTGGCGATATAATGGCAATGGGCGCACTTTGTGCTACTGATGAGATGGACCTGCGCGTCACGCCGGATGTTTCGCTGATCGGTTATGATAACTAGCGCAACGCGCGCTATTTTCCGCCGGCGCTGACCACGATCCATCAGCCAAACGCATCGCTGGGTGAAACAGCGCTCAACAGGCTGTTGTATCATAACGTCAACAAACGGCAAGGACCGAAGCCTAATGAAGGTCATCcgcgattactggaacg                                                     >  2:296550/1‑239 (MQ=255)                                                                                                                                                                                                    gATATCATGGCAATGGGAGCACTTGGTGCTGATGATGAAAAGGACCCGCGCGACCCGCAGGATGTTTCGATGATCGGTTATGATAACGAGCGCAACGCGCGCGATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGACACAGCGTTCAACATGCTGTTGGAACGTATCGTCAACAAACGTGAAGAAACGCAGTCAATTGAAGTGCATAGGAGCTTGATTGAAAGCCGCTGCGt                                            >  2:561405/1‑251 (MQ=255)                                                                                                                                                                                                       atCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGCAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTGGc                                         >  1:346426/1‑251 (MQ=255)                                                                                                                                                                                                       atCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTgg                                          <  2:107955/250‑1 (MQ=255)                                                                                                                                                                                                       atCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTgg                                          <  1:326156/250‑1 (MQ=255)                                                                                                                                                                                                                 tGGGCGCACTTTGTGCTGCTGATGAAATGGGCCCGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTGGCTGACGGCCCg                               >  1:246048/1‑251 (MQ=255)                                                                                                                                                                                                                    gcgcACTTTGTGCTGCTGATGAAAGGGGCCCGCGCGGCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATGTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTGGCTGACGGCCCGTTc                            <  2:527461/251‑1 (MQ=255)                                                                                                                                                                                                                     cgcACTGTGTGCTGCTGATGAAATGGGCCTTCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTACGTGGCTTACGGCCCGTTcc                           <  1:296550/251‑1 (MQ=255)                                                                                                                                                                                                                                  gctgATGAAATGGGCACGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTGGCTGACGGCCCGTTCCGCGAC‑TAtcgtcg             <  1:196626/251‑1 (MQ=255)                                                                                                                                                                                                                                  gctgATGAAAGTGGCCCGCGCGTCCCGGAGGATGGTGCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTAGTTTACGCCGGCGCGGACCACGATCCATCAGCCAAAAGAGTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCTTTTGATTGAACGCCGCTCCGTGGCTGACGGCCCGGTCCTCGAC‑TAtcgtc              <  2:246048/250‑1 (MQ=255)                                                                                                                                                                                                                                             gggCCCGCGCGTCACGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACCAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGCTTGAACGCCGCTCCGTGGCTGACGGCCCGTTCCGC‑ACATTACGTCGTTAATCACCCg  >  2:520360/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       CATCCCCGGCCCGCTGGAACGTAACACCGGCGCAGGCCGCCTTGCCGGTTTTATGAAGGCGATGGAAGAAGCGATGATCAAGGTGCCGGAAAGCTGGA‑TTGTGCAGGGTGACTTTGAACCTGAATCCGGTTATCGCGCCATGCAGCAAATCCTGTC‑GCAGCCGCATCGCCCTACTGCCGTCTTCTGTGGTGGCGATATCATGGCAATGGGCGCACTTTGTGCTGCTGATGAAATGGGCCTGCGCGTCCCGCAGGATGTTTCGCTGATCGGTTATGATAACGTGCGCAACGCGCGCTATTTTACGCCGGCGCTGACCACGATCCATCAGCCAAAAGATTCGCTGGGTGAAACAGCGTTCAACATGCTGTTGGATCGTATCGTCAACAAACGTGAAGAACCGCAGTCTATTGAAGTGCATCCGCGCTTGATTGAACGCCGCTCCGTGGCTGACGGCCCGTTCCGCGAC‑TATCGTCGTTAATCACCCG  >  NC_000913/1738392‑1738876

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