Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,623,911	C→T	A69A (GCG→GCA)	lplA ←	lipoate‑‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,623,911	0	C	T	97.5%	113.7 / ‑4.3	40	A69A (GCG→GCA)	lplA	lipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/25);  total (14/26)
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

GTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTC  >  NC_000913/4623671‑4624150                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 gtggAAGCCACGATCTTTGGTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCTTGCCAGGCGGa                                                                                                                                                                                                                                        <  1:175573/251‑1 (MQ=255)                tttGGTTTCGCGATAGGCCGATTCTGAGATGTTGCTGTCGCCTTCGGCGGTTTGCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTGAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCTCCGCTACTGCGCCGTTCCAGGCGTACGTTATCTTCTTcc                                                                                                                                                                                                                          <  2:170030/250‑1 (MQ=255)                  tggattCGCGATAGGCCGATCCTGAGACTTTGCGGTAGCCTTCGGCGGTTTTCAACACCAGATCGTTACGCCCGGACGCCTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCAAAATCGACGTGGGGGTGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGGATTGCCGGGGTCGTGGGACACCGCGCCGCCACCGCAACTGCGCCGTGCCAGGCGGGCGTTATCTTCTTCCATc                                                                                                                                                                                                                       >  2:193030/5‑251 (MQ=255)                  tcgTTTCGCGATAGGCCGTTCCTGATCCTTTGCGGTCGCCTTCGGCGGTTTTCCCCACCAGATCGGTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCTCAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCCCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCAt                                                                                                                                                                                                                        <  1:165983/248‑1 (MQ=255)                   ggTTTCGCGATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATcc                                                                                                                                                                                                                      <  1:155775/251‑1 (MQ=255)                        cgcgATAGGCCGATCCTGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAAACGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTAATGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCgg                                                                                                                                                                                                                 >  1:195214/1‑251 (MQ=255)                                  cgatcctgactcTTTGCGGTCGCCATCGGCTGTCTTCACCACCAGGTCTGTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCTTATCATTGAGCACAGTCCGCGTGGAGATGGTTTTATCGTACTCCAGCGTGTCAGCCATAAAGGTATAGCAGGTATTGCGTAGTTCGTGGCACCCCGCGCCTCCACCGCTACTCCGCCGTTCCAGGCATACGTTATCTTCTTCCATCCGCCGGGTATTACAct                                                                                                                                                                                                       <  1:117183/240‑1 (MQ=255)                                        tGAGACTTTGCGGTCGCCTTCGGCGGTTTTCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCCAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGGCGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTcc                                                                                                                                                                                                 >  2:33938/1‑251 (MQ=255)                                                           tcggcGGTTTTCACCACCCGATCGTTAAGCCCGGACGATTCGGCGCTGACGCCAAGCGCGTTCAGGGCATTGAGCACAATCGACGCGGAGATGGTTTTATCGTAATCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCCAGATCGTGGAACAACGCGCCGCCACCGCTACTGCGCCGTGCCAAGCGGGAGTTATCTTCTTACATCCCCCGGGGATTACCCTCCTTCCACGGGTTCTGCGCCCGAcc                                                                                                                                                                              >  2:130753/5‑251 (MQ=255)                                                                    ttCACCACCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATcacta                                                                                                                                                                      <  1:199488/250‑1 (MQ=255)                                                                      caccacCAGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTCCACTCTTTCCACGTGTTCTGCGCGCGAcccctcactaccg                                                                                                                                                                   >  1:226647/1‑240 (MQ=255)                                                                         caccaGATCGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCgtgt                                                                                                                                                                >  1:37029/1‑251 (MQ=255)                                                                                 cgtttcgccctgACGCTTCGGCCCTGACGCCAAGCGCGTTCAGCGCATTGAGCACAAGCGGCGTTTAGATGGTTTTTTCGTACTCCGGCTTGCCTGCCATAAAGGTAAAGCATGTATTGCCGAGATCGTGGAACACCTCGCCGCCACCGCTACTGCTCCGTTCCTGGCGGACGTGATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCTCTAGCGTGTCTGCATTg                                                                                                                                                        <  2:256006/240‑1 (MQ=255)                                                                                 cGTTACGCCCGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCAGTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGGCGGCAtt                                                                                                                                                         <  2:237684/250‑1 (MQ=255)                                                                                         acGGACGCTTCGGCGCTGACGCCAAGCGCGTTCAGCGCCTTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCCTTGCGCCag                                                                                                                                                  >  2:1576/2‑249 (MQ=255)                                                                                              cGCTTCGGCGCTGAGGCCAAGCGCTTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTATCCCGGCTGGCCAGCCGTAAAGGTAAAGCAGGTATTGCCGATATCGTGGAACACCGCTCCTCCACCGCTACTTCGCCGTGCCAGGCGGACGTGATCTTCTACCATCCTCCGGGTATTACACTCTTTCGATGTGTTCTGCGCGCTACCAAACACTACGGTGTCGGCATTTCGTCTGAGAAACa                                                                                                                                           <  2:123214/251‑1 (MQ=255)                                                                                                              ccaaGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATAGCCGAGATCGTGGAACACCGCGCCGCCACCGCTTCTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCAGCCGCCGGGTATTACACTCTTGCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGc                                                                                                                           <  1:162035/248‑1 (MQ=255)                                                                                                                          aGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCg                                                                                                               <  1:195386/251‑1 (MQ=255)                                                                                                                               aTTGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACGCTTTCCACGGGTTCTGCGCTCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGaaaaa                                                                                                          <  1:126916/251‑1 (MQ=255)                                                                                                                                ttGAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCGGGCGAAAAAt                                                                                                         >  2:189357/1‑251 (MQ=255)                                                                                                                                  gAGCACAATCGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATac                                                                                                       <  1:193030/251‑1 (MQ=255)                                                                                                                                           cGACGTGGAGATGGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAGTACACtcctccc                                                                                              >  2:256462/1‑250 (MQ=255)                                                                                                                                            gACGTGGAGATGGTTTTAGCGTACGCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTGTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGTCGAAAAATACTCTCCTCCAc                                                                                             <  2:63248/251‑1 (MQ=255)                                                                                                                                                        gTTTTATCGTACGCCTGCTGGCCAGCCAGAAAGGTAAAGCAGGTATGGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTaaa                                                                                 <  2:25027/251‑1 (MQ=255)                                                                                                                                                                  acgccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCg                                                                       <  1:256462/248‑1 (MQ=255)                                                                                                                                                                    cccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTa                                                                     <  1:223241/250‑1 (MQ=255)                                                                                                                                                                     cccgCTTGCCAGCCATAAAGGTAAAGCAGGGATTGCCGAGATCGGGGAACACCGCGCCGCCCCCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCAGCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACGCCTCCTCCTCCAGGTTAAACCACGGGTCGt                                                                      <  2:37029/246‑1 (MQ=255)                                                                                                                                                                     ccGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGt                                                                      <  1:33938/249‑1 (MQ=255)                                                                                                                                                                       ggCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACCCCGCGCCGCCACCGCTACTGCGCCGTGACAGGCGGACGTTATCTTCTTCCATCCGCCGGGCATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCTCTACCGTGTCGGCATTGCGCCAGCGCAACCGATCGCGCTGCGTGCCGGGCATCTGGCGAAAAATACACTCTTCCACCGCAAGGTTATACCACgtggctta                                                                     >  1:230059/1‑241 (MQ=255)                                                                                                                                                                        gCTTGCCAGCCAGAAAGGTAAAACAGGTAGTGCCGAGATCGTGGAACCCCGCGCCGCCACCGCTACGGCGCCGTGCCAGGCGGACGTTATCTTCTTCCCTCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAagag                                                                 <  2:151641/251‑1 (MQ=255)                                                                                                                                                                            gccagccaGAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTTCTGCTCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCa                                                              <  2:65861/250‑1 (MQ=255)                                                                                                                                                                                       aGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAgg                                                 >  2:205816/1‑251 (MQ=255)                                                                                                                                                                                         gTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCCCCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCATAGATG‑AGCAGGCg                                               <  1:201172/251‑1 (MQ=255)                                                                                                                                                                                                      tGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACAt                                   >  1:77225/1‑250 (MQ=255)                                                                                                                                                                                                       gCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATaa                                 <  1:130753/251‑1 (MQ=255)                                                                                                                                                                                                                       cccGCGCCGCCACCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTgt                 <  1:240085/250‑1 (MQ=255)                                                                                                                                                                                                                       ccagcgcCGCCACCGCTCCGGCGCCGTTCCAGGCGGACGTTAACTTCTGCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTGCCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCAGCTGGCGAAAAATACACTCCGCCACCGCTAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCATGCGTAATGTGGACATAACGATTTCCTTGCTTgt                 <  2:264601/248‑1 (MQ=255)                                                                                                                                                                                                                                  aCCGCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACCGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGACAGAGATGCACCAG‑CGTAATGTGGCCATAACGATTTCCTTTCTTGCAGTAGTGGGt       >  2:8567/1‑250 (MQ=255)                                                                                                                                                                                                                                     gCTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGGGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGTTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCAct   <  1:237355/251‑1 (MQ=255)                                                                                                                                                                                                                                      cTACTGCGCCGTGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCActc  >  1:166925/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 GTGGAAGCCGCGATCTTTGGTTTCGCGATAGGCCGAGCCTGAGACTTTGCGGTCGCCTTCGACGGTTTTCACCACCAGATCGTTACGTCCGGACGCTTCGGCGCTGACGCCGAGCGCGTTCAGCGCATTGAGCACAATCGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTC  >  NC_000913/4623671‑4624150

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