Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	504,319	A→G	G126G (GGT→GGC)	fsr ←	fosmidomycin efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	504,319	0	A	G	97.8%	162.7 / ‑3.6	46	G126G (GGT→GGC)	fsr	fosmidomycin efflux pump
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (17/28);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

GTGCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGA  >  NC_000913/504087‑504553                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           gtgCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTa                                                                                                                                                                                                                                      <  2:194847/239‑1 (MQ=255)  tgCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTAcc                                                                                                                                                                                                                                    >  2:40818/1‑240 (MQ=255)           ccAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAgcg                                                                                                                                                                                                                            <  2:567016/239‑1 (MQ=255)              aCGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGCGCCGGTACCGACCAGcgccgc                                                                                                                                                                                                                        >  1:542463/1‑240 (MQ=255)               cGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGcgccgcc                                                                                                                                                                                                                       <  1:370617/240‑1 (MQ=255)               cGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGcgccgcc                                                                                                                                                                                                                       <  1:52557/240‑1 (MQ=255)                            ccAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACt                                                                                                                                                                                                           >  1:524800/1‑239 (MQ=255)                                              gcagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAACCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGcc                                                                                                                                                                                         <  1:550160/239‑1 (MQ=255)                                                     cGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAag                                                                                                                                                                                 <  2:143174/240‑1 (MQ=255)                                                       cAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAag                                                                                                                                                                                 <  2:483383/238‑1 (MQ=255)                                                                          gTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGt                                                                                                                                                            <  1:493236/240‑1 (MQ=255)                                                                            tGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTg                                                                                                                                                                                                          >  2:147196/1‑192 (MQ=255)                                                                                   gCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGcaca                                                                                                                                                    >  2:205034/1‑239 (MQ=255)                                                                                    ccATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATg                                                                                                                                                  <  1:243202/240‑1 (MQ=255)                                                                                        aaGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCaa                                                                                                                                              <  2:206966/240‑1 (MQ=255)                                                                                                                gCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAgcg                                                                                                                                                                                     <  2:60160/177‑1 (MQ=255)                                                                                                                      gTCCCAGGGAACTGCCAAAGTTGCCGCCGACNTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGa                                                                                                                                                                       <  2:19943/185‑1 (MQ=255)                                                                                                                      gTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTAt                                                                                                                >  1:376218/1‑240 (MQ=255)                                                                                                                      gTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTAt                                                                                                                <  1:310721/240‑1 (MQ=255)                                                                                                                       tCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTAt                                                                                                                <  2:462065/239‑1 (MQ=255)                                                                                                                                aCTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGGAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGt                                                                                                       >  2:417481/1‑239 (MQ=255)                                                                                                                                                ccgACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTg                                                                                       <  1:564230/239‑1 (MQ=255)                                                                                                                                                            atagatTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCAt                                                                                                                           <  1:457970/191‑1 (MQ=255)                                                                                                                                                            atagatTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAagag                                                                          <  1:228501/240‑1 (MQ=255)                                                                                                                                                             tagatTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAgg                                                                         >  1:335661/1‑240 (MQ=255)                                                                                                                                                              agatTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAagag                                                                          >  2:579255/1‑238 (MQ=255)                                                                                                                                                                     gcCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTg                                                                  <  1:197404/239‑1 (MQ=255)                                                                                                                                                                      ccAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTgg                                                                 <  1:510753/239‑1 (MQ=255)                                                                                                                                                                      ccAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTgg                                                                 <  1:493575/239‑1 (MQ=255)                                                                                                                                                                      ccAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTg                                                                  <  1:147406/238‑1 (MQ=255)                                                                                                                                                                           ccATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGt                                                            <  1:566830/239‑1 (MQ=255)                                                                                                                                                                                   cccgccGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTa                                                    <  2:203638/239‑1 (MQ=255)                                                                                                                                                                                   cccgccGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTAt                                                   >  2:25103/1‑240 (MQ=255)                                                                                                                                                                                    ccgccgGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAg                                                                                                     >  1:489785/1‑189 (MQ=255)                                                                                                                                                                                    ccgccgGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAg                                                                                                     <  2:489777/189‑1 (MQ=255)                                                                                                                                                                                     cgccgGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCa                                                 >  1:396329/1‑240 (MQ=255)                                                                                                                                                                                     cgccgGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCa                                                 >  2:349685/1‑240 (MQ=255)                                                                                                                                                                                         ggAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACCGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGcc                                             <  2:490018/240‑1 (MQ=255)                                                                                                                                                                                             gCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAAt                                          <  2:563579/239‑1 (MQ=255)                                                                                                                                                                                                           gcgAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCa                            >  2:578738/1‑239 (MQ=255)                                                                                                                                                                                                                aagaTTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGaaa                      >  1:83640/1‑240 (MQ=255)                                                                                                                                                                                                                aagaTTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGaaa                      >  1:612286/1‑240 (MQ=255)                                                                                                                                                                                                                   aTTCCGGATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAAtt                   <  1:40819/240‑1 (MQ=255)                                                                                                                                                                                                                         gATGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTctgac              >  1:31406/1‑239 (MQ=255)                                                                                                                                                                                                                           tGAAAGACCGATGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAAt                                          <  2:359907/209‑1 (MQ=255)                                                                                                                                                                                                                                     aTGAGCCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGa  >  1:352604/1‑239 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                           GTGCCGAGTACCAACGGCTGATTTGCGCCAACACCACGATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGAAAGACCGATGAACCGGTACCGACCAGCGCCGCCGCCAGCAGAACTGCGCCAAAACTGCCCGCCAGCGCAAGCAGCACCAGACCGCTTAAGGTAAAGCACATGCCAATTGGCAACGACCATGGCATCGGATATTTATCGGTCCAGTAGCCGACCACTGGTTGCAGTAGCGAAGAGGCGAGCTGGAAGGTGAGGGTTATCATGCCAATCTGCATAAATGTCAGAGAAAATTCTGACTGAAGCAGCGGA  >  NC_000913/504087‑504553

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