Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,633	A→C	intergenic (+86/+222)	nadR → / ← ettA	trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor/energy‑dependent translational throttle A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,628,633	0	A	C	97.5%	119.4 / NA	41	intergenic (+86/+222)	nadR/ettA	trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor/energy‑dependent translational throttle A
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (26/13);  total (27/14)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACT  >  NC_000913/4628406‑4628859                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   tCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGt                                                                                                                                                                                                                             <  1:24979/235‑1 (MQ=255)  cGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGttt                                                                                                                                                                                                                           >  1:7424/1‑236 (MQ=255)               aGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTg                                                                                                                                                                                                                <  2:300993/234‑1 (MQ=255)                      gCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTATACAAGGCTTGAAACAACCTGCACgggct                                                                                                                                                                                                        >  2:554120/1‑233 (MQ=255)                                             accatatCGAATTTGTGCGGGTAGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTAAAACCCTCTCTGGAAGAGGCTGGCGAGATTACACAAGGTTTGAAAAAACCTGCCCGGGTGCCCCGCCACGaaa                                                                                                                                                                                           <  1:598635/222‑1 (MQ=255)                                                    cGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGAACCGCCACGAAATTGAGGa                                                                                                                                                                                    >  2:598635/1‑225 (MQ=255)                                                     gatttggtgCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGATTTGAAACCCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCCCGGGTGCCCCGCCACGAAATTGAGGATGCTAAgg                                                                                                                                                                            <  1:645305/226‑1 (MQ=255)                                                                ggTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAAGTTTGAAAAAACCTGCACGAGTGACCCGCCACGAAATTGAGGATGCTCAGGCTGTTCTAGACCggc                                                                                                                                                             >  2:364891/1‑235 (MQ=255)                                                                 gTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTg                                                                                                                                                                                                                >  2:766510/1‑184 (MQ=255)                                                                 gTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTg                                                                                                                                                                                                                <  1:766510/184‑1 (MQ=255)                                                                 gTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCAGGAAATTGAGGAGGCTAAGGCGGTTGTGAACCggga                                                                                                                                                            >  2:38960/1‑234 (MQ=255)                                                                           aaGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGaa                                                                                                                                                                                                                        <  1:406683/166‑1 (MQ=255)                                                                           aaGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGaa                                                                                                                                                                                                                        >  2:406683/1‑166 (MQ=255)                                                                                  cGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTAAAGCCGCa                                                                                                                                           >  2:68113/1‑236 (MQ=255)                                                                                      tgtcgtcTCGAGCGATGCGTGTAATGGGTGCGGAAGAGGATGGGGGAGCAGAG‑TAACCAGTGGTAAAACTGCTCAAAGGCGGGG‑TAAAAAGGGAGTTTTTGTGTAGAGTTAAAACCCTCTCGGGAAGAGGCTGCGGAGATTACAAAGGGTTTAAAAAACCCGGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATcag                                                                                                                                       <  2:747009/235‑3 (MQ=255)                                                                                        tCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCgg                                                                                                                                      >  1:10599/1‑235 (MQ=255)                                                                                                  cgctgcgTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGAACCGCCACGAAATTTAGGATGCTAAGGCTGTTGTTGACCGGCAAAGGTGTTCATGCCGCATCCAGCATTCGGTGCt                                                                                                                           >  2:587006/1‑236 (MQ=255)                                                                                                       cgTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGGGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGGGTTCACGCCGCATCCGGCATTCGGGGTTCa                                                                                                                         >  2:305849/1‑233 (MQ=255)                                                                                                                 gTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGc                                                                                                              >  1:373085/1‑234 (MQ=255)                                                                                                                      ggAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCtt                                                                                                       >  1:505128/1‑236 (MQ=255)                                                                                                                      ggAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCATGATGCGACGCtt                                                                                                       >  1:757203/1‑236 (MQ=255)                                                                                                                                    aGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGGGACCCGCCACGAAATTGAGGATGCTAAGGCt                                                                                                                                                                          >  1:370010/1‑155 (MQ=255)                                                                                                                                                c‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAATACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAAGTCCTGATGCGACGCTTACCGCGTCCTATCAGGCCTACAGttt                                                                             >  2:695270/1‑236 (MQ=255)                                                                                                                                                c‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGttt                                                                             >  2:631479/1‑236 (MQ=255)                                                                                                                                                     ataaaaCTGCTCAGAGG‑GAGGATATAAAATGAGTTTTTGTGACGAGTAGAAACCCTCTATGAAAGGCTCTGTCAAGTTTACAGAGGGTTAGAAAACCCCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACa                                                                         <  2:366460/233‑1 (MQ=255)                                                                                                                                                      tgaactggcggaaaGGCGAGG‑TATAAAATGTGTTTTTTGGAGGTGTGAAAACCCTCTGGGAAAGAGGCTGTCAAGATACAACAAGGTTTGAAAAACCCGGCACGGGTGCCCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTa                                                                                   <  2:314782/214‑1 (MQ=255)                                                                                                                                                              cTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGa                                                                >  1:99170/1‑236 (MQ=255)                                                                                                                                                                         agataaaagaTGATTTTTTTTGGTGATTAGAAACCCTCTGTGAAAGAGGCTGTCGAGATTACACAAGGTTGAAAAAAACCTGCACGGGTGACCCGCCTCGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGc                                                          <  1:554120/223‑1 (MQ=255)                                                                                                                                                                               aaaaTGATTTTTTTTGATGAGTTAAAACCCTCTATGGAGGGGGCTTGCGAGATTACAAAAGTTTAAAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTAt                                                                                           <  2:112402/193‑1 (MQ=255)                                                                                                                                                                                      gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCt                                                                                    >  1:112402/1‑193 (MQ=255)                                                                                                                                                                                      gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGa                                         >  1:74133/1‑236 (MQ=255)                                                                                                                                                                                      gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGGTTTCCAACGTATTGAAAATGCTGGGTTTGTAGGCCggc                                         >  2:33130/1‑235 (MQ=255)                                                                                                                                                                                      gTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTATTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGGTTTGTAGGctcgc                                         >  2:255693/1‑232 (MQ=255)                                                                                                                                                                                                              ctctGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGc                 >  1:704439/1‑236 (MQ=255)                                                                                                                                                                                                               tctGGAAGAGGCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCa                >  1:256927/1‑236 (MQ=255)                                                                                                                                                                                                                 gggAAGATGTTGTCAAGATTACACGAGTTGTGAAAAAGCCTGCCGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAtt              <  1:587006/235‑1 (MQ=39)                                                                                                                                                                                                                       agGCTGTCAAGATAACACAAGGTTAAAAAAACCCTGCAGGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACg         <  1:695270/235‑1 (MQ=255)                                                                                                                                                                                                                         gCTGTCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGc        >  2:500731/1‑234 (MQ=255)                                                                                                                                                                                                                           tGTCGAGATTACAAAGGGTTAAAAAACCCCTGCACGGGGAACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCgg                                                                                                                               <  2:57787/113‑1 (MQ=255)                                                                                                                                                                                                                             tCGAGATTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACt  >  1:227980/1‑236 (MQ=255)                                                                                                                                                                                                                                  aTTACACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCaa                                                                                                                        >  1:57787/1‑113 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                   TCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACC‑GTGATGAAACTGCTCAAAGGCGAGG‑TATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACT  >  NC_000913/4628406‑4628859

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