Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,624,055	T→C	E21E (GAA→GAG)	lplA ←	lipoate‑‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,624,055	0	T	C	97.4%	117.9 / ‑6.1	38	E21E (GAA→GAG)	lplA	lipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/23);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGT  >  NC_000913/4623809‑4624295                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  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)                                                                                                     cgtgccAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTtctc                                                                                                                                          >  1:247931/4‑251 (MQ=255)                                                                                                           aGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCttttaaaagagtgac                                                                                                                                    <  1:89182/251‑12 (MQ=255)                                                                                                                          tcttcCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCa                                                                                                                      <  1:205816/250‑1 (MQ=255)                                                                                                                             tcCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCAcc                                                                                                                    <  1:149131/249‑1 (MQ=255)                                                                                                                                   cgccgGGTATTACACTCTTTCAACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCCCCGGTTTGCt                                                                                                            >  2:59681/1‑251 (MQ=255)                                                                                                                                              acacacTTTCCACGGGTTCTGCGCGCGACCAATCACTACCGTGTCGGCATTGCGCCAGAGCAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAATTACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGCAATGGGTCACTCTTTTTTTTCGATTTCTTCCTCTTCCGGCACCGGGTGTCTGGCGGtttac                                                                                                  >  2:150506/1‑247 (MQ=255)                                                                                                                                                   ctTTCCACGGGTTCTGCGCGCGACCAATCACTGCCTTGTCGGCAGTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAATCCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGTGTCACTCGTTTTACTCGCTGTCTTCCGCTTCCGGCACCGGTTTGCTGGCGTTTACCAGGa                                                                                              <  2:136126/249‑1 (MQ=255)                                                                                                                                                                        gACCAATCACTTCCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTTGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGta                                                                       <  2:195214/251‑2 (MQ=255)                                                                                                                                                                                 ctacCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTTACAGGGGGGGCGATTGCTGCCAGCGGGGCCGTTTAccc                                                              >  2:219822/1‑251 (MQ=255)                                                                                                                                                                                   acCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAGAGGCGCTTGCTGCCAGCGGGGACGTTTATCCTg                                                            >  2:219282/1‑251 (MQ=255)                                                                                                                                                                                          cGGCATTGCGACAGAGAAACAGAACGCGCTGCGGGGCGGGCATCTGGCGAAAAATACACTACCCCACCGCCAGGGTAAAACACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTGACCCTGTagcag                                                      <  2:226647/250‑1 (MQ=255)                                                                                                                                                                                                 gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTTTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCgg                                               <  1:120080/250‑1 (MQ=255)                                                                                                                                                                                                 gcgcCAGAGAAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCgg                                               <  1:267869/250‑1 (MQ=255)                                                                                                                                                                                                   gcCAGCGAAACAGAACGGCCTGCCTGGCGGGCATCTGGCGTAAAATACACTCCTCCCCCGCCCGGTTAAACCACGGGTCGTGAGAGTCAGAGATG‑CGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGCAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACCGGGAAGGCGCTTACTGCCCGCGGGTAAGTTTAccccatcgccgcaggcgcgtt                                            >  2:20263/1‑233 (MQ=255)                                                                                                                                                                                                        agaAACAGAACGCGCTGCGTGGCGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGGGGGCGGGTccgccc                                       >  2:266490/1‑246 (MQ=255)                                                                                                                                                                                                                              cGGGCATCTGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGGACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGGTTGcccgggaa                    >  2:178166/1‑242 (MQ=255)                                                                                                                                                                                                                                   atctGGCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTTCCAGTGAAAgcagcac            >  1:245923/4‑250 (MQ=255)                                                                                                                                                                                                                                       ggCGAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCGTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATc        <  1:178166/251‑1 (MQ=255)                                                                                                                                                                                                                                          gAAAAATACACTCCTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATccgg     >  1:38959/1‑248 (MQ=255)                                                                                                                                                                                                                                             aaagacacTCCTCCTCCGCCAGGTTAAACCACGGTTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTCTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCATCCGGCAGCGGTGTGCTGGCGGTTAAAAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTTCCCTGTAGCAGCGTGCGGGGCAGCACTAAGCCGATTGCCAGTGATAGCAGCACCATTAGGCGt  <  2:230059/247‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  CGACGTGGAGATAGTTTTATCGTACTCCGGCTTGCCAGCCATAAAGGTAAAGCAGGTATTGCCGAGATCGTGGAACACCGCGCCGCCACCGCTACTGCGCCGCGCCAGGCGGACGTTATCTTCTTCCATCCGCCGGGTATTACACTCTTTCCACGGGTTCTGCGCGCGACCAATTACTACCGTGTCGGCATTGCGCCAGAGAAACAGAACGCGCTGCGTGGCGGGCATTTGGCGAAAAATACACTCTTCCACCGCCAGGTTAAACCACGGGTCGTAAGAGTCAGAGATG‑AGCAGGCGTAATGTGGACATAACGATTTCCTTTCTTGTAGTAATGGGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGT  >  NC_000913/4623809‑4624295

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