Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,888,554	T→C	7.1%	E398G (GAA→GGA)	fadD ←	acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,888,554	0	T	C	7.1%	191.4 / 5.0	56	E398G (GAA→GGA)	fadD	acyl‑CoA synthetase (long‑chain‑fatty‑acid‑‑CoA ligase)
Reads supporting (aligned to +/- strand):  ref base T (21/31);  new base C (2/2);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

TGCATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAA  >  NC_000913/1888320‑1888792                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               tgcATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCAtt                                                                                                                                                                                                                                            <  1:1085845/239‑1 (MQ=255)   cATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTa                                                                                                                                                                                                                                           >  1:248839/1‑238 (MQ=255)     tGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTCCATTatca                                                                                                                                                                                                                                        >  1:1289429/1‑239 (MQ=255)           cATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTCCATTatcatcatca                                                                                                                                                                                                                                  >  1:703910/1‑239 (MQ=255)               ttCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATcca                                                                                                                                                                                                                              <  1:279794/239‑1 (MQ=255)                tCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATcca                                                                                                                                                                                                                              >  2:279794/1‑238 (MQ=255)                         ttGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTg                                                                                                                                                                                                                      >  1:733446/1‑237 (MQ=255)                              aTAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGCCTTcc                                                                                                                                                                                                                >  2:1325059/1‑238 (MQ=255)                                      tAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGTAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACgg                                                                                                                                                                                                        <  1:968830/238‑1 (MQ=255)                                       aaaaCCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGc                                                                                                                                                                                                       >  1:311433/1‑238 (MQ=255)                                             ggAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCa                                                                                                                                                                                                <  1:67986/239‑1 (MQ=255)                                              gAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCATCTGCGGTACTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCa                                                                                                                                                                                                >  1:1236342/1‑238 (MQ=255)                                              gAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCaa                                                                                                                                                                                               >  2:1218055/1‑239 (MQ=255)                                                aaCCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTCCCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAAc                                                                                                                                                                                             <  2:1190140/239‑1 (MQ=255)                                                aaCCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAAc                                                                                                                                                                                             <  2:248839/239‑1 (MQ=255)                                                      aaTCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGAt                                                                                                                                                                                         >  1:1203680/1‑237 (MQ=255)                                                       aTCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGc                                                                                                                                                                                       <  1:1325059/238‑1 (MQ=255)                                                       aTCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGc                                                                                                                                                                                       <  1:649579/238‑1 (MQ=255)                                                       aTCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGc                                                                                                                                                                                       <  1:560936/238‑1 (MQ=255)                                                       aTCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCt                                                                                                                                                                                      <  1:1366667/239‑1 (MQ=255)                                                           tGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTAcc                                                                                                                                                                                   >  2:711063/1‑238 (MQ=255)                                                                tttttACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATg                                                                                                                                                                             <  2:1006813/239‑1 (MQ=255)                                                                tttttACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATg                                                                                                                                                                             <  2:191909/239‑1 (MQ=255)                                                                   ttACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGTGGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCAACGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATa                                                                                                                                                                          <  2:1027563/239‑1 (MQ=255)                                                                     tcGATCGACAATGCGCAGGAATCCGTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAAt                                                                                                                                                                        <  2:1234546/238‑1 (MQ=255)                                                                     aCGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAAt                                                                                                                                                                        <  2:544552/239‑1 (MQ=255)                                                                     aCGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAAt                                                                                                                                                                        <  2:733446/239‑1 (MQ=255)                                                                         tCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCaa                                                                                                                                                                     <  2:945242/238‑1 (MQ=255)                                                                         tCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCaa                                                                                                                                                                     <  2:25726/238‑1 (MQ=255)                                                                         tCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTCCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAat                                                                                                                                                                    <  2:703910/239‑1 (MQ=255)                                                                               aTGCGCAGGCATCCTTCTTCATCCATTCCCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCata                                                                                                                                                              <  1:456220/239‑1 (MQ=255)                                                                               aTGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCata                                                                                                                                                              <  2:311433/239‑1 (MQ=255)                                                                               aTGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCata                                                                                                                                                              <  1:945152/239‑1 (MQ=255)                                                                                                      cATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTCCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAgcggcg                                                                                                                                       <  2:1289429/239‑1 (MQ=255)                                                                                                          aCCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGAGcaca                                                                                                                                   <  2:1203680/239‑1 (MQ=255)                                                                                                                 tGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGt                                                                                                                             <  1:1218055/238‑1 (MQ=255)                                                                                                                                    gCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCcagca                                                                                                         >  1:412475/1‑239 (MQ=255)                                                                                                                                    gCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCcagca                                                                                                         >  2:275231/1‑239 (MQ=255)                                                                                                                                        tttttGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGAt                                                                                                      >  2:181939/1‑238 (MQ=255)                                                                                                                                          tttGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAat                                                                                                                                                                    <  1:1087214/174‑1 (MQ=255)                                                                                                                                          tttGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAat                                                                                                                                                                    >  2:1087214/1‑174 (MQ=255)                                                                                                                                                             gCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCa                                                                                >  2:203137/1‑239 (MQ=255)                                                                                                                                                                tCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAAc                                                                              >  2:1247056/1‑238 (MQ=255)                                                                                                                                                                                aaCCCAGCATCTCCTGCGGTCCTTTGACACAAAACCCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAACGCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTtg                                                              <  2:933787/238‑1 (MQ=255)                                                                                                                                                                                aaCCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTtg                                                              <  2:1196368/238‑1 (MQ=255)                                                                                                                                                                                   ccAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTtgctg                                                           >  2:1178058/1‑238 (MQ=255)                                                                                                                                                                                       cATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTg                                                      >  1:1030711/1‑239 (MQ=255)                                                                                                                                                                                         tCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGc                                                    <  1:711063/239‑1 (MQ=255)                                                                                                                                                                                                    ccTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGCCGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGc                                          <  1:208874/238‑1 (MQ=255)                                                                                                                                                                                                               aaGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGc                                                                            >  2:1293904/1‑193 (MQ=255)                                                                                                                                                                                                               aaGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGc                                                                            <  1:1293904/193‑1 (MQ=255)                                                                                                                                                                                                                               ccTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAAt               >  1:492382/1‑238 (MQ=255)                                                                                                                                                                                                                                   gtggtACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCa            >  2:1067895/1‑237 (MQ=255)                                                                                                                                                                                                                                        aCTTCATTATCATCATCATCCCCCAGTTTGGCTTCCGTCGACGGCACCGGCAACCCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGctgctg     <  1:95592/239‑1 (MQ=255)                                                                                                                                                                                                                                         cTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTAc                                                                                                                                                                                    >  1:1371475/1‑63 (MQ=255)                                                                                                                                                                                                                                         cTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTAc                                                                                                                                                                                    <  2:1371475/63‑1 (MQ=255)                                                                                                                                                                                                                                           tCATTATCATCATCATCCACCAGTTTGTCTTCCGTCGACGGCACCGGCAAGCCGATGCTACCTCTATGATAATCAATATCATATGGGTTAACGCTGACCGGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGaa  <  2:1102719/239‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                               TGCATGACGACATCTTCAATCTCGTTGGGATAGACGTTAAAACCGGAAACCAGAATCATGTCTTTTTTACGATCGACAATGCGCAGGAATCCTTCTTCATCCATTACCGCGATGTCGCCGGTGTGTAACCAGCCATTTTTGATGATTTCATCGGTAGCATCGGGACGCTGCCAGTAACCCAGCATCACCTGCGGTCCTTTGACACAAAGCTCACCCGGTTGACCTGGTGGTACTTCATTATCATCATCATCCACCAGTTTGGCTTCCGTCGACGGCACCGGCAAACCGATGCTACCACTATGATAATCAATATCATATGGGTTAACGCTGACCAGCGGCGCACACTCGGTAAGGCCATAGCCTTCCAGCAGATACTGTCCGGTCAGTTTCACCCAACGCTCTGCCACCACTTGCTGCACTGGCATCCCACCGCCTGCGGAAAGATGCAGACTGGAGAAATCCAGCTGCTGGAA  >  NC_000913/1888320‑1888792

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