Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,696,272	C→A	intergenic (‑201/+190)	uidA ← / ← uidR	beta‑glucuronidase/DNA‑binding transcriptional repressor UidR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,696,272	0	C	A	95.2%	117.7 / ‑0.3	42	intergenic (‑201/+190)	uidA/uidR	beta‑glucuronidase/DNA‑binding transcriptional repressor UidR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (25/15);  total (26/16)
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.54e-01

GAGTTTTTTGATTTCACGGGTTGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGACAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTG  >  NC_000913/1696027‑1696506                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            gAGTTTTTTGATTTCACGGGTTGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCACTACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGa                                                                                                                                                                                                                                          >  1:344828/1‑249 (MQ=255)   gTTTTTTGATTTCACGGGTTGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGcc                                                                                                                                                                                                                                      >  1:50357/1‑251 (MQ=255)          gATTTCACGGGTTGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAAt                                                                                                                                                                                                                               >  1:131792/1‑251 (MQ=255)                 cgggatagggTTTGTACGGGAAGTAACATAAGGGACTCCTAATTAAGAGAATAGTAGTGGTCAACCTGTAATCTGATCAGTCAAACAATTCAGAGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAACTAAAAAGAAATAAGACCGCGATAGAAAGGCCAAAAAATGAACTACAC‑AAAGCGTGA‑TATA‑GATCGCATTAAACGTTCACAAAGCGTTGAAAGAAGCAATAAAAAAAGAA‑Ac                                                                                                                                                                                                                       >  2:533656/8‑251 (MQ=255)                   ggTTGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCg                                                                                                                                                                                                                     >  1:164194/1‑251 (MQ=255)                      tGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGACAGAAGCCATAAAATAAGCA‑ACCGaa                                                                                                                                                                                                                   <  1:70211/250‑1 (MQ=255)                               tACAGGACGTAACATAAGTGACTCCGCATTAAGATAATAAGACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACtttcgttc                                                                                                                                                                                                          <  2:228297/250‑1 (MQ=255)                                                     tcctcATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATTGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCAATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAAct                                                                                                                                                                                   >  1:167115/1‑251 (MQ=255)                                                                              cAACCTTTAATCTGAATAGTTAACCAATTCTGTGTTTCACACAGTTATCTGTGAGAAACCAATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGTGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGaa                                                                                                                                                           <  2:290821/250‑1 (MQ=255)                                                                                  cTTTAATCTGACTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCGATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAACAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAACACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTACGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGacgattt                                                                                                                                                      >  1:528157/1‑245 (MQ=255)                                                                                      aaTCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTgg                                                                                                                                                   <  1:316702/250‑1 (MQ=255)                                                                                         cTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTgggggg                                                                                                                                               >  1:94773/1‑251 (MQ=255)                                                                                          tGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGTGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTgggggg                                                                                                                                               >  1:35379/1‑250 (MQ=255)                                                                                           gATTAGTTAACCAAGTCTGTGTTGCACACAGTTATCGGGGAGAAACCTATGGAAGAGAAATGGAAATAAAAAGAACTGTCACCGCTATCAAAAAACCAAAAACTGCACTACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCGTTAAACCAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACGCTCATGACATATGGCTACAGTGAATATTTTGGGGGGa                                                                                                                                              <  2:303176/250‑1 (MQ=255)                                                                                              tAGTTAACCAATTCTGTGGTGCACACAGTTATCTGTGAGAAACCCATTGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAAtg                                                                                                                                           <  2:102178/250‑2 (MQ=255)                                                                                                   aaCCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGa                                                                                                                                     <  1:70753/251‑1 (MQ=255)                                                                                                                  gCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATa                                                                                                                       >  1:307253/1‑250 (MQ=255)                                                                                                                    acacaGTTATCTGTGAGAAACCCATGGAAGAGAAATGGCAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACAGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGGGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAAtt                                                                                                                    >  2:98211/1‑251 (MQ=255)                                                                                                                     cacaGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAAttt                                                                                                                   >  1:36830/1‑251 (MQ=255)                                                                                                                            aTCTGAGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGCACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGcac                                                                                                             >  1:448893/1‑250 (MQ=255)                                                                                                                               ggtgAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATAAAAAAACCAAAAACTGCAATGCCA‑ACAGCGTGA‑TATA‑GAACGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCTAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACaga                                                                                                          <  2:420641/249‑1 (MQ=255)                                                                                                                                       aCCCATGGAAGATAAATGGAAATAAAAAGAACTGACACTGCGATCAAAAATCAAAAAACTGCAAGACCC‑GCAGCTGGA‑GATA‑GATCTCATTAATCTTTAAAACAGCTTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCTAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTGGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACg                                                                                                 <  2:306422/251‑1 (MQ=255)                                                                                                                                         ccagggAAGAAAAATGAAAATAAAAAGAACTAAAACCGCGATCAAAAAACAAAAAAC‑GCAAAACCCCACAGC‑TGA‑AAAAAGATCGCATAAACCTTTAAAAAAGCGTGAAAAGAAGCC‑TTAAATAAGCA‑ACCGAACGTTCGTTCTAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATGCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGata                                                                                              <  2:519888/247‑1 (MQ=255)                                                                                                                                          ccaggaaaggaAATGGAAATAAAAAGAACTGACAACGCTGTCAAAAAACCAAAAAATGCAATACCC‑ACAGCGTGA‑TATA‑AATCGCATTAATTGTTAAAACAGCGTTTAAAGAAGCCATAAAATAAGCA‑TCCAAACGCTCGTTCGAGCTTATTGGTAACGATAAACTCTCATGACATATGTCTTCAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGGTAATTTTTGCACAGATGAATTACGata                                                                                              <  2:336951/242‑1 (MQ=255)                                                                                                                                          cATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATCCCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAATCAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCATAGAGAAATTACGata                                                                                              <  1:523674/251‑1 (MQ=255)                                                                                                                                           aTGGAAGAGAACTGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGAGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGACAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTATCATGACATATGGCTACAGAGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTAAGatat                                                                                             >  2:321095/1‑251 (MQ=255)                                                                                                                                             ggctgccccagccgatgtaaaaGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGata                                                                                              <  1:399152/230‑1 (MQ=255)                                                                                                                                             ggAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACAAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTATCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATg                                                                                            >  1:198527/1‑250 (MQ=255)                                                                                                                                              gAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCGTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGaa                                                                                          <  2:512576/251‑1 (MQ=255)                                                                                                                                                  agaAATGGAAATAAAAAGAACTGACACCGCGATCCAAACACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGGAATTACGATATGAAAcgg                                                                                      >  2:330700/1‑250 (MQ=255)                                                                                                                                                            aTAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTTGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGtt                                                                            >  1:186567/1‑251 (MQ=255)                                                                                                                                                                         acacCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATAACCAGCccgc                                                               >  2:495753/1‑251 (MQ=255)                                                                                                                                                                                     aaaaaCCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAagag                                                     >  2:365896/1‑249 (MQ=255)                                                                                                                                                                                       aaaCCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAagagag                                                   >  2:490961/1‑249 (MQ=255)                                                                                                                                                                                         acaaaaaaCTGCAAAACCC‑AAAGC‑TGAATTTA‑GTTCGATTTAATCCTTAAAAACGCGTTGAAAAAAGACATAAAATAAGCA‑ACCGAACGTCCTTTCGAGCGTATGGGTAAATATAAACTCTCTTTACATATGGGTATAGTTAAGATATTGGGGGGAATTCTCCGAATCCGGGAAGGATAATTTTGGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCATCCCGCCGAAAAAGAGAGAGaa                                               <  2:398872/248‑1 (MQ=255)                                                                                                                                                                                           cAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATc                                             >  1:540937/1‑251 (MQ=255)                                                                                                                                                                                                     aaTACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGt                                   >  1:372125/1‑251 (MQ=255)                                                                                                                                                                                                     aaTACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGt                                   >  2:101102/1‑251 (MQ=255)                                                                                                                                                                                                                            gATCGCATTAATCTTTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATcg               >  1:330166/1‑251 (MQ=255)                                                                                                                                                                                                                                 cATTAATCTGTAAAACAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGATCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAAAGAATCAGGATGCGGTTAAGTTACCGCCAATCATCGCGCGt          >  2:344425/1‑251 (MQ=255)                                                                                                                                                                                                                                     aagatttAAAACAGCGTCGAAAAAAGCCATAAAATAAGGACACGGAA‑GTTCGTGCGAACGTAAGGGGAACGATAAACCCTCATGACATAGGTCGAAAGTGAATATTTTGGGGGGAAGTCTCCGAAGCCGGGAAGGATAATTTTTGCCCAGAGAAATTACGAGATGAACCGCGGGCACAGTTACCGCCAGCCCGCCGAAAAAGAGGGAGAACCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAAt       <  2:212763/246‑1 (MQ=255)                                                                                                                                                                                                                                         tttAAAAAAGCGTTGAAAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGCTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGCATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAAACATCGCGCGTAATCCCTg  >  1:39611/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            GAGTTTTTTGATTTCACGGGTTGGGGTTTCTACAGGACGTAACATAAGGGACTCCTCATTAAGATAATAATACTGGTCAACCTTTAATCTGATTAGTTAACCAATTCTGTGTTGCACACAGTTATCTGTGAGAAACCCATGGAAGAGAAATGGAAATAAAAAGAACTGACACCGCGATCAAAAAACCAAAAACTGCAATACCC‑ACAGCGTGA‑TATA‑GATCGCATTAATCTTTAAAACAGCGTTGACAGAAGCCATAAAATAAGCA‑ACCGAACGTTCGTTCGAGCGTATTGGGAACGATAAACTCTCATGACATATGGCTACAGTGAATATTTTGGGGGGAATTCTCCGAAGCCGGGAAGGATAATTTTTGCACAGAGAAATTACGATATGAAACGCGGGCACAGTTATCACCAGCCCGCCGAAAAAGAGAGAGAATCAGGATGCGGTTAAGATACCGCCAATCATCGCGCGTAATCCCTG  >  NC_000913/1696027‑1696506

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