Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	968,837	G→A	D158N (GAT→AAT)	lpxK →	tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	968,837	0	G	A	97.7%	144.2 / NA	44	D158N (GAT→AAT)	lpxK	tetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/21);  minor base C (1/0);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGG  >  NC_000913/968590‑969074                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              cgGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAAtg                                                                                                                                                                                                                                            <  1:85912/251‑1 (MQ=255)          cgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCCTTCTGGCGCAACACCCTGGTGTGCAGATCATCGTAAC‑CGACGACGGGTTACCGCATTACCGTCTGGCGCGGAATGTGGAAAttt                                                                                                                                                                                                                                   >  2:246542/1‑250 (MQ=255)          cgTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTTTCCCGTTCGTTCTGATGCGGGAAAAGCCCTTCTGGCGCAACACCCTGATGTGCAGCTCATCGTAAC‑CGACGACGGGTTTCAGCATTACCGTCTGGCGCGGAATGTGGGAAttt                                                                                                                                                                                                                                   >  1:283207/1‑250 (MQ=255)                 gTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAt                                                                                                                                                                                                                            <  1:246542/251‑1 (MQ=255)                    gTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGAAGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAA‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTTTTg                                                                                                                                                                                                                          >  1:319937/1‑250 (MQ=255)                                  tggtggTAAGGCTGCATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATTGACCTGTGTTTTTTTTTCAACGCACTGATGCGCCTGGTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAACCATTCTGGCGCCACACCCTGATTGGCAGGTCATCGTAACACG‑CGACGGGTTACAGCCTTACCGTCCGGCGCGTCACGTGGAAATTGTCGGTATTGATGGGGTGGGTCGCt                                                                                                                                                                                                           >  1:343710/1‑251 (MQ=255)                                      ggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACCACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCACCGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTgg                                                                                                                                                                                                       <  1:183905/251‑1 (MQ=255)                                      ggtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTgg                                                                                                                                                                                                       >  1:198806/1‑251 (MQ=255)                                                          gCTGTTATTGTCGGAAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCATCACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGGTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcgg                                                                                                                                                                                   <  1:272133/251‑1 (MQ=255)                                                             gTTATTGTCGGCAGATACCACAACTGCACAGGCGGGTGATGAACCTGTGTTGATTTAGCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGTGTAATGTGGAAATTGTCGGTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcggc                                                                                                                                                                                  <  2:43614/249‑1 (MQ=255)                                                             gTTATTGGCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATGACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCtggcgg                                                                                                                                                                                <  2:343710/251‑1 (MQ=255)                                                                                            gCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCAATGGGTGAGGGGGCGGGGCGCGTaa                                                                                                                                                  >  2:178724/1‑250 (MQ=255)                                                                                                      aaCCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCGAGGCGGGAGCGGGGGGGGTGCTTaacggaggttggg                                                                                                                                       >  2:15727/1‑240 (MQ=255)                                                                                                        ccTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTAGGTTGATGc                                                                                                                                     <  1:375263/251‑1 (MQ=255)                                                                                                          tgtgTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCACTGCGTGAGCGAGCGGGGGGCTTAAAGTCGGTTGATGCgg                                                                                                                                   >  1:244093/1‑251 (MQ=255)                                                                                                                taTTTATAAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTATGATGCGGTAAAAGCCATTATGGCGCCACACCCTGCTGCGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTACTGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCCTGCGTGGGCGCGCGGGGGGCTTACAGTCGGtttgtgtgtggaacc                                                                                                                             >  2:75816/2‑238 (MQ=255)                                                                                                                gATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATc                                                                                                                              <  1:429387/250‑1 (MQ=255)                                                                                                                 aTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGAGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGt                                                                                                                            >  1:275619/1‑251 (MQ=255)                                                                                                                    tATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGCGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCa                                                                                                                          >  1:435201/1‑250 (MQ=255)                                                                                                                     agcAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTGTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCTTTAATCGTCAAc                                                                                                                        <  2:283207/249‑1 (MQ=255)                                                                                                                        aaCGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGGAATCGGCAAcggc                                                                                                                     >  1:428418/1‑251 (MQ=255)                                                                                                                                      ccTGTTGCGGTTTCGCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAg                                                                                                       <  2:413207/251‑1 (MQ=255)                                                                                                                                           tGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGGTGCCGGCGGGGCCAATGCGTGGGCGAG‑GGGGCGCTTAAAGGCGGTTGATGTGGTAATCGTCACCGGTGGGGTCCCCCGccgaggtgg                                                                                                 >  2:284994/1‑243 (MQ=255)                                                                                                                                              ggTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCCATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGGGGCGGGAATCGTCAACGGCGGTGTCCCCCGCAGCggggtg                                                                                                >  2:98647/1‑246 (MQ=255)                                                                                                                                              ggTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGaaa                                                                                               <  1:316740/251‑1 (MQ=255)                                                                                                                                                 ttCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCtt                                                                                                                                                    <  2:43156/195‑1 (MQ=255)                                                                                                                                                 ttCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCtt                                                                                                                                                    >  1:43156/1‑195 (MQ=255)                                                                                                                                                     cccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCAt                                                                                        <  1:106418/251‑1 (MQ=255)                                                                                                                                                           cgttcTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATct                                                                                  <  1:221423/251‑1 (MQ=255)                                                                                                                                                                 gATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGcc                                                                            <  1:182214/251‑1 (MQ=255)                                                                                                                                                                           aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCg                                                                  <  2:347257/251‑1 (MQ=255)                                                                                                                                                                           aaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCCCGCAGCGGTGGAATTCCCATTCATCTGCTGCCGGGTCCGGCg                                                                  >  1:330294/1‑251 (MQ=255)                                                                                                                                                                             gCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTTAAGGCGGGTGATGGGGTAATCGTCCACGGCGGTGTCCCCCGCCGAGGGGAAATCCCCACGCaacttcctcccgggtaggcggg                                                                >  2:108498/1‑230 (MQ=255)                                                                                                                                                                             gCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGt                                                                <  2:135697/251‑1 (MQ=255)                                                                                                                                                                                       cgcAACACCCTGATGTGCAGATCATCGTAAC‑CGATGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAAGTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACg                                                       <  2:315243/250‑1 (MQ=255)                                                                                                                                                                                         cAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGGCAGGCGGGGGAATTACGTa                                                     >  1:49297/1‑250 (MQ=255)                                                                                                                                                                                          aaCACCCTGATGTGCAGATCAGCGTAAC‑CGACGACGGTTTACAGCATTACCGACTGGCGCGTAATGTGGAAATTGGCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCTATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTAcc                                                   <  1:233498/251‑1 (MQ=255)                                                                                                                                                                                               ccTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCCGCGGGGGAATTCCCCTGCATCTGTCGCCGGGTCCGGCGGGGGATTTTCGTACCGGTAc                                              >  1:230046/1‑251 (MQ=255)                                                                                                                                                                                                       gCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCACGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATTCGTGAGCGAGCGGGGCGCTTCACGTAGGTTGATGCGGGCATCGTCACCGGCGGTGGCCCTTGCGGCGGGGGAATCCCCATGCATCCGCTGCCGGGTCAGGCGGGTAATTTTCGTAcccggacgcctagtg                                      >  1:32090/1‑238 (MQ=255)                                                                                                                                                                                                           atcatcGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCTCGCAGCGGGGGAATCCCCCTGCCTCTGCTGCCGGGGCAGGCGGGGGATTTACGGACCCGCACGCGCtgtggagt                                  >  1:397143/1‑247 (MQ=255)                                                                                                                                                                                                                                     aGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCg         <  2:198806/251‑1 (MQ=255)                                                                                                                                                                                                                                         ttACAGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGGTGATGCGGGAATCGTCAACGGCGGGGTCCCCTGCAGCGGGGGAAAACCCATGCATCCGCAGCCGGGGCAGGCGGGGGATTTTCGCACCGGGACGCGTTGTGACGTTGCTccgacttaaaatggtgtggcgcggg     >  2:89626/1‑227 (MQ=255)                                                                                                                                                                                                                                          tACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGGGAAATCCCCCTGCATCTGCTGCCGGGGCCGGCGGGGGATTTTCGTACCGGGTCGCGGTGTGGCGGTGCTCAGCTTGGACCTTTAGTggggagggg    >  1:46205/1‑244 (MQ=255)                                                                                                                                                                                                                                            ccGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCgg  <  1:98647/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              CGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAAC‑CGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGG  >  NC_000913/968590‑969074

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