Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,608,681	T→C	11.7%	intergenic (+33/+69)	zntA → / ← tusA	zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,608,681	0	T	C	11.7%	154.0 / 2.3	58	intergenic (+33/+69)	zntA/tusA	zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (40/10);  new base C (2/6);  total (42/16)
Fisher's exact test for biased strand distribution p-value = 3.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCG  >  NC_000913/3608394‑3608957                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                         acaGACGTGGCGATGGAAACCGCCCACGCAGCATTAACCAATAACCACCTGCGCTGCATGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTCCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGt                                                                                                                                                                                                                                                                                          <  1:556036/288‑1 (MQ=255)  caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGGGCAGGAGATAAGGCAAACC‑GACCGCAACATTGAGCGG‑GATCGGT‑c                                                                                                                                                                                                                                                                                        >  2:232573/1‑288 (MQ=255)  caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGAAGGGGATAAGGCAAACC‑GATCGAAACATTGAGCGC‑GATCGGT‑c                                                                                                                                                                                                                                                                                        >  2:546936/1‑288 (MQ=255)  caGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAAGGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑c                                                                                                                                                                                                                                                                                        >  1:673585/1‑288 (MQ=255)      cGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑cc                                                                                                                                                                                                                                                                                       <  1:239238/285‑1 (MQ=255)       gTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCGGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCg                                                                                                                                                                                                                                                                                  <  2:663529/289‑1 (MQ=255)          gcgcTGGAAACCGCCGACGCGGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGcc                                                                                                                                                                                                                                                                                <  1:1061340/288‑1 (MQ=255)                     cgcctACGCAGCATTAACCCACAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCCAGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGAGCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATGGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGga                                                                                                                                                                                                                                                                      <  1:423918/287‑1 (MQ=255)                         gACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGt                                                                                                                                                                                                                                                               <  1:546936/289‑1 (MQ=255)                                          tAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAg                                                                                                                                                                                                                                            >  1:1193886/1‑289 (MQ=255)                                                ccTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTGAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGT‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GCGGGTTA‑GGGTGAGGGG‑AAAAGGCCGCa                                                                                                                                                                                                                                      >  1:970056/1‑289 (MQ=255)                                                 cTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCa                                                                                                                                                                                                                                      >  1:916865/1‑288 (MQ=255)                                                 cTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCa                                                                                                                                                                                                                                      >  1:81087/1‑288 (MQ=255)                                                   gcgcGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGGTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTTGCACCTCTGGGGG‑GAGGGCTA‑GGGTGAGGGG‑AAAAGGCGGCATCg                                                                                                                                                                                                                                   >  1:91470/1‑289 (MQ=255)                                                                   aTGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GCGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                                                                                                                   >  1:705853/1‑289 (MQ=255)                                                                   aTGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                                                                                                                   >  1:231801/1‑289 (MQ=255)                                                                                         cTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATAAGCCCCTATCAACCGCCTTTAGGAATCaa                                                                                                                                                                                             >  1:1185765/1‑289 (MQ=255)                                                                                           cACGCAACTAGCCGCCAGAACAGCACTTTGGCGTTGGGGGTGAAGGGGTTTGCCCCCGACCCGAGGCTTTTGGGGA‑GAGCGGGTTGGGGGTGGGAGGGG‑TTGCGAGGTGGGGGGGGGGCGGGTCTGGGGAGGGCGAAGGTTTTAAGATGTTTGGGAAGGAGAAGAGGAAAACG‑GCTCGAAAAATGGGGGGC‑GATCGGC‑CCCCCCGCCCCTCGGGGGA‑GGGGTTTG‑GGGTGGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                                                                                                                           <  2:969746/289‑1 (MQ=255)                                                                                           cACCCCAATATCCGCAGAAACATCACTGTTGCGCGGGGGGGTAAGGGGTTCTCCCTGTTACCCACGGT‑TTGGGGGAGGGCGGGGTGTTGCGGGGAAGGCGTGGC‑AGATACGGGGGCGCGGGTGCGGGAGACACGGATTCGGTAAAGTTTGTTGC‑CAGGAAGAAAGGAAAACCGGTTCAAAAATTTGGGGGC‑GCCCGGC‑CCCCCCCCCCCTCGGGGAA‑GGGGGTAG‑GGGGGGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                                                                                                                                           <  2:1185765/289‑1 (MQ=255)                                                                                              gCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                                                                                                                                            >  1:301250/1‑269 (MQ=255)                                                                                                             aaCATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGga‑ga                                                                                                                                                                                                                                                                   <  1:574429/201‑1 (MQ=255)                                                                                                                       ttGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGCTGACCGGGTTGTGGCTGGCGGTGC‑TGGC‑AGATTCGGGGGCGACGGTGCTGGTGACAGCGGATGCGTTAAGATTGTTGCGCGGGAGATTAGGCAGACC‑GAACGGCAAATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGGTC‑GGGt                                                                                                                                                                                                                                                        >  2:574429/1‑201 (MQ=255)                                                                                                                               gggCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCGAATAACGATAAGGCAGGCAATCCGTCTCTTTAGCAACCa                                                                                                                                                       >  1:712971/1‑289 (MQ=255)                                                                                                                                               tCCTGGGCACCGCGCTGTTGGGTAGGACGGGGTTGGGGCTGCGAGTGT‑TGGC‑AGATAGGGGGGGGGCGGTGTGGGTGACAGCAGAGGAGTTAAGATTGTTGCGCGGGAGCTAAGGCAAACC‑GATCGAAACATTGGGCGC‑GATCGGT‑CCCCTCCCCCCTCGGGGGG‑GGGGGTTAGGGGGGAGGTG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                                                                                                                        <  2:370349/289‑1 (MQ=255)                                                                                                                                                        ccacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGCGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCACCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAAGCATCAAATAACGATAAGGCAGTCCATCCGTCCCTTTAGCAACCAGTTCGTGTTCCCTAAAGGTACaaac                                                                                                                              >  2:946913/1‑288 (MQ=255)                                                                                                                                                                      aTGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCg                                                                                                                                                                                                                                                                                  <  1:316182/130‑1 (MQ=255)                                                                                                                                                                      aTGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCg                                                                                                                                                                                                                                                                                  >  2:316182/1‑130 (MQ=255)                                                                                                                                                                      aTGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                                                                                                                >  1:1195043/1‑289 (MQ=255)                                                                                                                                                                        gACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTACCAACCGCCTTCACGAATCAAATAACGATAAGACAGTCCACCCGTCTCATTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGAgg                                                                                                              >  1:1162057/1‑289 (MQ=255)                                                                                                                                                                        gACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCATCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg                                                                                                              >  1:284912/1‑289 (MQ=255)                                                                                                                                                                        gACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATAGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAACCAGCCCCTATCAACCGCCCTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaactgtatcttatacacatccccgagc                                                                                                               >  1:98296/1‑263 (MQ=255)                                                                                                                                                                          ccGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                                                                                                          >  1:958764/1‑211 (MQ=255)                                                                                                                                                                                  tggctggcAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg                                                                                                    >  1:448495/1‑289 (MQ=255)                                                                                                                                                                                        gcAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAAAGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGC‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTAACACCCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCCAGTTCGTGTTCCATAACGGTACAAAACCCAGCAATCTCGCGGGTAGTGGCCGGAtcg                                                                                               >  2:969556/1‑288 (MQ=255)                                                                                                                                                                                            tGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGc                                                                                          >  1:551887/1‑289 (MQ=255)                                                                                                                                                                                                        tACGGGGGCGACGGTGCTGGTTACAGCGAATGCGTTAAGATTGTTGAGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGATA‑GGGTGAGGGG‑AAAAGGCGGCATCGAATCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGTTAAGGAAGTCCTTCCGTCTCTTTATCAATCAGTTCGTGTCCCATCAACGTACAAGACCCAGTCATATCCCGTGAGTTTGCGGGATCATCGGCGATAGTCaac                                                                                >  2:431836/1‑287 (MQ=255)                                                                                                                                                                                                                    ggtgcgggtgagaggaaactgttagaGATTGTGGGGCAGGAGAAAAGGAAACCG‑GATCGAAACTGTGGGCGC‑GAGGGGC‑CCCCCCCCCCCGCGGGGGA‑GAGGGTGA‑GTGTGGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGcc                                                                     <  2:551887/263‑1 (MQ=255)                                                                                                                                                                                                                         tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATcctgattatcgccgacgatcc                                                               >  1:825409/1‑269 (MQ=255)                                                                                                                                                                                                                                aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTATGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCAAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCAATAAAGGTACAAAACCAAGGAATATCGCGGGTAGTGGACGGATCGTCGGCGATAAACAGCAAGGGTTCGCAAGGCTGACTAtt                                                         >  2:29813/1‑288 (MQ=255)                                                                                                                                                                                                                                aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGC‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGACGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAACTAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTAGTGTTCCATAAAGGTACAAAACCCAGCAATATCGCGGGTAGTGGCCGGACTGTCGGCGATAACCAGCAACGTTTCGCCAGGCTGCATAtt                                                         >  2:865831/1‑288 (MQ=255)                                                                                                                                                                                                                                    aaggCGTTAAGTTTGGGGCGAAGGAGATAAGGCAAACC‑GA‑CGGAAACTTGGGGGCGGATCGCC‑CCCCCCGCCCTGTGGGGGA‑GGGGGTTA‑GGGTAGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTgcgc                                                     <  2:1156314/285‑1 (MQ=255)                                                                                                                                                                                                                                     aTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCAc                                                   >  1:1156314/1‑289 (MQ=255)                                                                                                                                                                                                                                                  gttgCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                                                                                                                    >  1:774697/1‑131 (MQ=255)                                                                                                                                                                                                                                                   ttgCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGTTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCCTATTGCGCACGGTTTTGCGCAc                                       >  1:615248/1‑287 (MQ=255)                                                                                                                                                                                                                                                   ttgCGAAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGA‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGCGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTATATTTAGCAAACAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATAGGGGGTAGAGGGCGGACCGTCGGCGATAATCCGGAACGTTTCGCCAGGCGGAATATTGGGAGCGGGTTTGCGAACca                                     >  2:1015131/1‑289 (MQ=255)                                                                                                                                                                                                                                                      cgcAGGGGAAAAGGCAACCT‑GACCGAAACAGTGAGCGC‑GCTCCGC‑CCCCTCTCCCCTCTGGGGA‑GGGGGTTA‑GGGTGGGGGC‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACcatc                                   <  1:865831/288‑1 (MQ=255)                                                                                                                                                                                                                                                         aGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCatcatc                                >  1:636233/1‑288 (MQ=255)                                                                                                                                                                                                                                                           acggaaaGGCAAACC‑CATGCTAACGTTGAAGTC‑GCTCGGT‑CCCC‑CGCCCCTCCGAGGAGGAGGGTAG‑GGGTGGGGGA‑AAAAGGCGGGATCGAAGCCAATCAGCCCATATCAACCGCCTTTACGAATCAAAGAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGGTCGTGTTCCATAAATGTACAAAACCCAGGAATATCGCGGGTAGTGTCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCAc                              <  2:290578/283‑1 (MQ=255)                                                                                                                                                                                                                                                                aaGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCACTCGCCCCTCTGGGGA‑GAGGGCTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGtt                         >  1:504009/1‑288 (MQ=255)                                                                                                                                                                                                                                                                     aaaCC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa                                                                                                                       >  1:36790/1‑189 (MQ=255)                                                                                                                                                                                                                                                                       aCC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGTG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCCACCGCCTTTACGAATCCAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAACACCCAGGAATATCGCGGGAAGTGGCCGGATCGTCGGAGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCCCCCGGTTTTGCGCACCATCCTCACCGGTTCCggtct                  >  2:806327/1‑285 (MQ=255)                                                                                                                                                                                                                                                                       aCC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTCCAGGGCAg                 >  1:325763/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                   cATTGAGCGC‑GATCGGT‑CCCCCCGCCCCTCTGGGGA‑GAGGGGTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCACGGCGAGTGGCCGGATCGTCGGAGATAATCAGCAACGTTTCGCCAGGCTGCATACTGCGCAGTGTTTTACGCCACAACATCACCGGTTCCGGGCAgcgcggcgcaag      >  2:939969/1‑281 (MQ=255)                                                                                                                                                                                                                                                                                    agttggggggggtCGGC‑CCCTTCGCCCCTGGGGGGA‑GAGGGTTA‑GGGTGGGGGA‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGCCAAg      <  2:705853/277‑1 (MQ=255)                                                                                                                                                                                                                                                                                    aatggaggc‑gGTCGGTCCCCCCCGCCCCGCGGGGGA‑G‑GGGTTT‑GGGGGCGGGGAAAAAGGCGGCATCGAAGCAAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGCCAAg      <  2:636233/282‑1 (MQ=255)                                                                                                                                                                                                                                                                                    aTTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCata                                                           >  1:564566/1‑235 (MQ=255)                                                                                                                                                                                                                                                                                        aGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCACTCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGTCCGGATCGTCGGCGATAATCAACCACGTTTCGCCAGGCCGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGCCAGCGCAGGCCAAGcgcg  >  2:70020/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                        aGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGCCAAGcgcg  >  1:452759/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                        aGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGCCAAGcgcg  >  1:107110/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                         ACAGACGTGGCGCTGGAAACCGCCGACGCAGCATTAACCCATAACCACCTGCGCGGCCTGGTGCAAATGATTGAACTGGCACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGC‑TGGC‑AGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGC‑GATCGGT‑CCCCTCGCCCCTCTGGGGA‑GAGGGTTA‑GGGTGAGGGG‑AAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACGTTTCGCCAGGCTGCATATTGCGCACGGTTTTGCGCACCATCATCACCGGTTCCGGGCAGCGCAGGCCAAGCGCG  >  NC_000913/3608394‑3608957

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