Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,823,162	C→T	H255Y (CAC→TAC)	spoT →	bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,823,162	0	C	T	97.8%	153.9 / NA	46	H255Y (CAC→TAC)	spoT	bifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/17);  minor base G (0/1);  total (28/18)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA‑GATGGCGGAGATGGGT  >  NC_000913/3822920‑3823401                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  gTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATc                                                                                                                                                                                                                                           >  1:42665/1‑251 (MQ=255)                      aaCTCGAAGAGCTGGGTTTTAAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAACGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGAGGTCGCGCGAAGCATCTTTATTCGATTTACTGCAAAATGGGGCTCAAAGAGCAGCGTTTTTACTCGATCATGTACATCTACGCTTTCcgc                                                                                                                                                                                                                      >  2:438929/1‑251 (MQ=255)                                 cTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑Ca                                                                                                                                                                                                           >  1:372966/1‑250 (MQ=255)                                          gAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGAACCAGAAGATTCTTTCTGAAAGCGAAGGGCGTTTGCAGGAAGCGGGAATTCCGTGCCGCGTTAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATTGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTCCGCCTTCCGCGTGATCGT‑CAATGATTCTg                                                                                                                                                                                                  <  1:359486/250‑1 (MQ=255)                                                      cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATc                                                                                                                                                                                      <  1:156623/250‑1 (MQ=255)                                                            cGTTAGCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCt                                                                                                                                                                               <  2:372966/251‑1 (MQ=255)                                                                 tCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCCTCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                                                                                                                                                                          >  2:245032/1‑251 (MQ=255)                                                                     ggaaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTAGCGCGTGCTGGGcc                                                                                                                                                                          >  2:159546/3‑247 (MQ=255)                                                                     gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGAt                                                                                                                                                                      <  1:353191/251‑1 (MQ=255)                                                                      tAATCAAAGAAGGGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATg                                                                                                                                                                     >  1:141645/1‑251 (MQ=255)                                                                        agcAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTTAAGAGCTGATCCAGAAGATTCTTACTGAAAGCGAAGGGCGTTTGTATGAAGCGGGAATACCGTGCCGCTTCAGTGGTCGCGAGAAGCATCGTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTGTTACGCGATCATGGACTTCTGCGCTTTCCGCGTGAGCGT‑CAATGATTCGGACACCTGGTATCGCGTGCTGGGCCAGATGca                                                                                                                                                                   <  2:246857/249‑1 (MQ=255)                                                                            aagaagCGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGc                                                                                                                                                               >  1:117168/1‑251 (MQ=255)                                                                               aagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGCTGATCCAGAAGCTTCTTTCTGCAATCGACGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCACAATGGTGCTCAACGAGAAGCGTTTTTACTCGATCATGGACCTCTACGCTTTACGCGTGATCGT‑CAATGATTCTGACTACTGTTATCGCGTGCTGGGCCAGATGCACAGCCTg                                                                                                                                                            >  2:340539/1‑251 (MQ=255)                                                                                 gtattgAAATCCGCGCGCGTCAACCGTAAAGAGATGATCCAGAAGAGTCTTTCTGAAATCGAAGGGCGGTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTGCCGCGGGTTCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTa                                                                                                                                                          <  2:59248/247‑1 (MQ=255)                                                                                  tggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc                                                                                                                                                         >  1:406414/1‑251 (MQ=255)                                                                                       aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTACTGGGCCAGATGCACAGCCTGTACAAGcc                                                                                                                                                    <  1:226399/251‑1 (MQ=255)                                                                                                cgtggCAAACGTAAAGAGATGAACCAAAAGATTCTTCCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAGGGTGCTCAAAGAGCAGGGTTTTTACTCGATCATGGACATCTACGCTTGCCGCGTGATCGT‑CAATGATTCTGACCCCTGTTATCGCGTGCGGGGCCTGATGCACCGCCTGTACAAGCCGCGTCCggg                                                                                                                                           <  2:99343/248‑1 (MQ=255)                                                                                                     cAACAGTAAAGAGATGAGCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGATCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAAATGCACAGCCTGTACAAGCCGCGTCCGTGCcgcg                                                                                                                                      <  2:137701/251‑1 (MQ=255)                                                                                                            aaaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGCAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa                                                                                                                                  >  2:198307/1‑248 (MQ=255)                                                                                                               gagaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTAAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACta                                                                                                                            <  1:187654/251‑1 (MQ=255)                                                                                                                 gatgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGAGACCTGTTATCGCGTGCTGGTCCAGATGCACAGCGGGTACAAGCCGCGTCCGGGCCGCGTGAAATACtata                                                                                                                          <  1:385067/251‑1 (MQ=255)                                                                                                                   tgatcacgacacttctttctGACCTCGAAGGGCGAATGCCGAAAGCGGGCATACCGTCACGCGACAGTGGAAGCCAGCACCATATTTATACGATATACTGCAAAATGGTGCTCAAAGAGCAGCGTCTTTACTCCATCATGGACATCTCCGCTTTCCGCGAGCTAGT‑CAATGATTATGCCACCTGTAATCGCGTGCTGGGCCCGCTGCACAGCCTGTACAAGCGGCGACCGGGCCGCGTGAAAGACTATATc                                                                                                                        >  2:178858/13‑251 (MQ=255)                                                                                                                      tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGcc                                                                                                                     >  1:410928/1‑251 (MQ=255)                                                                                                                             gaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCaa                                                                                                              <  1:89945/251‑1 (MQ=255)                                                                                                                                  cttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                                         >  1:95766/2‑251 (MQ=255)                                                                                                                                   ttctGCAATCGAAGGGCGTTTGCAGGAAGCGGCAATCCCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTCTTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATACACAGCTTGAACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCCTTCCAAAAGCg                                                                                                         >  1:251885/1‑250 (MQ=255)                                                                                                                                    tctGAAAGCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCGGCGTCGGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGTGTTTTGACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCGGACACCTGTTATCGCGTGCTGGGCCAGATGTACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGTCTATATCGCCATTCCAAAAGCGa                                                                                                        <  2:97032/250‑1 (MQ=255)                                                                                                                                     ctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                                      >  1:59804/1‑251 (MQ=255)                                                                                                                                     ctGAAATCGAAGCGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGacc                                                                                                      >  2:133068/1‑249 (MQ=255)                                                                                                                                                 ggCGTTTGCAGGAAGCGGGCATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGCACGGCTATCAGTCt                                                                                          >  1:252774/1‑251 (MQ=255)                                                                                                                                                        gCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTAGATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcaca                                                                                   <  1:416434/251‑1 (MQ=255)                                                                                                                                                          aGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACAcc                                                                                 >  1:14033/1‑251 (MQ=255)                                                                                                                                                              aGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCg                                                                              >  2:161302/1‑250 (MQ=255)                                                                                                                                                                            gCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCAc                                                               >  1:366561/1‑251 (MQ=255)                                                                                                                                                                                    aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGcc                                                       >  1:396778/1‑251 (MQ=255)                                                                                                                                                                                            ggggaAGCATTTGTATTCGATTTACTGCAAAATGGTGCCCAAAGAGCCGCGTTTTTACGCGATCATGGACATCTTCGCTGTCCGCGTGATCGG‑CAATTATTTTCACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATAGCGCCATTCCAAAAGCGAACGGCTATTAGTCTTTGCACACCTCGATGATCTGCCCGCACGGTGTGCCGGTTGAgg                                               <  2:262977/248‑1 (MQ=255)                                                                                                                                                                                                           ttCGATTTACTGCAAAATGGTGCTCAAAGAGCATCGTTTTTACTCGATAATGCACATCTACGCTTTCAGCGTGATCGT‑CAATGATTCTGACAACTGATATCGCGTGCTGGGCAAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGCCTTTCTCGCCATTCCAAAAGCGCACGGCTATCAATCTTTGAACCACTCGATGATAGGCCCGCACGGAGAGCCGGTTGAGGTCCAGATCCGTCCCg                                >  2:250023/1‑251 (MQ=255)                                                                                                                                                                                                              gATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                               >  1:346287/1‑249 (MQ=255)                                                                                                                                                                                                                 ttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGat                           >  2:228185/1‑250 (MQ=255)                                                                                                                                                                                                                    cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg                       >  1:16860/1‑251 (MQ=255)                                                                                                                                                                                                                     tGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg                       >  1:410479/1‑250 (MQ=255)                                                                                                                                                                                                                              gggCTCAAAGAGCAGGTTTTTTACGCGATCATGGACTTCAACTCTTTCCGCGT‑TTCCTACAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA‑GATGGc            <  2:182579/249‑1 (MQ=255)                                                                                                                                                                                                                                   cAAAGAGCAGCTTTTTTACCTGATCATGTGCATCTGCGCTTGCCTCGTCATCGT‑CAATTATTCTGACACCTGTTATCGCTTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGGGAAAGACTATAGCGCCATTCCAAAAGCGAACGGCTATCAGTCGTTGCACACCTCGATGATCGTCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA‑GGTGGCGgaga       <  2:171526/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaaGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTAAGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCATCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA‑GATGGCgg          >  1:416187/1‑247 (MQ=255)                                                                                                                                                                                                                                       gagCAGCGTTTTTACTCGATCATGCACATCTACGCTTTCCGCGTGATCGT‑AAATGATTCGGCCACCTGTTATCGCGTGCGGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGTCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA‑GATGGCGGAGATggg   <  1:306453/251‑1 (MQ=255)                                                                                                                                                                                                                                        agcagcGTTTTTACTCGATCATGGCCATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCCGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGG‑CCACGATGGCGGAGATgggg  >  2:217828/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  GTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGT‑CAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCA‑GATGGCGGAGATGGGT  >  NC_000913/3822920‑3823401

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