Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,189,980	A→G	T156T (ACT→ACC)	phoP ←	DNA‑binding transcriptional dual regulator PhoP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,189,980	0	A	G	97.7%	140.5 / NA	44	T156T (ACT→ACC)	phoP	DNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/14);  minor base T (0/1);  total (29/15)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACC  >  NC_000913/1189738‑1190220                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  cAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTc                                                                                                                                                                                                                                          >  1:510702/1‑251 (MQ=255) cAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTc                                                                                                                                                                                                                                          <  1:576904/251‑1 (MQ=255)  aGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCACTTCCAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATCTTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGATTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCg                                                                                                                                                                                                                                         >  1:400750/1‑251 (MQ=255)   gCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGCCCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCg                                                                                                                                                                                                                                         >  2:207522/1‑250 (MQ=255)      agagCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAAcg                                                                                                                                                                                                                                     <  1:29465/251‑1 (MQ=255)           gggAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCGGGCCGCGAACGGTGGTAATCATTTCGTGGTGATATTGTGCCTTAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTAGGGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTc                                                                                                                                                                                                                                <  2:642166/251‑1 (MQ=255)             gAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAg                                                                                                                                                                                                                              <  1:426096/251‑1 (MQ=255)               aaaaaaGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCAGCGAACGGTGGTAATCACTTCTTGGGGCTATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGGG‑TATTCGAACGCGGTc                                                                                                                                                                                                                                >  1:444521/1‑247 (MQ=255)                  aaaGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTGGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGa                                                                                                                                                                                                                         <  2:137470/251‑1 (MQ=255)                     gACGCAGTAATTTTTTCATCAGCGCACTTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATc                                                                                                                                                                                                                       >  1:623059/1‑250 (MQ=255)                      aCGCAGTAATTTTTTCATCAGCGCAATTCGAACAGCTAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCGTCAT‑GTGGTATTCGAACGCGGCCAGTTTgtttcc                                                                                                                                                                                                                     >  2:290856/1‑246 (MQ=255)                              aTTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCAt                                                                                                                                                                                                             <  1:657042/251‑1 (MQ=255)                                       tcaGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTTGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGACTCCGGTAAGAGTTGGAGCATTAACGAATCTTTGCTGCCCATTTTGCCATTATAGCGTATCAACGTTTCCATCATGGTG‑TATTCGAACGCGGTCAGTTTTATCACGTCTTCTTTAATTGAt                                                                                                                                                                                                     >  2:56873/1‑250 (MQ=255)                                           cgcAAGTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTc                                                                                                                                                                                                <  2:248038/251‑1 (MQ=255)                                             cAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGAAATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCACTGGTATGTCTTTCCCGCAGCTCCGCATCCGGATAGTGTTGGAGCATTAACGCATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCCGTTTGATCACTTCGTCATTAATAGATAATTctc                                                                                                                                                                                              >  2:162455/1‑249 (MQ=255)                                                ttCGAACAGATAGCCCTGGCCGCGAACGGTGGAAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGc                                                                                                                                                                                           >  1:244427/1‑251 (MQ=255)                                                    aaCAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATACTGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATACTTCACGGCg                                                                                                                                                                                          >  1:481785/1‑248 (MQ=255)                                                                        aaCGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACgggg                                                                                                                                                                   >  1:191504/1‑251 (MQ=255)                                                                                  aTCACTTCTTGGGGATATTGAGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCACTGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAGCGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATAAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGCAAt                                                                                                                                                         >  1:340392/1‑251 (MQ=255)                                                                                         cttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGGGCGAAATGACAtgt                                                                                                                                                  >  2:519386/1‑251 (MQ=255)                                                                                                 tactgtgCCTGAATTTTTTTGCGCAGACGTCCCATCAGTATATTAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGACATTATTGAGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTAGTCATTAATAGATAATTCACGGCGAGAGTGCTCATACTGAAACGGGGGGAGCGAAATGACCtgtgcagaca                                                                                                                                           >  2:185931/4‑244 (MQ=255)                                                                                                                             gTCCCATCAGTAGATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACAACTTTGACATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGAGAAATGACCTGTGCAGCAAGACCGATATTAAGACGAATTAATGCCTg                                                                                                               >  2:314252/1‑250 (MQ=255)                                                                                                                                ccATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCAtt                                                                                                           >  1:447526/1‑251 (MQ=255)                                                                                                                                      gTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGACATGACCTGTGACGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGcc                                                                                                     >  2:312552/1‑251 (MQ=255)                                                                                                                                                gttATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTTTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCtt                                                                                           <  2:118011/249‑1 (MQ=255)                                                                                                                                                    tGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCAGCATTAATGCCTGCATTCGCGCCGTCACCTCTTCAat                                                                                       >  1:410595/1‑251 (MQ=255)                                                                                                                                                               cAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGGGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCATCCTGAAACGGGGGGAGCGAAATGACCTTTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGtt                                                                            <  1:302156/251‑1 (MQ=255)                                                                                                                                                                aGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTAc                                                                                                                               <  1:381317/199‑1 (MQ=255)                                                                                                                                                                aGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTAc                                                                                                                               >  2:381317/1‑199 (MQ=255)                                                                                                                                                                    ccGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTc                                                                       >  1:575791/1‑251 (MQ=255)                                                                                                                                                                       cATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTcac                                                                     >  2:137636/1‑250 (MQ=255)                                                                                                                                                                                 gagTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAg                                                          <  1:518954/251‑1 (MQ=255)                                                                                                                                                                                     tGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCAcc                                                      >  1:202778/1‑251 (MQ=255)                                                                                                                                                                                      ggAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGCCACATAATAATCAGCACCg                                                     >  1:565328/1‑251 (MQ=255)                                                                                                                                                                                       gAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCCTCACCTCTTCAATATGAAACGGTTTTGTCACATact                                                                >  2:118055/1‑237 (MQ=255)                                                                                                                                                                                       gAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAat                                                                <  1:118055/239‑1 (MQ=255)                                                                                                                                                                                                 gAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGACACAGGGGGAGAGAAAAGACCTGTGAAGACAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATAACCTCTTCAATATGAAACGGTTTAGTCACATAAACATCAGCACCAGCACTTAATa                                           >  2:512306/1‑250 (MQ=255)                                                                                                                                                                                                          cTGACCAATTTGCCATTATTGCGTGTCACCGTTTCCATAATGGTG‑TATTCGACCGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGGTCAACCTGAAACGAGGGAAGCGAAATGCCCTGTGAAGCCAGACCGCTATTACGCCGCAGTAATGCCTGCATTCGCGCCATCACCTCTTCAATAAGAAACGATTTCGTCACATCATCATCTGCACCGGCACTTAATCCTTCCACttt                                 >  2:621459/1‑251 (MQ=255)                                                                                                                                                                                                            gACCTCTTTGCCATTATTGCGTATCAACGTTTCCATAATGTTG‑TAGTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATGCACGGCGAGAGAGATCAACCTGAAACTGGGGGAGCGAAATGACCTGTGAAGCCAGATCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCGCATAATCATCAGCACCGGCACTTAAGACTTCGACTTTg                                <  2:673020/250‑1 (MQ=255)                                                                                                                                                                                                              ccACTTTGACATTAATGCGAATCAAAGTTTCCATAATGGGT‑TATTCGAACGCGGTAAGTTTGCTAACTTCGTAATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGCGCGAAATGACCTTTGAAGACAGACCGCTATGACGCCGCATTAAAGCCTTCATTCGCGCCATCCCCTCTTCAATATGTAACGGTTTAGTCACAAAATCATCACCACCGGCACTTATTAATTAGACTATGTcc                             >  2:442883/1‑251 (MQ=255)                                                                                                                                                                                                                         ttattGCGTATCAACGTTTCCATAATTGTG‑TATTCGAACGCGGTCAGTTTGAACACTGCGTCATTAATAGATAATTCACGGCGAGATAGATCAACCTGAAACGGGGGGAGCGAAATGTCCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCttg                   <  2:615491/250‑2 (MQ=255)                                                                                                                                                                                                                           attGCGTATCAACGTTTCCATAAGGGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCa                 <  2:191504/250‑1 (MQ=255)                                                                                                                                                                                                                                     aaCGTTTCCATAATGGTA‑TATTCGAACGCGGGCAGTTTGATCAATTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAGCATGAAACGGCGGGAGCGAAATGACCTGTGAAGCCATACCGCTATTACGCCGCATTAATGCCTGCATTCGAGCCCTCACCTCTTCACTATGAAACGGTTTAGTCACATAATCGTCAGCACCGGCCATTAATACTTCGACTTTGTCCTGCCAGATTTCACGGGCGGTTa       >  2:613838/1‑250 (MQ=255)                                                                                                                                                                                                                                           tCCATAATGGTG‑TATTCGAACGCGGTAAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATAcc  >  1:493506/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTG‑TATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACC  >  NC_000913/1189738‑1190220

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