Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,078,268	A→T	F205L (TTT→TTA)	putA ←	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,078,268	0	A	T	94.8%	163.9 / ‑4.5	59	F205L (TTT→TTA)	putA	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base T (22/33);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 7.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGT  >  NC_000913/1078028‑1078505                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               atCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGAc                                                                                                                                                                                                                                       <  1:614036/250‑1 (MQ=255)              cGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGCCCAAt                                                                                                                                                                                                                         >  2:340336/1‑251 (MQ=255)               ggAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGTTAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCCGCCCCCCAGTGGCGGCATTAACTAACAGTGACGGGCTACGACCaa                                                                                                                                                                                                                          >  2:672897/1‑249 (MQ=255)               ggAAACCTTTCTCTTACAGCTTGCGGACATTGACTACCGCTTCCGCGCTGGTTTCGACAGAGGCGATCTGCTCACCCATCAGGCGCAACGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAACTAGCGGGAGAGGCTGGCTTCGTTATGGGGGGAAACCCGTTTGCACGTAAACAGCAGCCCCCCGGGGGCGGCATTAACACCCAGTGACGGGcatcgacacatg                                                                                                                                                                                                                        >  2:262433/1‑240 (MQ=255)               ggAAACATTTCTCGCCCAGCTTGCGGGCATTGGCTCACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTGTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCGGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCACCCAGGTGGCGGCATTAACTAACAGTGACGGGTTACGACCAAtg                                                                                                                                                                                                                        <  2:57864/251‑1 (MQ=255)                       ttCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTAAGACCAATGTGTGACTg                                                                                                                                                                                                                >  1:40630/1‑251 (MQ=255)                         ctctTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGcc                                                                                                                                                                                                              >  1:208203/1‑251 (MQ=255)                           cgtCCAGCTGGCGGGCTTTGGCTAACGCTTCCGCGATGGTTCCGCCAGTGACAAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACttcc                                                                                                                                                                                                              <  2:539809/247‑4 (MQ=255)                                    ggCGGGCATTGCCTAACTGTTCCGCGATGTTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCCCACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGCCCAATGGGTGACTGCCAGT‑TACCGTTg                                                                                                                                                                                                  <  2:484820/250‑1 (MQ=255)                                     gcgcgCATTGGCTTACGCTTTCGCGATGGTTTCGCCAGGTCCGAACTGCTCTCCCATCATGCGCATCGACATATCCTCACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCTAGCGTGAGAGTCGGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTTGCGACATTAACTAACAAGGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTg                                                                                                                                                                                                  <  2:524843/246‑1 (MQ=255)                                     gCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTTCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGc                                                                                                                                                                                                 >  1:224589/1‑251 (MQ=255)                                     gCGGGCATTGGCTAACGCGTCCGCTATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTg                                                                                                                                                                                                  <  1:426695/250‑1 (MQ=255)                                      cGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCt                                                                                                                                                                                                >  1:428504/1‑251 (MQ=255)                                                  aacgcattcgcgATGTTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCACCGCCATATCCACACCTTTGCGGAGCAGCGGTTCACCTCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTTACGGGCGACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTcg                                                                                                                                                                                      <  2:56934/241‑1 (MQ=255)                                                  aaCGCTTCCGCTATGGTTGCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACCCCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAAGGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGTGCTTCGACCAATGTGTGACTGCCGGT‑TACCGTTGCTGATTTTGTcgcg                                                                                                                                                                                    <  2:245368/251‑1 (MQ=255)                                                    cGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGa                                                                                                                                                                                   >  1:471479/1‑250 (MQ=255)                                                     cctcccGCGATGTTTTCGCAAGTGACGAACTGCTCCCCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAAt                                                                                                                                                                                 <  2:386639/247‑1 (MQ=255)                                                      cTTCCGAGATGGTTGCGCCCGTGACGAACTGCTCACCCATCAGGCGCCTCGCCATATCCACACCTTTGCGCATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGGTATGGGTGGAAAACAGTTTGCCAGTAAACAGCAGCCCCAAGGGAGGGGCATGAACAAACAGAGACGGGCTAAGACCACTGAGAGACTGCCAGT‑TACCGTTGATGGttttagagcgaatt                                                                                                                                                                                >  2:323568/1‑241 (MQ=255)                                                           gcgATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATACGGTTCAACGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTCAACAGCAGCACCCAGGTGGCGGCATTAACTAACAGTGAAGGGCCACAACCAATGGGTGACTGCCCGTGTAC‑GTTACTGATTTTGGCGCGAATTAAcac                                                                                                                                                                           >  2:559442/1‑249 (MQ=255)                                                                     cGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTg                                                                                                                                                                 <  1:330825/251‑1 (MQ=255)                                                                          gTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCttt                                                                                                                                                            <  2:70632/251‑1 (MQ=255)                                                                               gAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCgg                                                                                                                                                       >  1:151054/1‑251 (MQ=255)                                                                                    gCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATa                                                                                                                                                  <  1:433145/251‑1 (MQ=255)                                                                                       cACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATAc                                                                                                                                                 <  1:361720/249‑1 (MQ=255)                                                                                       cACCCAACAGGCGCTTCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTGCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGTAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCTAATTAACGCGTCGCTGGTGGCTTTGTCGGGAATAc                                                                                                                                                 <  2:551859/249‑1 (MQ=255)                                                                                                ggCGCTTCGCCATATCCTCACCTTTGCTGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCt                                                                                                                                      <  2:471479/251‑1 (MQ=255)                                                                                                   gcATCGCCATATCCACACCTTGGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTcc                                                                                                                                   <  1:130075/251‑1 (MQ=255)                                                                                                               ccACACCTTGGCGGATCAGCGGTTCCCCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATc                                                                                                                       <  2:151054/251‑1 (MQ=255)                                                                                                                 ccaccgtgtcGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGGAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAg                                                                                                                     <  2:285880/242‑1 (MQ=255)                                                                                                                           ggATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGCCCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACACGTCCCGGGTGGCTTTGTCGGGAATACGCAACACCGCTTCCGCCAGACACATCCGCGcccagcct                                                                                                           >  2:505278/1‑245 (MQ=255)                                                                                                                           ggATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCt                                                                                                           >  1:395000/1‑251 (MQ=255)                                                                                                                             aTCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTc                                                                                                         <  1:191771/251‑1 (MQ=255)                                                                                                                               cAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAAAGCTTCCGCCAGACACATCAGCGCCACGCCTTccc                                                                                                       >  2:600067/1‑250 (MQ=255)                                                                                                                               cAGCGGTTACCCGCTTTTACCGATAAGACGGGTCAGCGAGCGGGGGAGGCTTGCTTCGTTAGGGGTGGAAACCAGTATGCCAGGAAACAGCATCCCCCATGTGGCGGCATTAACTAACAGGGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGGCtttat                                                                                                       <  2:571829/251‑4 (MQ=255)                                                                                                                               cAGCGGGTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCt                                                                                                       <  2:395000/251‑1 (MQ=255)                                                                                                                                     ttCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATg                                                                                                 >  1:31873/1‑251 (MQ=255)                                                                                                                                          cGCTTTAACCGATAATGCGGTTCAGAGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGACAGTAAAAAGCAGCCCCCAAGTGGCGGAATTAACTAACAGTAACGGGCTACGACCAATGTATGACTGCCAGT‑TACCGTTGCTGATTTTGTCGAGACTTAACGCATCGCGGGTGGCTTTGTCGGGAATAAGCACCAACGCTTCCACCAGACCCATCCACGCCACGCCTTccgcccatggcacc                                                                                            >  2:582638/1‑239 (MQ=255)                                                                                                                                                     atcaTGCGATTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGGTTGCCAGTAACCAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCCATGTGTGACTGCAAGT‑TAACGTTGCTGATTTTGTCTCGAATTAACGCGTCGCGGGTGGCTTTTTCGGGAATACGCAACAAAGCTTCCGCCCTCCACCTCACCGCCACGACTTCCAGCGGTAACAGC‑AAAAcagcctg                                                                                >  2:364119/4‑245 (MQ=255)                                                                                                                                                     aTAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGACACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAAATAACAGTGCCGGGCTAAGACCAATGTGTGAATGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGGGGCTTTGTCGGGGATACGCACCAAAGGCTCCGCCAGGCACACCAGCGCCAAGCCTTCCTGCGACGACAGCGCAAACTCCtc                                                                                 >  2:525010/1‑250 (MQ=255)                                                                                                                                                       aaggcggtgccgcgagcgGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGGAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTCCGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTga                                                                                <  2:518666/239‑2 (MQ=255)                                                                                                                                                          gCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATa                                                                            <  1:469291/251‑1 (MQ=255)                                                                                                                                                           cGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATa                                                                            >  1:38737/1‑250 (MQ=255)                                                                                                                                                                         ggAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATAc                                                             >  1:459722/1‑251 (MQ=255)                                                                                                                                                                               gctggctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTgcgc                                                       >  1:284649/1‑251 (MQ=255)                                                                                                                                                                                  ggctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGAGAc                                                    >  1:23697/1‑251 (MQ=255)                                                                                                                                                                                           aTGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCa                                            <  1:114552/250‑1 (MQ=255)                                                                                                                                                                                           aTGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCa                                            <  1:416709/250‑1 (MQ=255)                                                                                                                                                                                           aTGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCAt                                           <  1:71767/251‑1 (MQ=255)                                                                                                                                                                                                 ggAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCAtttttt                                      <  1:340336/250‑1 (MQ=255)                                                                                                                                                                                                       cAGTTTGCCCGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGAtgac                               <  2:57018/251‑4 (MQ=255)                                                                                                                                                                                                                aGTAAACAGCACCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGCCTGCCAGT‑TACCGATGCTGATTTTGTCGCGAATTAACGCGTCGCGGGGGGCTTTGTCGGGAATACGCAACAACGCATCAGCCAGCCACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGAAGTTTAt                      >  2:54668/1‑251 (MQ=255)                                                                                                                                                                                                                   aaaCAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCg                    <  2:208203/250‑1 (MQ=255)                                                                                                                                                                                                                       agcagcCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCa                                            <  1:327641/222‑1 (MQ=255)                                                                                                                                                                                                                       agcagcCCCCAGGCGGCGGCATTAACTAACAGCGCCGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGAAACCTCAGCGCCACGCCTTCCTGCGATGCAAGCGAAAACTCCTGCACTAACCCCTGGCCCATACATGCGCGACCACTGGCa                                            >  2:327641/1‑222 (MQ=255)                                                                                                                                                                                                                            cccccAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGAt          <  1:672897/251‑1 (MQ=255)                                                                                                                                                                                                                             ccccCGGTGGCGGCATTAACTAACAGTGACGGTCTACGACCATTGGGTGACTGCCGGT‑TACCGTTGTTGATTGTGTCGCGAATTAACGCGGCGCGGGTGGCGTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCAGACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATa         <  2:38737/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggcggcATTAACGAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCTGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACATCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGt  <  1:92440/250‑1 (MQ=255)                                                                                                                                                                                                                                      gcggcATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCa                <  1:237228/235‑1 (MQ=255)                                                                                                                                                                                                                                      gcggcATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCa                >  2:237228/1‑235 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               ATCGTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGT‑TACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGT  >  NC_000913/1078028‑1078505

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