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	95.3%	130.0 / NA	43	T156T (ACT→ACC)	phoP	DNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (19/22);  minor base T (0/2);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCCGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAAT  >  NC_000913/1189734‑1190225                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      cccGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCGCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGgtgt                                                                                                                                                                                                                                                   <  1:421878/251‑1 (MQ=255)  ccGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTa                                                                                                                                                                                                                                                  <  2:547733/251‑1 (MQ=255)   cGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTAt                                                                                                                                                                                                                                                 >  1:223974/1‑251 (MQ=255)         gagaGCAGGAAACAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGGATCAACGTTTCCATAATGGTGAATTCGAAc                                                                                                                                                                                                                                           >  2:512789/1‑251 (MQ=255)                    aaaaaGACGCAATAATTTTTTCATCAGCGCAATTCGACCAGATAACACTGGCCGCGAACGGTGGTAATCACTTCTTGGGGAAATTGTGCCTGACTTTTTTTGCGCAGACGTCCCATCAGTACATAAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGATGGAGCATTAACGAATCATTGCTCGCCCCTTTGCCCTTATTGAGTATCAACGTTTCCATAATGGTTTATTCGACCGCGGTCAGttt                                                                                                                                                                                                                                >  2:536153/1‑251 (MQ=255)                      aaaGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTg                                                                                                                                                                                                                               >  1:582854/1‑250 (MQ=255)                           cGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACt                                                                                                                                                                                                                         >  1:462490/1‑251 (MQ=255)                            gCAGTAATGTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACtt                                                                                                                                                                                                                        <  1:619236/251‑1 (MQ=255)                                    tttttCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTTGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGCTAGAGTTGGAGCATTAACGCATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTAtt                                                                                                                                                                                                                                                >  2:103060/1‑219 (MQ=255)                                    tttttCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTAtt                                                                                                                                                                                                                                                <  1:103060/219‑1 (MQ=255)                                                                      gacttgAACGGTGTGAATCTCTGCTTGTGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTTGATCAATAGTTGGGCTTTCCCGCAGCGTCGTATCCTGATAGAGGTGGAGCATTAACGAATCTTTTCTGTCCACTTTTCCATTATTGCGTATCAACGGTTTCATAATTGTGTATTCGAACGCGGTGAGTTTGATCTCGTCTTCAGTAATAGATAATTCATGGCGATATAGAACAACCCGaa                                                                                                                                                                              <  2:678535/246‑1 (MQ=255)                                                                         gagaacggcaggatgcACTTCCTGGGGATATTGTGCATGAATTCTTTCTCGCAGCCGTCCCATCAGTACAACAATGGTATGGCGTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAGCGAACCTTTGCTGACCACTTTGCCATTATTGCGTAGCAACGTGTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCTTCACGAATAGAGAAGTCACGGCGAGAGACATCAACTTGTAACg                                                                                                                                                                           <  2:250007/236‑1 (MQ=255)                                                                           gAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCCATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATc                                                                                                                                                                                                                            >  1:120450/1‑200 (MQ=255)                                                                           gAACGGTGGTAATCACTACTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGGCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATc                                                                                                                                                                                                                            <  2:120450/200‑1 (MQ=255)                                                                            aaCGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACgggg                                                                                                                                                                        <  1:660622/251‑1 (MQ=255)                                                                                            tcttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATAAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGCTAGAGTTGGAGCATAAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTACATAATGGTGTATTCGAACGCGTTCAGTTTGATCACTTCGTCATTAATAGATAATTCCCGGCGAGAGAGCTCAACCTGAAACGGGGGTAGCGACATGCCCtg                                                                                                                                                        >  2:64963/1‑251 (MQ=255)                                                                                               tGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCtgtg                                                                                                                                                      >  1:432426/1‑250 (MQ=255)                                                                                                         gtgCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGAcc                                                                                                                                            <  1:240610/250‑1 (MQ=255)                                                                                                         gtgCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAAGTCACGGCGAGAGAGATCATCCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGAcc                                                                                                                                            <  1:204037/250‑1 (MQ=255)                                                                                                          tgCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACggggg                                                                                                                                                                       >  1:469108/1‑222 (MQ=255)                                                                                                          tgCCTGAATTTTTTTGCGCAGACGTCCCAACAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTTCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGGCAGTTTGATCACTTCGGCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACggggg                                                                                                                                                                       <  2:469108/222‑1 (MQ=255)                                                                                                                   ttttttGCGACGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATATAGTTGGAGCATTAACTAATCTTTGCTGACCACTGTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGGCATTAATAGATAATTCACGGCGAGAGAGATCTACCTGAAACGGGGGGAGCGAAATGTCCTGTGAAGCCAGACCTCTATTAcgc                                                                                                                                  <  2:396381/250‑1 (MQ=255)                                                                                                                                              aTCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGCTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATc                                                                                                       >  1:170824/1‑250 (MQ=255)                                                                                                                                                cAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCAc                                                                                                     <  1:347355/250‑1 (MQ=255)                                                                                                                                                  aTGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACctc                                                                                                  <  1:599298/251‑1 (MQ=255)                                                                                                                                                    ggTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAGTGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCtt                                                                                                <  1:42377/251‑1 (MQ=255)                                                                                                                                                       aTGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCaa                                                                                             <  1:252479/251‑1 (MQ=255)                                                                                                                                                               cccGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGAATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGAACGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAAc                                                                                     >  2:600094/1‑251 (MQ=255)                                                                                                                                                                        ccGCATTCGGATATAGTTGGAGCATTAACTAATCTTTGCTGACCTCTTTGCCATTATTTCGTATCAACGTTCCCATAATGGTGTATTCGAACGCGGTCGGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTTTGAAGCCAGACCGCTTTTACGCCGCATTAATGCCTGCATTCGTGCCATCACCTCGTCAATACGGAACGGTTTAGTc                                                                            <  2:652011/251‑1 (MQ=255)                                                                                                                                                                        ccGCATCCGGATAGAGGTGGCGCATTAACGAATCTTTGCTGACCTCTTGGCCATTATTGCGTATCAACGTTGTCATAATTGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGAAAATTCACGGCTAGAGAGATCAACCTGAAACGGGGGGAGCTAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATGCGCGCCATCACCTCTTCAATATGGACCGGTTTAGTc                                                                            <  2:254353/251‑1 (MQ=255)                                                                                                                                                                          gCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCgagaga                                                                                                                                                                                          >  2:115058/1‑139 (MQ=255)                                                                                                                                                                          gCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCgagaga                                                                                                                                                                                          <  1:115058/139‑1 (MQ=255)                                                                                                                                                                                 gATAGAGTTGGAGCATTAACGAATCTTTGCTGATCACTTTGCCATTATTGCGGATCAACGTTTCAATACTGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTCATCGATAATTCACGGCGAGAGAGATAAACCTGAACCGGGGGGAGCGAAATGACCTGTAAAGCCAGACCGCTATTCCCACACATTAATGCCTGAATTCGCGCCATAACCTCTTCAATATGAAACGGTTTAggaccataatc                                                                    >  2:657419/1‑240 (MQ=255)                                                                                                                                                                                     gagTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAatcatc                                                                 >  1:347183/1‑249 (MQ=255)                                                                                                                                                                                              cATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCAATTCGTCATTAATAGATAATTCACGGCGCGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTAAGACGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGACAAGGTTTAGTCACATAATCATCAGCACCGTCAc                                                      >  2:248743/1‑251 (MQ=255)                                                                                                                                                                                                            tGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACt                                        >  1:604400/1‑251 (MQ=255)                                                                                                                                                                                                             gCTGACCCCTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCGCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACt                                        <  1:334192/250‑1 (MQ=255)                                                                                                                                                                                                             gCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATAACTTCGTCATTAATAGATAATTCACGGCGAGCGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGAATTAATGCCTGCATTCGCGCCACCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACtt                                       >  2:625523/1‑251 (MQ=255)                                                                                                                                                                                                             gCTGACCACGTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGGATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACtt                                       <  2:347183/251‑1 (MQ=255)                                                                                                                                                                                                              cggACCACTTTGCCATTATTGCGTATCAACGTTTTCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGCCAATAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGAACGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAAGATGAAACGTTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACtta                                      <  2:49618/249‑2 (MQ=255)                                                                                                                                                                                                                              tattGCGAATCAACATTACCATAAGGGTGTATTCGACCGCGGTATGTTTGATCACTTCGTCATTAATAGATAATTACCGGCGAGAGCGTTAAACCTGAAACGGGGGAAGCGAAATGCCCTGTGCAGCCAGACCGCTATTAAGCCGAATTAATGCCTGCATTCGCGCCATCACCTCTTCATCATGAAAAGGTTTAGTCAAATACACATCCGCACCGGCAATTAAAAACACGACTTTGTCCTGCAAACTTTc                       >  2:364193/1‑250 (MQ=255)                                                                                                                                                                                                                                    gTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCg                <  1:521600/251‑1 (MQ=255)                                                                                                                                                                                                                                          aCGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCTCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAAt          <  1:600094/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ataatgGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAAt  >  1:180251/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      CCCGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAAT  >  NC_000913/1189734‑1190225

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