Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,035	C→T	G188D (GGC→GAC)	yjcS ←	linear primary‑alkylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,035	0	C	T	94.9%	106.3 / ‑2.1	39	G188D (GGC→GAC)	yjcS	linear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (19/18);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

ACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACT  >  NC_000913/4305792‑4306272                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              aCGATATTTTTCGTCGGTGCAGTAAGGCCCGGGCCGCCCGTTGCCAGCGTCCCGCAAAGGCCATTCCCCACATTGCCTTGCGCGTTGGGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATa                                                                                                                                                                                                                                        <  1:348073/251‑1 (MQ=255)  cGATAGTTTTCGTCGGGCCAATAATGCTCGGGCCCCCCGTTGCCAGCGTCCCGCCAAGGCCATTCCCCACATTGCCTTGCGCGTTGTGCGGCAGAAACAGACCTTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAg                                                                                                                                                                                                                                       <  1:39255/251‑1 (MQ=255)         tttCGTCGGGGCAATAATGCTCGGGTCGCCCGTTGCCAGCGCCACGCAAAGGCCATTACCTACATTGCCTTGCGCGTGGTGCGGCATTAACAGACCGTAAGAGTAGAGCGCACGGGGGCTAATGATTTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCg                                                                                                                                                                                                                                 <  1:63137/250‑1 (MQ=255)             gtcgGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGtgtg                                                                                                                                                                                                                            >  1:274084/1‑251 (MQ=255)              gcggggCAATAATGCCCGGGTCGCCCGTTGCCAGCGCCACGCCAAGGCCATTACCCCCATGGCCTTGCGCGTTGGGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCTCCATAGTGGTCGGTGTgg                                                                                                                                                                                                                           <  1:108157/246‑1 (MQ=255)                  tGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtgt                                                                                                                                                                                                                       >  1:104076/1‑251 (MQ=255)                   gCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCtgtg                                                                                                                                                                                                                      >  1:180634/1‑251 (MQ=255)                                 gTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCCTGATGTTACCCGCCAGCACGGTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCCGAGATAATGCCTTTCACGTCCCCATAGTGGTCGGTGTGGCTGGGGGTGTAGATAACgg                                                                                                                                                                                                        >  2:380362/1‑251 (MQ=255)                                    gCCCGTTGACAGCGTCACGCCAAGGCCATTACCCACATTACCTTGCGGGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGATACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCCAGAGCACTCACCTGAACTTTGCGGGGTTTAACAT‑GGGCTTCAGCGATAATGACTTTCACTTCCCAATAGTGGACGGTGAGGCtggaggtgtaggtcaacgca                                                                                                                                                                                                      >  2:291957/1‑232 (MQ=255)                                                      gCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGAc                                                                                                                                                                                   >  1:271945/1‑251 (MQ=255)                                                       ccAAGGCCATTACCCCCATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACg                                                                                                                                                                                  <  1:193334/251‑1 (MQ=255)                                                                            ccTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCGTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCa                                                                                                                                                               >  1:253128/1‑249 (MQ=255)                                                                                                  aCAGACCGTAAGAGTAGAGCGCACGGCGGATCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAAGGCAAAAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAggcggg                                                                                                                                       >  1:281218/1‑250 (MQ=255)                                                                                                     gACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTc                                                                                                                                    >  1:2550/1‑251 (MQ=255)                                                                                                                      gcACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCACGGGTCGATAAc                                                                                                                   >  1:64746/1‑251 (MQ=255)                                                                                                                         gggcggcTCATGAGGTCCTCCGCCAGCACGTTTTCGCGGATGGATTCGTCCTTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGAGAACGGCAACAATCGGGTTTTGCGGACGATGCTGGAAGTAAAGGTCAATTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaat                                                                                                                <  2:354656/250‑1 (MQ=255)                                                                                                                               cTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGt                                                                                                                                                                                             <  1:20305/168‑1 (MQ=255)                                                                                                                               cTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGt                                                                                                                                                                                             >  2:20305/1‑168 (MQ=255)                                                                                                                               cTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGGTAACAATAATGcc                                                                                                          >  1:161350/1‑251 (MQ=255)                                                                                                                                       gTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGc                                                                                                  >  1:263268/1‑251 (MQ=255)                                                                                                                                         tccccGCCAGCACATTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTGTGTGGCAGGCGGCATCACCAGCGGGTCGATAACAATAATGCCTTTCTCGcc                                                                                                 <  1:323115/248‑1 (MQ=255)                                                                                                                                            ccGCCAGTACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCa                                                                                             <  1:11357/251‑1 (MQ=255)                                                                                                                                              gCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCTGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCTGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaac                                                                                           <  1:18487/251‑1 (MQ=255)                                                                                                                                                   cACGTTTTCGCTGATTGCTTCGTCCATAAAGCCTGTAGGGACAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTGTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaacg                                                                                      <  2:104076/251‑1 (MQ=255)                                                                                                                                                   cACGGTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaacg                                                                                      <  1:380362/251‑1 (MQ=255)                                                                                                                                                    aCGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCCACCGGTCCCTAACAATAATGCCTTTCTCGCCCTCaacgaacgg                                                                                     >  1:380077/1‑250 (MQ=255)                                                                                                                                                                  ggCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCtt                                                                                                        <  1:9298/218‑1 (MQ=255)                                                                                                                                                                      tCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCTCCATAGTGGTCGGTGTGGCTGTGATTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCTGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATcc                                                                     <  1:377286/249‑1 (MQ=255)                                                                                                                                                                            aTAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTTTCGCCCTCCACGAACGTAATGTTTGAAATATCCTGGCcgcg                                                             >  1:81042/1‑251 (MQ=255)                                                                                                                                                                                      cAGGGGTAAACACCTGAACGTTGCCGGACTTAA‑ATGCTGCTTTAGAGATAAGTCCTTGCACGTCACCATTGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCcgcgaaacttatac                                                   <  2:168125/251‑11 (MQ=255)                                                                                                                                                                                      cAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCCACGAACGTACAGTTAGAGATATCCTGGCCGCGCACCTGATAc                                                   >  1:213939/1‑251 (MQ=255)                                                                                                                                                                                            cAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATACTGCCTTTAACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGGTAAAAATAATGCCTTTCTCGCCCTCAACGAACGTAATTTTAGAGGTATCCCGGCCGCGCAACTGCTACAttttt                                             >  1:223460/1‑250 (MQ=255)                                                                                                                                                                                              aTCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGTTAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGACGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCCCCCTCCACGAACGTAATGTTAGAGGTATCCTGGCCGCGCAACTgcttaattttatt                                           >  1:345006/1‑239 (MQ=255)                                                                                                                                                                                                                gATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCCGTGACTTTGAACAGGcc                         >  1:78702/1‑251 (MQ=255)                                                                                                                                                                                                                   ttAACAT‑CGGCTTCAGAGATAATGCCTTTCACGTCACCATAGTGGGCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAg                       <  2:180634/250‑1 (MQ=255)                                                                                                                                                                                                                               ttCAGAGATAATGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTagag                                                                           <  1:89169/187‑1 (MQ=255)                                                                                                                                                                                                                                  agagATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGAt                                                                                                                      <  1:107334/141‑1 (MQ=255)                                                                                                                                                                                                                                  agagATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGAt                                                                                                                      <  1:134106/141‑1 (MQ=255)                                                                                                                                                                                                                                        aaTGCCTTTCACGTCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACt  <  1:313182/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              ACGATAGTTTTCGTCGGTGCAATAATGCTCGGGTCGCCCGTTGCCAGCGTCACGCCAAGGCCATTACCCACATTGCCTTGCGCGTTGTGCGGCAGTAACAGACCGTAAGAGTAGAGCGCACGGCGGCTCATGATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACAT‑CGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACT  >  NC_000913/4305792‑4306272

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