Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,625	G→A	S70F (TCC→TTC)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,625	0	G	A	98.3%	190.0 / NA	59	S70F (TCC→TTC)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/37);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGA  >  NC_000913/3176381‑3176870                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      cAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGAt                                                                                                                                                                                                                                                 >  1:325742/1‑251 (MQ=255)  aGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATc                                                                                                                                                                                                                                                >  1:57539/1‑251 (MQ=255)     atacactcccgcccatCAAGCTGAAACGCGCGCAGGTCACCGTGCGGCACGCCAAACACCGTGCTATCCAGCAACAGGATTTGCCCTTGCTCACCAATAAACACGCTCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGTAAATC‑GTGGTTGCCCGGCAGCCAGACTCAGGGCGCACGAAAACTTGCGATGCCTTTAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCt                                                                                                                                                                                                                                              <  2:355218/235‑1 (MQ=255)       aCGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTgc                                                                                                                                                                                                                                            <  1:430000/250‑1 (MQ=255)           tCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTa                                                                                                                                                                                                                                        <  1:484730/250‑1 (MQ=255)            ccAGCCCCTCAATCAGAAACTCGCCCAGCTCACCGGGCGGCACGCCAAACCCCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGTGCTGGAAAGC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTaaa                                                                                                                                                                                                                                      <  2:374141/251‑1 (MQ=255)            ccAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTaaa                                                                                                                                                                                                                                      <  1:34902/251‑1 (MQ=255)             cAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGAGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAt                                                                                                                                                                                                                                     >  1:449280/1‑251 (MQ=255)              aGCCAATCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTAGCCAGCAACAGGATTTGCCATTGCTCACCAATAATCACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATc                                                                                                                                                                                                                                    <  2:321459/251‑1 (MQ=255)                     cAAGCTGAACCTCGCGTAGCTCACCGGGCGGCCCGCCACACTGCAGGCTATCCAGCAACAGGATTTGCCAGTGCTCAACAATAAACACGCGGTGCGCGGGGGATATACCCGCATCCTGTAACGCTAGGGTCACCGCGGGCTGGAAATC‑GTGGTGGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCGGATAGGCCGGCGAGAATTGATCCTGCGCTAAATCACCTGta                                                                                                                                                                                                                             <  2:303570/251‑2 (MQ=255)                     cAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTg                                                                                                                                                                                                                               <  1:120999/249‑1 (MQ=255)                     cAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTc                                                                                                                                                                                                                             <  2:130637/251‑1 (MQ=255)                     cAAGCGGAAACTCGCGCAGCTCACCGGGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCGTCGCCGGGGAGATACCCGCAGCCTGTAACGCGCTGTACAGCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTg                                                                                                                                                                                                                               <  2:427566/249‑1 (MQ=255)                           gAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGAc                                                                                                                                                                                                                        >  1:24978/1‑250 (MQ=255)                           gAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACa                                                                                                                                                                                                                       <  1:241629/251‑1 (MQ=255)                                         tcaCCGTGCGGCGCGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAACTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGACGCGACAATCAGGACGAATTc                                                                                                                                                                                                         >  2:137656/1‑251 (MQ=255)                                          caCCGTGCGGCCCGCCAAACCCCTGGCTATCCAGCAACAGGATTTGCCAGTGCTCACCAATAAACCCGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTTGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCg                                                                                                                                                                                                        <  2:325742/251‑1 (MQ=255)                                            ccGTGCGGCACGCCAAACGCCTGGGTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCGTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTg                                                                                                                                                                                                      <  1:107467/251‑1 (MQ=255)                                            ccGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTg                                                                                                                                                                                                      <  1:267238/251‑1 (MQ=255)                                                        ccAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTcg                                                                                                                                                                                                                            >  2:125787/1‑217 (MQ=255)                                                        ccAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTcg                                                                                                                                                                                                                            <  1:125787/217‑1 (MQ=255)                                                            actcctcgCTATCCAGCAACAGGATTGGCCATTGCTCACCAATAAACACGCGCGTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATc                                                                                                                                                                                      <  2:101011/244‑1 (MQ=255)                                                                   gCTATCCAGCAACAGGATTTGCAATTGCTACCCAATAAACACGCGCTTCTCCGGGGAGATTCCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAATCTTGCGATGCCTGCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGGGCTGGTGTGGCCGAATCGCCTCCa                                                                                                                                                                               <  2:372262/251‑1 (MQ=255)                                                                            cAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCt                                                                                                                                                                      >  1:427275/1‑251 (MQ=255)                                                                            cAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCt                                                                                                                                                                      <  1:474101/251‑1 (MQ=255)                                                                                      tGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAActct                                                                                                                                                             <  1:102528/250‑1 (MQ=255)                                                                                            tGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGAACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGCCAATCAGGTCGATTTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAggg                                                                                                                                                      >  2:451619/1‑250 (MQ=255)                                                                                                       aaaCACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGTCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTGCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGGTTACCCCTa                                                                                                                                           <  2:295651/251‑1 (MQ=255)                                                                                                            cgcgCTTCGCCGGGGAGATACCCGCATACCGTAACGCGCTGTACAACGCTGGCTGGAAACC‑GTGGTTGCCCGGAAGCCAGACGGAGGGCGAAAGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGGCGCAGAGAAGTGATCCTGCGCTAAAGCACCGGGCGCGAAAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGCTTACCCCTAACAgg                                                                                                                                      <  1:361185/251‑1 (MQ=255)                                                                                                             gcgcTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCCCCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAggg                                                                                                                                     <  1:19184/251‑1 (MQ=255)                                                                                                                   cGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCg                                                                                                                                <  1:406249/250‑1 (MQ=255)                                                                                                                       ggggAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGATTACCCCTAACAGGGCTTCGTGCtt                                                                                                                           >  1:387308/1‑251 (MQ=255)                                                                                                                        gggAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCttt                                                                                                                          >  1:76411/1‑251 (MQ=255)                                                                                                                        gggAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCttt                                                                                                                          <  1:214340/251‑1 (MQ=255)                                                                                                                           agaTACCCGCCTCCTGTAACGCGCTGTACAACGCGGGCTGGAAATC‑GTGGTGGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGACGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGGTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                                       <  2:146371/251‑1 (MQ=255)                                                                                                                             aTACCCGCATCCTGTAACGCGTTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACGTGCGATGCCGTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAAGCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                                     <  2:480652/251‑1 (MQ=255)                                                                                                                             aTACCCGCATCCTGTAACGCGCTGTACAGCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                                     <  1:224443/251‑1 (MQ=255)                                                                                                                                 ccGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAAc                                                                                                                 <  1:451619/251‑1 (MQ=255)                                                                                                                                            aaCGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTc                                                                                                      <  1:192800/251‑1 (MQ=255)                                                                                                                                              ggcggttgtcATCGCGGGCTG‑TAATCTGTGGTTGGCCGGCGGCCAGGCGAAGGGCGATCGAAAACTGGCGATGTCTTCAGCGAAATACGGATAGGGCGCAGAGAATTGAGTCTGCGCTAAATCACCGGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGGGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGGGTCgg                                                                                                    <  2:428823/242‑1 (MQ=255)                                                                                                                                                  cTGTACATCGCGGGCAGGAAATC‑GTGGTTGCCCGGCAGCCAGCCGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAAATCGTGCTGGGGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTGTACCCCTACCAGGGCTTCGTGCTTTTGTGAAAACAGGAGAGTGTCGGGAAt                                                                                                >  2:571739/1‑251 (MQ=255)                                                                                                                                                   tGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTACATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                                     >  2:49475/1‑229 (MQ=255)                                                                                                                                                   tGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGc                                                                                                                     <  1:49475/229‑1 (MQ=255)                                                                                                                                                      aCATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGGAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGt                                                                                            >  2:560124/1‑251 (MQ=255)                                                                                                                                                        aTCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATGGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTa                                                                                           >  1:420255/1‑250 (MQ=255)                                                                                                                                                          cgcgGGCTGGAAAGC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGGGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                         <  2:427275/250‑1 (MQ=255)                                                                                                                                                                cTGGAACTC‑GTGGTTGCCCGGCAGCCAGACGCAGGGAGCACGAAAACTAGCGATGCCTTCAGCGACATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCATCAAGAACCGCCTGGTAACTCTCCCCGGTGTTTAGCCCTACCAAGGCTTAGTGCATTTGTGCAACCAGGTGAGAGTAGGTAAtttagaaaaacctg                                                                                   >  2:495300/1‑239 (MQ=255)                                                                                                                                                                             gttccccGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAc                                                                      <  2:404112/247‑1 (MQ=255)                                                                                                                                                                                      cAGCCAGACGCAGGGCGCGCGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGACGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGGGTTTACCCCTAACAGGGCTTCGTGCTTTTGGGCCAACAGGTGAGTGTCGGTAATTTATAAAATCATGCTTCTTGCCTCCccagcaca                                                               >  2:347458/1‑246 (MQ=255)                                                                                                                                                                                            gACGCAGGGCGCACGCAAACTTGCGATGCCTTCAGAGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGATCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAggg                                                       >  1:435868/1‑251 (MQ=255)                                                                                                                                                                                                cAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTa                                                    >  1:425061/1‑250 (MQ=255)                                                                                                                                                                                                           caaaCTTGCGTTGCCTTCAGCTAAATGCTTATAGGCCCCATAGAATTGATCCTGCGCTAAATCTCCTGTCGCGACAATCAGTGCTAAGTCGTGCTGGTGTGGTCAAATCGCCTCCAGCCCCGCCTGGTATCCCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGGGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTc                                          <  2:1283/248‑1 (MQ=255)                                                                                                                                                                                                                    gATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTc                               <  1:93271/251‑1 (MQ=255)                                                                                                                                                                                                                       gCCTTCAGCGAAAGGATGATAGGCTGAAGGTAATTGATACTGCGCTAAAACGCCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTATCTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAAGCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCtg                            <  2:449280/251‑2 (MQ=255)                                                                                                                                                                                                                                  aaTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacg                  >  1:268694/1‑250 (MQ=255)                                                                                                                                                                                                                                       tGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGACATTTGTAAAATCCTGACTCTGGCCTCAACAGCCAGAGGAAGGGTTAACAGGCTATCCAAATGGTGTCCTTAGGTTTCacgccgat             >  2:13924/1‑245 (MQ=255)                                                                                                                                                                                                                                         aTAGGCCGCAGCGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACACCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAACTCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCATTCCAAATTGAGTCCTTAGGTTTCACGACGCtaat          >  2:536197/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccAGAGAATGGATCCTTCGCCATATCACCTGTCCCGAAAATCAGTGCGAATTCGTGTTGGTTTGTCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACATGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGGCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGTGTTAACAGGCTTTCCAGATGGTGTCCTTAGGTTTCACGACGCTAATAAACCgg   <  2:201049/250‑1 (MQ=255)                                                                                                                                                                                                                                                 cAGAGAATTGAGCCTGCGCTAAATCACCGGTCGCGACAATCAGGGCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGACTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGa  <  1:536197/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATC‑GTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCGGA  >  NC_000913/3176381‑3176870

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