Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,318,215	G→A	L184L (CTG→TTG)	phnK ←	carbon‑phosphorus lyase subunit PhnK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,318,215	0	G	A	86.8%	82.0 / 2.4	40	L184L (CTG→TTG)	phnK	carbon‑phosphorus lyase subunit PhnK
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (15/18);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

GCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTG  >  NC_000913/4317970‑4318460                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       gCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCAGCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCggg                                                                                                                                                                                                                                                   <  1:151021/251‑1 (MQ=255)   cGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATAAGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGcc                                                                                                                                                                                                                                                 >  1:213150/1‑251 (MQ=255)     ccctaacacgacgccctgccgatctcACCGCCTCAATTCTGCAAAACCGATGCCCCCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTGAACCCACTCTCCACCACGTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTg                                                                                                                                                                                                                                               <  2:113475/226‑1 (MQ=255)        aaGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCAc                                                                                                                                                                                                                                            >  1:75617/1‑251 (MQ=255)          ggCGTTTGCGCCTCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTTTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCCCCCCTTG‑CCCCTGTTGCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGTGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCg                                                                                                                                                                                                                                          <  1:236431/251‑1 (MQ=255)             gTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGaca                                                                                                                                                                                                                                       <  1:327502/251‑1 (MQ=255)                     cgcACCCGGCTCCGGCCAAATTCTGCAAACCCGTTGACCCCATCAGCTGTGTAAACTGAGGAGGCGGGTCTTCGAGCACGCGGTCGTTTAACCCACTCGCCACCACTTG‑CCCCCTCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACTCCTAAATCATGGGTGACAATCCCCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGc                                                                                                                                                                                                                               <  2:368309/251‑1 (MQ=255)                           cGGCACCGGCCCAATTCGTCAAACCCGAAGACACCATCAGCTGTGTATACGGATGATGCGGGCCGCCGAGCCCCCGGTCGGTTAACCCACTCTCCCCCACTTG‑CCCCTGCTTCACCACCAGCAAACGTTCCGCCGGCAGGCGGGCGACGCCTAAACCATGGGTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGGGCAAGCGGGCCTGCACCGACACATCCAGCccgccg                                                                                                                                                                                                                         <  2:379818/251‑1 (MQ=255)                              cACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTg                                                                                                                                                                                                                                               >  1:113475/1‑226 (MQ=255)                                                             cagcagCTGTGTATCCGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCg                                                                                                                                                                                       <  1:130190/251‑1 (MQ=255)                                                                cagcTGAGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGGTCAGCTCCCCCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCggctgagtc                                                                                                                                                                                     >  2:348065/1‑243 (MQ=255)                                                                  ggtgtgTATCCGGATGAGGCGGCTCGTCGAGCATGCGGTCGGTTACCCCACTCGCCACC‑CTTGGCCCCGCCTTCATCACCAGCAAACGGTTCGCCAGAAGGCGCGCGCCGCCTAAATAATGGGTGACAATCACCACCGCGAGGTTCCGCTCCACCACCAGGCCGCGCAGCAGGTCGATCAAGCGGGCTTGCACAGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCAcc                                                                                                                                                                                  <  1:126416/249‑1 (MQ=255)                                                                   cTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCAcc                                                                                                                                                                                  >  1:269915/1‑250 (MQ=255)                                                                       gtATTCGGGTGATGCGGGCCGTCGTGCACCCGGTCGGTTAACCCACTCCCCACCCCTTG‑CCCCTGTTTCATCCCCAGCAAACGTTCCCCCAGCGGGCGGGCGCTGCCTAAACCTCGGTTGACAATCCCCACCGCTAGTTACATCTCCTCCCCCAGGCCGCCCACCAGTCCGACCAAGCGTTCATTCACCGTCACGCCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCATCGGATTCGTCACCAGGtt                                                                                                                                                                             <  2:237140/251‑1 (MQ=255)                                                                            cGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGGGACAATCACCACCGCGAGGTTCAGCACCCCCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACCCATCCAGCCCGCCGGTCGGTTCATCCATAAACAACAACTTAGGATGCGTCACCCGGGTGCggg                                                                                                                                                                        >  2:294795/1‑251 (MQ=255)                                                                                       ggTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCCCCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCaaa                                                                                                                                                             <  1:241366/251‑1 (MQ=255)                                                                                          cgtcgAGCACGCGGTCGGTTAACCCACTCCCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACgc                                                                                                                                                          <  1:166978/251‑1 (MQ=255)                                                                                                         cgggtAACCCACTCTCCCCCACTTG‑CCCCTGCTTCATCACCAGCAAACGGCCCGCCACCAGCCGGGCGACCCCTAAACCAGGGTTGACAATCACCACCGCTGGGTTCAGCTCCACCACCAGGCCGCGCAGCAGTTCGATCAAGCGGGCCTGCCCCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCATCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATAccgcc                                                                                                                                            <  1:234982/246‑1 (MQ=255)                                                                                                            ttAACCCACTCTCCACCACTTG‑CCCCTGCTTCCTCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAACTCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCTGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGa                                                                                                                                         >  1:74160/1‑250 (MQ=255)                                                                                                                  cACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAgggg                                                                                                                                  >  1:261689/1‑249 (MQ=255)                                                                                                                           ccacTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCACGCGGGCCTGCACCGACACATCCAGCCCGCCGGGCGGTTCCTCCCTAAACACCAGCTTCGGAGGCCTCCACAGGTTGCGGGCAATCTGCATACCCTGCTGCATACCCCCGGAACAGGGGGTGGGCAgg                                                                                                                         >  2:193434/1‑251 (MQ=255)                                                                                                                            cacTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGgt                                                                                                                        <  1:101211/251‑1 (MQ=255)                                                                                                                               ttG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCCCCGCCACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCCTCGGATGCGTCCCCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGCAACGGGGGTCGGcggggcg                                                                                                                      >  2:121083/1‑244 (MQ=255)                                                                                                                                        gcttaatcccaagaaacggctctcccAGCAGGGGGGCTACGCCTAAACCAGGGGGGACAATCACCACCTCGGGGTACAGCCCCACCACCAGGCCGCGCAGCGGGTTGGGCAGGCGGTCCTGCCCCGACACACCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATcc                                                                                                               <  2:124238/225‑1 (MQ=255)                                                                                                                                             atcaCCAGCAAACGGTCCGCCAGAAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCACGCGGGCCTGCACCGACACATCAATCCCGCCCGTAGGTTCATCCATAAACACCAGCTTCGGCTGCGTCACCAAGATGAGGGCAGTCTGCAAACCCTGCTGCCTACAGCCGGCAAAGGGGGTCGGCAGGACGTCGAACAGGTTggg                                                                                                        >  2:136390/1‑250 (MQ=255)                                                                                                                                                      aagcagcccGCCATCAGGCGGGCGGTGCCGAAATCATGGGTGACAATCGCCACCGGGAGGTTCAGCTCCACCGCCAGGCCGCGCAGCAGGTCGGGGAGGGGGGCCTGCCCCGACAGATCCTTCCCGCCGGTCGGTTCATCCATAATCACCAGCTTCGGATGCGTCCCCAGGTTGCGGGCAATCTGCAAACGCTTCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAAtctc                                                                                               <  1:347633/244‑1 (MQ=255)                                                                                                                                                          ggTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCCCCACCGCGAGGTTCAGCTCCCCCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCAcc                                                                                           <  1:282280/251‑1 (MQ=255)                                                                                                                                                               gCCAGCAGGCGGGCGACGCCTAAATCATGGGGGACAATCACCACCGCGAGTTTCAGCTCCCCCCCCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTc                                                                                      <  1:294795/251‑1 (MQ=255)                                                                                                                                                                ccAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTcc                                                                                     >  1:165283/1‑251 (MQ=255)                                                                                                                                                                    gagtggggcgACGCCAAAATCATGGGTGACAATCCCCCCCGCGAGTTTCAGCTCCCCCCCCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACACCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGcc                                                                                 <  1:336870/246‑1 (MQ=255)                                                                                                                                                                                      aTCATGGGAGACAATCACCACCGCGCGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCATCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGCATCTCCACCTCTTCCAGCCCCTTCTGCGCGGGGGCAc                                                               >  1:298259/1‑251 (MQ=255)                                                                                                                                                                                        cATGGGAGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACCTCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCAGATGCCACACCAGGTTGCGGGCAATCTGCAAACCCTGCTGCATACCGCCGGCAAAGGTGGTCGGCAGGTCGTCGCTCCGGTTGGCGGGGATCTCCACCCCATCCAGCCCCTTCTGCGCCGTgtccac                                                              >  2:184726/1‑245 (MQ=255)                                                                                                                                                                                             gTGACAATCACCCCCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATc                                                       <  1:348065/251‑1 (MQ=255)                                                                                                                                                                                                         accGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATGACg                                           >  1:304735/1‑251 (MQ=255)                                                                                                                                                                                                              gAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAAAGCTGCTGCATACCGCCGGAAAAGGAGGTCGGCAGGGAGCCGATCCGGTTGGCGGGAAGCACCAACTCTTACACCCACGTCTGCGCGGAGGAACGACAT‑TCGCCGTAAAGAcatgaac                                      >  2:325786/1‑245 (MQ=255)                                                                                                                                                                                                                          caccaccaGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATGACGTGCCCCGGTCGCCATCa                          <  1:235781/251‑1 (MQ=255)                                                                                                                                                                                                                           accaccaGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGTAGGCACGA‑ATATCGCCGTAATGACGTGCCCCGGTCGccctc                           >  1:289704/1‑246 (MQ=255)                                                                                                                                                                                                                                 aGGCCGCGCAGCAGGTCGAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTc                   >  1:332628/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGGTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGGCTGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATTACGTGCCCCGGGCGCCATCAGCCGCGCGCCGATATTGCCGCCTg  >  1:171031/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gAGCAAGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTg  >  1:381011/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GCCGGATAAGGCGTTTACGCCGCATCCGGCACCGGCTCAATTCTGCAAAACCGATGACACCAGCAGCTGTGTATACGGATGATGCGGGTCGTCGAGCACGCGGTCGGTTAACCCACTCTCCACCACTTG‑CCCCTGCTTCATCACCAGCAAACGGTCCGCCAGCAGGCGGGCGACGCCTAAATCATGGGTGACAATCACCACCGCGAGGTTCAGCTCCACCACCAGGCCGCGCAGCAGGTCGAGCAGGCGGGCCTGCACCGACACATCCAGCCCGCCGGTCGGTTCATCCATAAACACCAGCTTCGGATGCGTCACCAGGTTGCGGGCAATCTGCAAACGCTGCTGCATACCGCCGGAAAAGGTGGTCGGCAGGTCGTCGATCCGGTTGGCGGGAATCTCCACCTCTTCCAGCCACTTCTGCGCGGTGGCACGA‑ATATCGCCGTAATGACGTGCCCCGGTCGCCATCAGCCGCTCGCCGATATTGCCGCCTG  >  NC_000913/4317970‑4318460

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