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	98.0%	145.7 / ‑3.9	49	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 (0/1);  new base T (22/26);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

GTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTG  >  NC_000913/1078031‑1078510                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   gTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGATAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAAACCGTTTGCCAGTTAACCGCAGCCCCCAGGTGGCGGCATTAACTACCCGTGACGGGc                                                                                                                                                                                                                                       >  2:26942/1‑251 (MQ=255)  aaaGAGTAACGGAAACCTGTCCCTTCCAGCTTGCGGTCATTGGCTAACGCTTCCGCGATTGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCAGCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGTGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCt                                                                                                                                                                                                                                      <  2:139574/250‑1 (MQ=255)    agagtacaggAAACCTTTCTCGCCCAGCTTGCGGGCATTGGCTACCGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTAc                                                                                                                                                                                                                                    <  2:392087/243‑1 (MQ=255)                         ttCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGAGACGGGCTACGACCAAtgtgtacatgcccgg                                                                                                                                                                                                               >  2:458418/1‑241 (MQ=255)                           ccAGCTTGCGGGCTTTGGATAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCTCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTa                                                                                                                                                                                                             <  1:330623/251‑1 (MQ=255)                                                  gCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGGTTTTGTCGCGAAt                                                                                                                                                                                      >  2:135433/1‑251 (MQ=255)                                                     tCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTa                                                                                                                                                                                    <  1:418430/250‑1 (MQ=255)                                                       cgcgATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAAcg                                                                                                                                                                                 <  1:444000/251‑1 (MQ=255)                                                              gTTTCGCCAGGGACGAAATGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTTACAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCCATGTATGACTGCCCGTTACCGTTGCTGATTTTGTCGCGAATTAccgcgccgcg                                                                                                                                                                          >  2:417548/1‑251 (MQ=255)                                                                  cGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTg                                                                                                                                                                      <  1:10608/251‑1 (MQ=255)                                                                   gCCAGTGACGAACTGCTCACCCAGCAGGCGCCTCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTgg                                                                                                                                                                     <  2:137038/251‑1 (MQ=255)                                                                      aGTGACGAACTGCTCCCCCATCAGGCGCATCGCCATATCCACACCTGTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCtt                                                                                                                                                                  <  2:399328/251‑1 (MQ=255)                                                                         gACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGt                                                                                                                                                               <  1:20370/251‑1 (MQ=255)                                                                                   ccACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACg                                                                                                                                                     <  1:428654/250‑1 (MQ=255)                                                                                                    cGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCa                                                                                                                                    <  1:250510/251‑1 (MQ=255)                                                                                                           tCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGG‑‑GGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCCACAACGCTTCCGCCAGACCCATcc                                                                                                                           >  2:75215/1‑250 (MQ=255)                                                                                                             cacaCCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCa                                                                                                                           <  1:175859/251‑1 (MQ=255)                                                                                                                      gCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACg                                                                                                                   <  1:421473/250‑1 (MQ=255)                                                                                                                          aTCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGCGCGGGAGAGGCTGGCTTCGTTTTGGGTGGAAACCAGTTTTCCCGTAAACAGCAGCCCCCAGGGGGCGGCATTAACTAACAGTGACGGGCTACGCACACTGTGTGACTGCCAGTTACCGTTGCTGATTTTTTTGCGAATTAACTCGTCGCGGGTGGCTTTGTCGGGCATACGCAACAACGCTTCCGCCAGacaaccaagcgccccgccttc                                                                                                              >  2:36738/1‑233 (MQ=255)                                                                                                                               gggTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCGTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCg                                                                                                         <  2:357815/250‑1 (MQ=255)                                                                                                                                    cACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCGGCGCCACGCCTTCCTGCGagtac                                                                                                    >  2:229030/1‑247 (MQ=255)                                                                                                                                              aCCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGG‑‑GGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTcc                                                                                        <  1:75215/251‑1 (MQ=255)                                                                                                                                               ccGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTc                                                                                         >  1:75703/1‑251 (MQ=255)                                                                                                                                                    aaTGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCCTTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCCCCACGCCTTCCTGCGATGCCAGCGAAAACTCCTGc                                                                                     >  2:371631/1‑250 (MQ=255)                                                                                                                                                     aTGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCa                                                                                    >  1:453986/1‑250 (MQ=255)                                                                                                                                                      tGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAAt                                                                                  >  1:91552/1‑251 (MQ=255)                                                                                                                                                       gCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATa                                                                                 <  1:327856/251‑1 (MQ=255)                                                                                                                                                               gcgagcgGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGTTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTgg                                                                         <  2:242423/251‑1 (MQ=255)                                                                                                                                                                     gggAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGGTGACAGCGAAAACTCCTGCAATAACCCCTGGACCAtt                                                                   >  2:216573/1‑250 (MQ=255)                                                                                                                                                                      ggAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGGGTCGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCATACTCATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATAc                                                                  <  2:69909/251‑1 (MQ=255)                                                                                                                                                                         gtggttggctTCGTTATGGGTGGAAACCAGTTGGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCGGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCGGGACCATACCTg                                                               <  2:91552/246‑1 (MQ=255)                                                                                                                                                                               ggctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGGGGAGGCATTAACTAACAGTGACGGGCTACGACCAATGGGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGAc                                                         >  2:310952/1‑251 (MQ=255)                                                                                                                                                                                gctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGAcc                                                        <  1:145159/251‑1 (MQ=255)                                                                                                                                                                                 ctTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGAccc                                                       >  2:13796/1‑250 (MQ=255)                                                                                                                                                                                   tCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCa                                                       <  1:135433/249‑1 (MQ=255)                                                                                                                                                                                    cGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTg                                                    <  1:458418/251‑1 (MQ=255)                                                                                                                                                                                        aTGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCAt                                                <  1:455034/251‑1 (MQ=255)                                                                                                                                                                                          gtgTGGAAACCAGGTTGCCAGTCAACCGCAGCCCCCAGGTGGAGGCATTAACTAACAGTGTCGGGCTACGACCAATGTGTTACTGCCAGTTGCCGTTGCTGATTGTGTCGCGAATTAACGCGCCGCGGGTGGCTTTGTCGGGAATTCGCAACAACGCTGCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCTTGCAATAACCCCTTGACCATACCTTCGCGACCACTGGCAttt                                              <  2:454987/249‑1 (MQ=255)                                                                                                                                                                                           ggTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCAtttt                                             >  1:166502/1‑251 (MQ=255)                                                                                                                                                                                               gAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCAtttttttt                                         >  1:139187/1‑250 (MQ=255)                                                                                                                                                                                                    cAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTAc                                    >  1:94812/1‑251 (MQ=255)                                                                                                                                                                                                              gTAAACAGCAGCCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATACGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGAACTTACCTGCGCGACCACTGGCCTTTTTTTGTTTACGCAGTTTATc                          >  2:123970/1‑251 (MQ=255)                                                                                                                                                                                                                    agcagcCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAg                    <  2:70998/251‑1 (MQ=255)                                                                                                                                                                                                                    agcagcCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAg                    <  1:165881/251‑1 (MQ=255)                                                                                                                                                                                                                        gcCCCCAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTg                 >  1:345175/1‑250 (MQ=255)                                                                                                                                                                                                                             cAGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGAAGTTTATCGGCCAGCTGAAACGc           >  2:412906/1‑251 (MQ=255)                                                                                                                                                                                                                              aGGTGGCGGCATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGcc          <  2:94812/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggcATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCAGCCAGCTGATACGCTAGTTtgg   >  1:136110/1‑250 (MQ=255)                                                                                                                                                                                                                                      gcATTAACTAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCCGTTtgtg  >  2:205358/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   GTAAGAGTAACGGAAACCTTTCTCTTCCAGCTTGCGGGCATTGGCTAACGCTTCCGCGATGGTTTCGCCAGTGACGAACTGCTCACCCATCAGGCGCATCGCCATATCCACACCTTTGCGGATCAGCGGTTCACCGCTTTTACCGATAATGCGGTTCAGCGAGCGGGAGAGGCTGGCTTCGTTATGGGTGGAAACCAGTTTGCCAGTAAACAGCAGCCCCCAGGTGGCGGCATTAACAAACAGTGACGGGCTACGACCAATGTGTGACTGCCAGTTACCGTTGCTGATTTTGTCGCGAATTAACGCGTCGCGGGTGGCTTTGTCGGGAATACGCAACAACGCTTCCGCCAGACACATCAGCGCCACGCCTTCCTGCGATGACAGCGAAAACTCCTGCAATAACCCCTGGACCATACCTGCGCGACCACTGGCATTTTTTTGATTACGCAGTTTATCGGCCAGCTGATACGCCAGTTTGTG  >  NC_000913/1078031‑1078510

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