Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,175,517	A→G	T31A (ACC→GCC)	lolC →	lipoprotein release complex ‑ inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,175,517	0	A	G	95.2%	137.1 / ‑1.7	43	T31A (ACC→GCC)	lolC	lipoprotein release complex ‑ inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (22/18);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

CGGCAAGATTGATGATAAAACATGGCAATTTAGCCGATTGATTTACGGGGGCTTTTCAGATTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTACCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTG  >  NC_000913/1175275‑1175753                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             cggcAAGATTGATGATAAAACATGGCAATTTAGCCGATTGATTTACGGGGGCTTTTCAGATTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTACCATCGGc                                                                                                                                                                                                                                      >  1:489793/1‑251 (MQ=255)           gatgatAAAACATGGCAATTTAGCCGATTGATTTACGGGGGCTTTTCAGATTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCg                                                                                                                                                                                                                            >  1:190464/1‑251 (MQ=255)                 aaaaCATGGCAATTTAGCCGATTGATTTACGGGGGCTTTTCAGATTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTg                                                                                                                                                                                                                       <  1:295385/250‑1 (MQ=255)                         gCAATATAGCCGATTGATTTACGAGGGCTTTTCAGATTAGCCCTGACGATCACTTAAAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTTTCTTTGCTACAGCAACCAGCCGGATTTCATGTACCACCCGGTCTCTCTATTTATTGGCATG‑CGTTACATGCGTGTGCGAGCAGCGGTTCGCTTCGGTCGTTTCTTCTCCTGGATTTCTGCCATCGGCATTACCCTCGGGGGGATGTCGCTg                                                                                                                                                                                                              >  2:418206/1‑251 (MQ=255)                               tAGCCGATTGATTTACGGTGGCTTTGCAGATTAGCCCTGACGATCACTTACAGTTCAGTCTTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCGTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTTGTcac                                                                                                                                                                                                         <  2:248001/250‑1 (MQ=255)                                        gattTACGGGGGCTTTTCAGCTTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGTCATTACCCTCGGGGTGATTGCGCTGGTCACCGTTTTGTc                                                                                                                                                                                                >  2:319637/1‑250 (MQ=255)                                              acccgacgctcgttcgatccccTGACTATCACTTTCAGGTCATACGTTTTCCCATCTTGCGTTCGCTTATATACTCGTGTCTTGGCTACAGCAACCAGACGGATGTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGTGTGATGGCGCTGGTCACAGTATTGTCAGtgatg                                                                                                                                                                                         <  2:297111/232‑1 (MQ=255)                                                  ggCTTTTCAGATTAGCCGTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCACTTATATACTCGAGACTTTGCTACAGCAACCAGACGGATTTCATATAACAACATGTCGCACTATTTGTTGGCATG‑CGTTACATGCGTGGGCGTGCAGCGGATCGATTCGATCGTTTCGTCTCCTAGCTTTCTGCGTTCGGCATTACCCTCGGGGTTATGGCTCTTGTACCAGTATTGTCAATGATGCACg                                                                                                                                                                                     >  2:204032/1‑251 (MQ=255)                                                  ggCTTTTCAGATTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCAATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGAAACAAGACGGATTTCATGTACCACCCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGtataagtattggcggtggtgaacg                                                                                                                                                                                     >  2:508112/1‑228 (MQ=255)                                                            agtAGCCCTGACGAGCGCTTACAGCTTAGACGGTTACCCTTCTTGCTTTCGCGTATATCCTCGTGGCTTTGATACAGCAACCAGACGGATTTCATGTGCCAACCTGGTGCTCTATTTATTGGACTG‑CGTTACATTCGTGGGCGTGCAGCGTAGCGCTTCGGTCGTTTCGGCTCCTGGCTGTCTGCCATCGGCATTACCCGAGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAgcg                                                                                                                                                                           <  2:4282/249‑1 (MQ=255)                                                              tAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGcgcg                                                                                                                                                                         >  1:254823/1‑251 (MQ=255)                                                                  ccTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGtgatg                                                                                                                                                                                         >  1:297111/1‑231 (MQ=255)                                                                          tCACTAACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTCTATCCTCGTGTCTTTGCTACAGAAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGCTGGCGCTGGTCACAGTATTGTCAGTGATGGACGGCTTTTAGCGCGAGCTGCAAaacc                                                                                                                                                             >  2:549587/1‑250 (MQ=255)                                                                               tACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGGACCAACCTGTCGCTCTATTTAGTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACAt                                                                                                                                                         <  1:359314/250‑1 (MQ=255)                                                                                               aCCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCACCCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGGGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGGGCTGCAAAACCAAATCCTTGGCCTGATGcccc                                                                                                                                        >  2:320058/1‑249 (MQ=255)                                                                                                   aTCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTCTTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAACACAACATCCTTGGCCTGATGCCACAGGc                                                                                                                                    >  1:361118/1‑251 (MQ=255)                                                                                                    ccTTGCTTTCGCTTATATACGCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGGTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTACCAGCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCa                                                                                                                                   <  2:489793/250‑1 (MQ=255)                                                                                                        gCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAAtt                                                                                                                                <  1:320058/250‑1 (MQ=255)                                                                                                             cGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTctctct                                                                                                                          <  1:191970/251‑1 (MQ=255)                                                                                                               cTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTTCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCtcttc                                                                                                                        >  1:238330/1‑251 (MQ=255)                                                                                                                   gataCTCGTTTCTTTGCTACATCAACCAGACGTATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTTGCCTGATGCCACAGGCAATTCTCTCTTCTGAg                                                                                                                    <  1:151948/250‑1 (MQ=255)                                                                                                                      aCTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACCGGCAATTCTCTCTTCTGAGCAt                                                                                                                 >  2:53814/1‑251 (MQ=255)                                                                                                                        tCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATgg                                                                                                               >  1:544710/1‑251 (MQ=255)                                                                                                                            gtCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCtt                                                                                                                                                                                 >  1:21120/1‑181 (MQ=255)                                                                                                                            gtCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCtt                                                                                                                                                                                 <  2:21120/181‑1 (MQ=255)                                                                                                                                  gCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGACTGGCTCCCTTAAc                                                                                                     >  2:492102/1‑251 (MQ=255)                                                                                                                                         cAACCAGACGGATTTCATGTACCAACCTTTCGCTCTATTTATTGGCCTC‑CGTTTCTTGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTTCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGATCATGGCTCTCTTATCCCgcagc                                                                                              <  2:418047/251‑1 (MQ=255)                                                                                                                                                gttgttttCATGTACCAACCTAGCGCTCTGTTTGTTGTCCCG‑CGTATCAGTCGTTGGCGTGCAGCGTATCTCTTCGGTCGATTCGTCTCCTGGCTTTCTGCCGTCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCGTTGAGCGCGAGCTGCGAAACAGCATCCGTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTccc                                                                                       <  2:105853/246‑1 (MQ=255)                                                                                                                                                 cGGATTTCATGTTCCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCa                                                                                      <  2:238168/251‑1 (MQ=255)                                                                                                                                                    aTTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGaa                                                                                   <  1:184322/251‑1 (MQ=255)                                                                                                                                                           gTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCa                                                                             <  1:319637/250‑1 (MQ=255)                                                                                                                                                              ccAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCCTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCAGATCATGGCTCTTTTAACCCGCCACAACTCCCAAATACGACAGTcc                                                                         >  2:312102/1‑250 (MQ=255)                                                                                                                                                               cAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGACCGGCTTTTAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTATTAACCAGCAGCAACTCCCAGAAACAGCAGtaa                                                                         >  2:205044/1‑248 (MQ=255)                                                                                                                                                                aaCCTGTCTCTCTATTTATTGGC‑TGTCGTTACATGCCTGGTCGTGCAGCGGATCGCGTCGGTCTTCTCGTCGCCTGGCTTTCTGCCAGCGGCATTGCCCTCGGGGTGATTGCGCTGGTCACAGGATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAATACAACATCCTTGGCCTGATCCCACAAGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGAAGCAACTCCCAGATACGTCAGTCaaa                                                                       <  2:363470/251‑1 (MQ=255)                                                                                                                                                                       cGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGCTGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCAAGAAACGGCAGTCAAACTGGACg                                                                >  2:136603/1‑251 (MQ=255)                                                                                                                                                                             atttattGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGCGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGtt                                                           >  2:112490/1‑250 (MQ=255)                                                                                                                                                                                           cGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCAcc                                             <  1:6517/251‑1 (MQ=255)                                                                                                                                                                                            gTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCt                                            >  1:321155/1‑251 (MQ=255)                                                                                                                                                                                                    cGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATtactac                                    >  1:181516/1‑251 (MQ=255)                                                                                                                                                                                                    cGTGGGCGTGCAGAGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTATTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCAAGAAACGGCAGTCCAACTGGACGGCGTTAATCGAGTCGCACCTATtactac                                    >  1:358643/1‑251 (MQ=255)                                                                                                                                                                                                     gTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGGCTCCGGGCTTTCTGCCATCGGCATTACCCTCGGGGGGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCAGGGCTCTCTTAGCCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCGCCGATTACTAcc                                   <  2:254823/251‑1 (MQ=255)                                                                                                                                                                                                                               cGTTTCGTCTCCTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGTACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAgcgcgc          <  1:549587/250‑1 (MQ=255)                                                                                                                                                                                                                                         ccTGGCTTTCTGCCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTg  <  1:53814/248‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             CGGCAAGATTGATGATAAAACATGGCAATTTAGCCGATTGATTTACGGGGGCTTTTCAGATTAGCCCTGACGATCACTTACAGTTCAGACGTTTACCCATCTTGCTTTCGCTTATATACTCGTGTCTTTGCTACAGCAACCAGACGGATTTCATGTACCAACCTGTCGCTCTATTTATTGGCCTG‑CGTTACATGCGTGGGCGTGCAGCGGATCGCTTCGGTCGTTTCGTCTCCTGGCTTTCTACCATCGGCATTACCCTCGGGGTGATGGCGCTGGTCACAGTATTGTCAGTGATGAACGGCTTTGAGCGCGAGCTGCAAAACAACATCCTTGGCCTGATGCCACAGGCAATTCTCTCTTCTGAGCATGGCTCTCTTAACCCGCAGCAACTCCCAGAAACGGCAGTCAAACTGGACGGCGTTAATCGCGTCGCACCTATTACTACCGGTGATGTGGTACTGCAAAGCGCGCGCAGCGTG  >  NC_000913/1175275‑1175753

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