Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,570	A→C	intergenic (+538/‑64)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,570	0	A	C	98.5%	209.1 / NA	66	intergenic (+538/‑64)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (33/32);  total (34/32)
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.63e-01

CTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATT  >  NC_000913/89348‑89812                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     cTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAAAGTTAATTTGATGCGCACGTGGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGAATGGCTTGACAAGCTTTTCCTCCCCTCCGTAAACTCCTTTCAGTGGGGAAt                                                                                                                                                                                                                          >  2:553947/1‑251 (MQ=255)      aTGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTgg                                                                                                                                                                                                                     >  1:100247/1‑251 (MQ=255)                cgTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTgg                                                                                                                                                                                                           >  1:198872/1‑251 (MQ=255)                   cGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGa                                                                                                                                                                                                         <  2:198872/250‑1 (MQ=255)                                    tatTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGGAATTGTGGGGCAAAGTGGGAATAAGGGGGTAGGCTggg                                                                                                                                                                                       >  2:636127/1‑250 (MQ=255)                                        tgatttgGAAAATAGTGAGTCATTTTAAAGCGGTGATGACGATGCGGGATTTTCTCTTTCAGTTATTCATAACGATAATTTTCTGCGCACGTTGGACGTAAAATAAACAACGCTGATAGTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTGGATAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGTGTAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATg                                                                                                                                                                                    <  2:2213/245‑1 (MQ=255)                                         gCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGtt                                                                                                                                                                                  <  2:363228/251‑1 (MQ=255)                                                         aGTCATTTTAAAACGGTGAGGTCGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGTTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCGCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCGAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTa                                                                                                                                                                  <  1:277474/251‑1 (MQ=255)                                                           tCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGaa                                                                                                                                                                                                        >  1:267290/1‑211 (MQ=255)                                                           tCATTTTAAAACGGTGAGGATGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGaa                                                                                                                                                                                                        <  2:267290/211‑1 (MQ=255)                                                            cATTGTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCTGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCGTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGGGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTc                                                                                                                                                               <  1:395060/251‑1 (MQ=255)                                                                ttAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAAt                                                                                                                                                            <  1:332361/250‑1 (MQ=255)                                                                 tAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGTCGTAAAATAAACAACGCTGATATTAGCCGTAATCATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGTCAAGCTTTTCCTCCGCTCCGTAAACGCCGTTCAGTGGGAAAGTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCAGGTTCCGGGGAGCAACGTTAGTCAAtc                                                                                                                                                           <  2:605733/250‑1 (MQ=255)                                                                        gTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGc                                                                                                                                                   <  1:432902/251‑1 (MQ=255)                                                                           atgaCGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCaa                                                                                                                                                 >  1:39414/1‑250 (MQ=255)                                                                              aCGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAggg                                                                                                                                             <  1:124246/251‑1 (MQ=255)                                                                                gATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTACGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGGGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGGGCAACGGTAGTCAATCTCGTAAGCAAAggggg                                                                                                                                           >  2:378222/1‑249 (MQ=255)                                                                                       gATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑Ca                                                                                                                                   <  1:669736/251‑1 (MQ=255)                                                                                              tCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGTAATTTTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCtt                                                                                                                            >  1:34511/1‑250 (MQ=255)                                                                                                 tACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGCTATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCCAAGTGGGAATAAGGGGTGTGGATGGGATGTTCCGGGGAGCAACGATAGTCAATCTCGACAGCAAAGGGCGCT‑ATACAGTGCCAAccc                                                                                                                         >  2:169423/1‑251 (MQ=255)                                                                                                  aCAGCTATTCCTAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGCGGGAATACGGGGTGAGGCTGGACTGTTCCGGGGGGCAACGTTAGTCACACTCGACAGCAAAGGGCGCGTAT‑CAGGGCCtccagg                                                                                                                        >  2:329799/1‑246 (MQ=255)                                                                                                  aCAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATACACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCCGTGGGGAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCGGGCATGTTCCGGGGAGCAACATTAGTCAATCTCGACAGCAAAGGGCGATTAT‑CCGTGCCTAcccc                                                                                                                        >  2:348869/1‑250 (MQ=255)                                                                                                  aCAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTAc                                                                                                                           <  1:329799/248‑1 (MQ=255)                                                                                                  aCAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCg                                                                                                                        <  1:281986/251‑1 (MQ=255)                                                                                                    aGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGAAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCGTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGtt                                                                                                                      <  2:367687/251‑1 (MQ=255)                                                                                                        aTTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTTGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTaccg                                                                                                                  >  1:626130/1‑248 (MQ=255)                                                                                                                cGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAgc                                                                                                          <  1:360482/251‑1 (MQ=255)                                                                                                                    aaTTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGTAAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCCAAGTGGGAATAAGGGGGGAGGCTGGCATGTTCCGGGGGGCCACGTTAGTAAATCTTGGAAGCAAAGGGGGCTTAT‑AAGTGCCTACCCGTTACCGGGGACAtctgct                                                                                                      >  2:519408/1‑251 (MQ=255)                                                                                                                      tttgctgagCACGTTGGACGTAAAATAAACACCGCTGTTATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAACTTTTTCTCCGCTCCGTCACCTCCTTTCAGTGGGAAATTATGGTGCAAAGTGTGAATCAGGGGTGTGGCTGGCATATTTCGGGGAGCAACGTATGTCAATCCCGACAGCACAGGGCGATTTA‑CAGTGCCTACCCGTTACCGGGAActaacgcttg                                                                                                    >  2:424452/9‑242 (MQ=255)                                                                                                                       ttGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTg                                                                                                    <  2:592363/250‑1 (MQ=255)                                                                                                                        tgttgcGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCGCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTga                                                                                                   <  2:521643/245‑1 (MQ=255)                                                                                                                         gctgcGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTAGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGACATTGGGGGGCAAAGGGGGAATAAGGGGTGTGGCTGGCATGTTCCGGCGAGCAACGGTAGTCAATCTTGACAGCAAACGGAGCTTAT‑CGGTGCATACCCGTTaccgtgcaccactgcttgatc                                                                                                 >  2:77531/5‑231 (MQ=255)                                                                                                                          ctgcGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTgaga                                                                                                 >  1:589785/4‑250 (MQ=255)                                                                                                                               cACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCAt                                                                                                                                                                                     <  2:139801/162‑1 (MQ=255)                                                                                                                               cACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCAt                                                                                                                                                                                     >  1:139801/1‑162 (MQ=255)                                                                                                                                aCGTTGGACGTAAAATAAACAACGCTTATATTAGCCGTAAACATCGGGCTTTTTACCTAGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTCAACTACTTTCAGTGGGAAATTGTGGGGCAAAGGGGGAATCAGGGGTGAGGCTGGCATGTTCCGGGGAGCACAGATAGACTATCACGGCAGCAAAGGGGGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAActcctc                                                                                          >  2:436400/1‑246 (MQ=255)                                                                                                                                 cGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGc‑c                                                                                        <  1:553199/251‑1 (MQ=255)                                                                                                                                  gTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑Cg                                                                                       >  1:286949/1‑251 (MQ=255)                                                                                                                                       aCGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCaa                                                                                  <  1:328901/251‑1 (MQ=255)                                                                                                                                       aCGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCa                                                                                   >  1:601246/1‑250 (MQ=255)                                                                                                                                           aaaaTAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCCGTGGGAAATTGTGGGGCAAAGTGGGAATACGGGGTGAGGCTGGCATGGTCCGGGGAGCAACGTTAGTCAATCTCGACAGCACAGGGCGATTAT‑CAGTGCCTACCCGTTAACGGGAACAGCTGCTTGAGACCGCTGC‑CGGACAAATgg                                                                              >  1:51446/1‑251 (MQ=255)                                                                                                                                           aaaaTAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATgg                                                                              <  1:82690/251‑1 (MQ=255)                                                                                                                                           aaaaTAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCACGGAATGCCTTGTGACTGGCATGACAAGCTTTTCTTCCGCTCCGTAAAATCCTTTCAGTGGGAAATTGTGGGGCAAAGGGGGAATAAGGGGGGAGGCTGGCAGGTTCCGGGGAGCAACGTTAGTCAATCTAGACGGCAAAGGGCGGTTAT‑CAGGGCCTAAACGTTACCGGGAACAGCGGCTTGGGAACGC‑GCACGGTCaaaggg                                                                              >  2:293548/1‑248 (MQ=255)                                                                                                                                              aTAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGttt                                                                           >  1:106116/1‑251 (MQ=255)                                                                                                                                                aaaCAACGCTGAGTTTAGCCGTAAACATCGGGTTTTTTTCCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTGCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCGCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGc                                                                         <  2:464507/251‑1 (MQ=255)                                                                                                                                                        cTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTg                                                                  >  1:267712/1‑250 (MQ=255)                                                                                                                                                         gtatatTAGTCTTAAACAGCGGTTTTTTTACCTCGGTATGCCTTGTGACTTGTTTGACATGCTTTTCCACCGCGCCGTAAACTCCTTGCAGTGGGAAAGTGGTGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCGCGACAGCAAAGGTCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTACTTGAGAACGCTGC‑CGGTCAAATTGTTTGCACCAtttac                                                                <  2:280303/249‑4 (MQ=255)                                                                                                                                                                aGCCGTAAACATCGGGATTTTTACCTCGGTATGCCTTGAGACTGGCTTGACAAGCTTTTTATCAGCTCCGTAAACTCCTTTCAGTGGGAACTTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTAAAACTCGAAAGCAAGGGGCGCTTAT‑CAGTGCCTATCCGTTACCGGGCACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTgacccattatcattcatc                                                         >  2:586146/1‑234 (MQ=255)                                                                                                                                                                    gTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGGGGGGCAAAGTGGGCATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGAAAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCACCAccc                                                     >  2:259898/1‑251 (MQ=255)                                                                                                                                                                           tCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCctg                                              >  1:63195/1‑251 (MQ=255)                                                                                                                                                                                   tttACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCttt                                       <  1:362749/250‑1 (MQ=255)                                                                                                                                                                                    ttACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAAGTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTa                                      <  1:665496/250‑1 (MQ=255)                                                                                                                                                                                     tACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTAc                                     <  1:169423/250‑1 (MQ=255)                                                                                                                                                                                      aCCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTAccc                                   >  1:258255/1‑251 (MQ=255)                                                                                                                                                                                       ccTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCCCCGTAAACTCCTTTCAGTGGGGAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTAcccc                                  >  2:173483/1‑251 (MQ=255)                                                                                                                                                                                           ggTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCcctgc                              <  1:519408/251‑1 (MQ=255)                                                                                                                                                                                            gTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCCCCGTAAACTCCTTTCAGTGGGAGATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCCAAGGGCGCTTAT‑CAATGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACCTTCATCACCCGGGCCTCATGCTTTACCCcccgcc                             >  2:510429/1‑251 (MQ=255)                                                                                                                                                                                             gaTGCCTTGTCGCTGGCCTGAAAAGCTTTTTCTACGCGCCCTAAACTCCTTTCAGTGGGAAATTTTGGGGCAAAGGGGGGATAAGGGGTGAAGCTGGCATGTTCCGGGGAGCAACGTTAGGCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCTGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATTGTCTGCACCATTGACATTCATCACCCGTTCCTGCTGTTTTACCCcctgcct                            <  2:336939/250‑1 (MQ=255)                                                                                                                                                                                              aTGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTc                                                                                    >  2:216827/1‑194 (MQ=255)                                                                                                                                                                                              aTGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTc                                                                                    <  1:216827/194‑1 (MQ=255)                                                                                                                                                                                                     gtgACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATggg                     <  2:63195/250‑1 (MQ=255)                                                                                                                                                                                                            gCTTGACAAGCTTGTCCCCCGCGCCTGAAGCTCGTTTCAGTGGGAAAGTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCGCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGGCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATc             <  2:106116/251‑1 (MQ=255)                                                                                                                                                                                                             cTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCg            >  1:44104/1‑251 (MQ=255)                                                                                                                                                                                                                gACAAGCTTTTCCTCCTCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAACGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCCAATGGTTTGCACCATTGCAATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGATTGGGAAATTATCGAGc         >  2:219076/1‑251 (MQ=255)                                                                                                                                                                                                                gACAAGCTTTTCCTCCGCTCCGTAAACTCCTTGCAGTGGGAAATTGGGGGGCAACGTGGAAATAAGGGGTGAGGCTTGCATGTTCCTGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTTCCTACCCTTTCCCGGGAACAGCTGCTTGAGAACTCTGC‑CGTTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTTCTTTACCCCCTGCCTGAATGGGAATTTATCGAGc         <  2:282088/251‑1 (MQ=255)                                                                                                                                                                                                                       ttttCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGATGCTTTACCCCCTGCGTGAATGGGAAATTATCGGGCAAAAAtt  >  2:474455/1‑251 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     CTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTAT‑CAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGC‑CGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATT  >  NC_000913/89348‑89812

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: ‑