Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,341,542	G→A	A329T (GCG→ACG)	lapB →	lipopolysaccharide assembly protein B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,341,542	0	G	A	95.5%	125.7 / ‑4.2	46	A329T (GCG→ACG)	lapB	lipopolysaccharide assembly protein B
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (23/19);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTGGTCAGCGAAACGCTGGAAATGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATG  >  NC_000913/1341296‑1341774                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         cTGGTCAGCGAAACGCTGGAAATGTTGCAAACCTGCTATAAGCAGTTGGGTAAAACTGACGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGGGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACACTGCGTGAGTTCCATAAGTAAATGGATACCCACTTAAATGAAACgg                                                                                                                                                                                                                                       >  2:608543/1‑250 (MQ=255)                      aTGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCCGCGTAATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAagagag                                                                                                                                                                                                                 <  1:165860/251‑1 (MQ=255)                       tGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAACCGGAAGAAGGGCGGGCCAAAGGGAGc                                                                                                                                                                                                                >  2:8655/1‑251 (MQ=255)                           gCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTg                                                                                                                                                                                                             <  1:15229/250‑1 (MQ=255)                                  tGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGtgct                                                                                                                                                                                                     >  1:572018/1‑251 (MQ=255)                                          gcagTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTTATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACa                                                                                                                                                                                             <  2:344172/251‑1 (MQ=255)                                                  ggaaaaCTGCCGAATGTGCAGAATTCCTGCTGCGTGCGGTGGAAGAGAACCCCGGGGCCGATGGTGAATTGATGCTTGCTGATATGATCGAAGCGCGCGACGGTAGTGAGGCCGCACATGTATATATTACTCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGCGCCTGATGGTGCTGCGTGACATTGTTgga                                                                                                                                                                                     <  2:214144/249‑2 (MQ=255)                                                     aaaCTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGGTGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGCCCATGCGGGTGTTCCATAAGTTAATGGATTACCACTTAATTGAATCGGAAGAAGGGCGTGCCAAAAAGAGCCTGATGGTGCTGCGTGACATGGTTggag                                                                                                                                                                                    >  2:503021/1‑247 (MQ=255)                                                              aaTGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACgt                                                                                                                                                                         >  2:157667/1‑251 (MQ=255)                                                                tGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACgt                                                                                                                                                                         >  1:318920/1‑249 (MQ=255)                                                                              gCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTAt                                                                                                                                                          <  1:180554/250‑1 (MQ=255)                                                                                aGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAACGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGt                                                                                                                                                       >  2:212495/1‑251 (MQ=255)                                                                                 gCGTGCGGGGGAAGAGAACACCGTTGCCGAGGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGGCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGtt                                                                                                                                                      <  2:239804/251‑1 (MQ=255)                                                                                     gCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGcc                                                                                                                                                   <  1:258876/250‑1 (MQ=255)                                                                                             agagAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCATAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGAAAGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgtg                                                                                                                                          >  2:19447/1‑251 (MQ=255)                                                                                             agagAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGCGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTCCGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgtg                                                                                                                                          >  2:665969/1‑251 (MQ=255)                                                                                                gaACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtt                                                                                                                                       <  2:318920/251‑1 (MQ=255)                                                                                                 aaCACCGGTGCCGATGCTGAATTTATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTGATCGTTGCCAGAAATGTGGttt                                                                                                                                      <  2:572018/251‑1 (MQ=255)                                                                                                  acacCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtttt                                                                                                                                     <  1:212282/251‑1 (MQ=255)                                                                                                     ccGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTAcc                                                                                                                                  <  1:665969/251‑1 (MQ=255)                                                                                                             gATGCTGAATTGATGCTTGCTGATATCCTCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATacac                                                                                                                          >  1:273134/1‑251 (MQ=255)                                                                                                              aTGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCAATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGACTGATGGTGCTGAGTGACATGGTTGGCGAGACGGTAAGTAGTAAGCAACGTTATCGTTGCCAGCACTGTGGTTTTACCGACTACAcc                                                                                                                         >  2:614699/1‑251 (MQ=255)                                                                                                                    aaTTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTGCGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCATGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTTTCGTTACCAGAAATGTGGTTTTACCGCATACACCCTCTAc                                                                                                                   >  1:530579/1‑251 (MQ=255)                                                                                                                                 tGATATCATCGAAGCGCGCGCCGGTAGTGAAGCCGCACAGGACTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGCAGGGCGTGCCAAAGAGTGCCTGATGGTGCTGCGAGACATGGTTGGCGAAAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAtttacagt                                                                                                      >  2:310575/1‑245 (MQ=255)                                                                                                                                    tatCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCTACCATGCGTGTGTTCCATAAGTTAATGGCTTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGATGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCAAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCtt                                                                                                   >  2:626547/1‑251 (MQ=255)                                                                                                                                                  gcgACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTg                                                                                                                                                                                       >  1:176075/1‑153 (MQ=255)                                                                                                                                                  gcgACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTg                                                                                                                                                                                       <  2:176075/153‑1 (MQ=255)                                                                                                                                                   cgACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCa                                                                                     >  1:591037/1‑250 (MQ=255)                                                                                                                                                    gACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAAc                                                                                   >  1:569124/1‑251 (MQ=255)                                                                                                                                                     aCGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAAcc                                                                                  >  2:572366/1‑251 (MQ=255)                                                                                                                                                                                    tcgcCAGCTTCAGCGTCATCTGACCATGCGTGTGTTCCATAAGTTACTGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGCGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTATACCGATTCGCGGTCTTGATGGACTGt                                                   >  2:179163/2‑251 (MQ=255)                                                                                                                                                                                      tccAGCTTCAGCGTCATCCTACCATGCGAGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTAGGCGAGAAGGTACGAAGTAAGCCTCGTTCTAGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGACATTGTCCGTCTTGTCGTGCCTGGTCAGCCATTAAACCGATTAGCGGTATTGATTGCCTGTa                                                  >  2:548531/2‑250 (MQ=255)                                                                                                                                                                                      gcCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTaa                                                 <  1:572366/251‑1 (MQ=255)                                                                                                                                                                                        cagcgtcagcgtCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTGAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAAtt                                               <  1:617504/251‑1 (MQ=255)                                                                                                                                                                                          gCTTCATCGTCATCCGTCCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATTGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTTTCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTGGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAAtttt                                             <  2:645312/251‑1 (MQ=255)                                                                                                                                                                                          gCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAAtttt                                             >  1:520875/1‑251 (MQ=255)                                                                                                                                                                                            ttCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAAttttt                                            >  1:236162/1‑250 (MQ=255)                                                                                                                                                                                                  gTCAACCGACCATGCGTTTGTTCCAAAAGTTAATGTATGACCTCTTAAAGGAAGCGGAAGAAGGGCGGGCCAAAGAGAGCCTGATGGTGCTGCGTGTCATGGTTTGCGAGAAGGTACGGAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTGTACCTCATACACCCTCGACTGGCATTTTCCGTCTTGTCGGGCCTTGTCAACCATAAAACCGAGTCGCTGTCGTGATGGCCTGTAATTTTTaaaaaaa                                     <  2:634969/251‑1 (MQ=255)                                                                                                                                                                                                        cGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAt                                                                                                             >  1:137815/1‑173 (MQ=255)                                                                                                                                                                                                        cGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAt                                                                                                             <  2:137815/173‑1 (MQ=255)                                                                                                                                                                                                        cGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGGCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGa                               <  2:445087/251‑1 (MQ=255)                                                                                                                                                                                                             aTGCGTGTGTGCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTa                          <  1:657385/251‑1 (MQ=255)                                                                                                                                                                                                                   gtgtTCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAACGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTac                     >  1:634580/1‑250 (MQ=255)                                                                                                                                                                                                                      ttCCATAAGTTAATGGATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTacaac                  >  1:310414/1‑250 (MQ=255)                                                                                                                                                                                                                                    ggATTACCACTTAAATGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTaa    <  1:212495/250‑1 (MQ=255)                                                                                                                                                                                                                                     gATTACCACTTAAAGGAAACGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATg  <  1:520196/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         CTGGTCAGCGAAACGCTGGAAATGTTGCAAACCTGCTATCAGCAGTTGGGTAAAACTGCCGAATGGGCAGAATTCCTGCAGCGTGCGGTGGAAGAGAACACCGGTGCCGATGCTGAATTGATGCTTGCTGATATCATCGAAGCGCGCGACGGTAGTGAGGCCGCACAGGTCTATATTACGCGCCAGCTTCAGCGTCATCCGACCATGCGTGTGTTCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGTCAACCATTAAACCGATTCGCGGTCTTGATGGCCTGTAATTTTTAAAAAAAATCCGACTTTAGTTACAACATACTAATTATTAATG  >  NC_000913/1341296‑1341774

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