Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,937,789	G→A	A108T (GCA→ACA)	fucK →	L‑fuculokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,937,789	0	G	A	98.1%	151.7 / NA	53	A108T (GCA→ACA)	fucK	L‑fuculokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/21);  minor base C (0/1);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATGGCCGGAACCAGCCAGATGCTGGATATCC  >  NC_000913/2937546‑2938029                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 aTTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGt                                                                                                                                                                                                                                           >  1:302792/1‑251 (MQ=255)   tgttgCCCGCGCCTCAACGCCTAATGCCAGCGATAGCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTGAGACGCCATTTTGCAACGCTTTGCTGAGTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCTCCTTTGGGGTGGATG‑GCGCTGTGGTATATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACACCAGCGGtt                                                                                                                                                                                                                                          <  2:361371/250‑1 (MQ=255)       gCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATgg                                                                                                                                                                                                                                      <  2:429882/250‑1 (MQ=255)                         aaTGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAAt                                                                                                                                                                                                                   <  1:136179/251‑1 (MQ=255)                             ccAGCGATATCGCGATGGAAAACAACACCTGGCCCCAGTGGTCTATAGACGCCATTACGCAACGCTTTGCTGATTGCTGCCGGCACATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGCCCGCGAACAACAGCGATTATGGACAATATTGAACGGTTCATCTcc                                                                                                                                                                                                               >  1:261957/1‑251 (MQ=255)                             ccAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACATCGGTTATGGACAATATTGAACTGTTAATCTcc                                                                                                                                                                                                               >  2:186187/1‑251 (MQ=255)                                        gcgATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCTCATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGTTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTTGAAATGTCCGCGAACAACATCGGTTATGGACAATATTGAACGGTTAATCTCCGCCCAGCGGtt                                                                                                                                                                                                    <  2:369927/251‑1 (MQ=255)                                         cgATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATCGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGAAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCAGCGACCAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTg                                                                                                                                                                                                   >  1:441126/1‑251 (MQ=255)                                               aaaaCAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCt                                                                                                                                                                                             >  1:160392/1‑251 (MQ=255)                                                  acaacaCCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTAtt                                                                                                                                                                                          >  2:95017/1‑251 (MQ=255)                                                     acacCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCt                                                                                                                                                                                       <  2:73787/251‑1 (MQ=255)                                                      cacCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTg                                                                                                                                                                                      >  1:19900/1‑251 (MQ=255)                                                        ccTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGc                                                                                                                                                                                    >  1:302161/1‑251 (MQ=255)                                                           ggctccAGTTGTCTTTAGACTCCATTTTGCAACGCTTTGCTTATTGCTGTCGGCAAATCAATAGTGAACTGTCGTAATGCCACAGCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTc                                                                                                                                                                                 <  2:436212/247‑1 (MQ=255)                                                           ggCACCAGTTGTGTTTATACGCCATTTTGACACGCTTTGCTGATTGCTGGCGGCAAATTAATAGTGAACTGACTGAATGCCCCAGCCGCGGTTTCGCCGTCACCACCTGTGGTGTGGATG‑GCGCTCTGGTATATAAGCAAGGCATTCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTGATGGACATTATTGATCGGTTAATCTCCTCACTGCTGTTCCCGCCTATTTCTGGCGt                                                                                                                                                                                  <  2:455212/250‑1 (MQ=255)                                                                 aGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGAAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGcc                                                                                                                                                                           >  1:48411/1‑251 (MQ=255)                                                                       cTTTAGAAGCCATTTTGCAACGCTTTGCTGATTGCTCTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTCTGGTGTGGACG‑GCGCTCTGGTAGATAAGCACGGCAATCTGCTCTATCCGATTATTAGCTGGAACTGTCCCAGCACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGATGCAGGCTATTTCTGGCGTCGGAGCCTTTAGt                                                                                                                                                                     >  1:39399/1‑251 (MQ=255)                                                                         ttAGACGCCAGTTTGCAACGCTTTGGTGATGGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGGATCGCCGGCGCCGCCTTTGGTGTGGATG‑GCGCTCTGGTAGAGAAGCAAGGCAATCTGCTCTAGCCGAGGATTAGCGGGAAATGTCCGCGAACAACAGCGGTTATGTACGATATTGAACGGTGAAGCTCGGCGCAGCGGTTGCAGGGTATTTCTGGCGTCGGAGGCtttggt                                                                                                                                                                     <  2:176440/249‑4 (MQ=255)                                                                          tAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCg                                                                                                                                                                                                                                             <  1:44486/176‑1 (MQ=255)                                                                          tAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCg                                                                                                                                                                                                                                             >  2:44486/1‑176 (MQ=255)                                                                              cGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGAGTGGATG‑GAGCTCTGGTAGATAAGCAAGGCAATCTTCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTAGCAGGCTATTTCTGGCGTCGGAGCCTTTAGtttacat                                                                                                                                                               >  2:293612/1‑246 (MQ=255)                                                                                  aTTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCTCCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGtt                                                                                                                                                          <  2:132425/251‑1 (MQ=255)                                                                                             cattgctggtttggtgcGGCAAATCTATAGTGAAATTACTGAATGCCACCTCCGCGGTATCGCCGTCGCCACCTTTTGTGTTGATG‑GGGCTCTTGTAGATAAGCAATGCAATCTGCTCTAACCGAATATTAGCTTGAAATGGCCGCGAACAACATGGGTTATGGACAATATTTAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTgg                                                                                                                                               <  2:35310/235‑1 (MQ=255)                                                                                                    gATTGCTGTCGGCAAATCAATAGTGAACTGACTGCATGCCACATCAGCAGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGtt                                                                                                                                        >  1:200341/1‑251 (MQ=255)                                                                                                      ttGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGAGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGACATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATCCGTTATATAAGTTGGTGTGGTTg                                                                                                                                       >  2:91478/1‑250 (MQ=255)                                                                                                            tCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTgaaagaa                                                                                                                                 >  1:374031/1‑250 (MQ=255)                                                                                                                 aaaTCAATAGTGAACTGACTGAATGCCACATCCGCGGCATCACCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGCCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAatcat                                                                                                                           >  1:67443/1‑251 (MQ=255)                                                                                                                          gTGAACTGGCTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTGTGTGTGGATG‑TCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCATTAGTTTCAATACGTTATATAAGTTCGTGTGGTTGAAAGAAAAACAGCCACAActa                                                                                                                  <  2:43290/251‑2 (MQ=255)                                                                                                                              actgactgAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGTTGTGGTTGAAAGAAAATCATCCACAACTGCTgg                                                                                                              >  2:391932/1‑251 (MQ=255)                                                                                                                                     gaATGCCACATCCGCGGTATAGCCGTCACCACATTTCGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGAAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACCGCGGTTATGGAAAATATTGAACGGTTAATCTACGCACAGCGGTTGCAGACTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTAAATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGATGGCAcgcgc                                                                                                       >  2:250990/1‑251 (MQ=255)                                                                                                                                      aaTGCCACATCCGCGGTATCGCCGTCACCCCCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAAcgcgcg                                                                                                      <  2:9768/251‑1 (MQ=255)                                                                                                                                        tGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACgcgcgc                                                                                                     >  1:109678/1‑250 (MQ=255)                                                                                                                                        tGCCACATCAGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCCATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACgcgcgc                                                                                                     >  2:453597/1‑250 (MQ=255)                                                                                                                                                ccGCGGTATCGCCTTCACCGCCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCTGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCAGCCACAACTGCTGGAACGCGCGCACGCCTGGc                                                                                            <  1:223618/251‑1 (MQ=255)                                                                                                                                                ccGCGGTATCGCCGTCACCACCTTTGGTGTTGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGTCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGc                                                                                            <  1:376151/251‑1 (MQ=255)                                                                                                                                                ccGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGc                                                                                            <  1:208779/251‑1 (MQ=255)                                                                                                                                                     gtaacGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCttt                                                                                       <  1:90868/247‑1 (MQ=255)                                                                                                                                                         cGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTAttt                                                                                   >  1:388196/1‑251 (MQ=255)                                                                                                                                                               caccacCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTtcgtcg                                                                              <  2:142775/250‑1 (MQ=255)                                                                                                                                                                accaccTTTGGTG‑GGGTGTGCGATTTGGTAGAGAAGCAAGGCAATTTGCTATATCCGATTATTAGCTGGAAATGGCCGCGAACAACAGCGGTTATGTACAATATTGAACGGTTAATCTCCGCACAGCGGTTGGATGCGATTTCTGGCGTCGGAGCCTTGAGTTTCAATACGTTATAAATGTTTGTGTGGTTGAAAGATAATCATCCGCAACTTCTGGTGCGCGCGCACGACTGTCTCTTTAGCTCGTCGCt                                                                            <  2:423418/251‑1 (MQ=255)                                                                                                                                                                                        tctGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCCTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCg                                                     >  2:153124/1‑251 (MQ=255)                                                                                                                                                                                          tGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGACCAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGa                                                    >  1:218615/1‑250 (MQ=255)                                                                                                                                                                                               gATAAGCGAGGCAAACTGCTCTCTCCGATTATTAGCTGGAAATTTCCGAGAACAACAGCGGTTATGGACCTTATTGAACGGATAATCTACGCACAGCGGTTGCAGGCTATTGCTGGCGTAGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGAGTGGTTGAAAGAAAATCATCAACAACTGCTGTAACGCGCGCACACCTGGAGCTTTATCTCGTCACTGATTAACCACCGTTTAACCGACGAATTccc                                              >  2:414645/1‑249 (MQ=255)                                                                                                                                                                                               gATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTCTAACCGGCGAATTCac                                              >  2:414672/1‑251 (MQ=255)                                                                                                                                                                                                           aaccTGCTCTACCCGATGATTAGCTGGAAATGTCCTCGAACAACAGCGGTGATGGACAATATTGAACTGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTTGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCTTCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCTTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCAc                                  <  1:187941/248‑1 (MQ=255)                                                                                                                                                                                                           aaTCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCAc                                  >  1:231771/1‑251 (MQ=255)                                                                                                                                                                                                                  tctATCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATGGcc                           >  1:20947/1‑251 (MQ=255)                                                                                                                                                                                                                   ctATCCGATTATTAGCTGGAAATGTCCGAGCACAACAGCGCTTATGGCCCATATTGAACGGTTAATCTCCGCACAGCGGTTCCAGGCTATTTCTGGCGCCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGAGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGCTTAAACACCGTTTAACCGGCGAATTCACTACTGATATCACggatgccg                          >  2:438768/1‑244 (MQ=255)                                                                                                                                                                                                                     aTCCGATTATTAGCTGGAAATGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGTTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATGGcc                           <  1:309664/248‑1 (MQ=255)                                                                                                                                                                                                                          attattAGCTGGAAATGGCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGTTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATGGCCGGAAcc                     <  2:302161/249‑1 (MQ=255)                                                                                                                                                                                                                                        aTGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGATATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGCATTCACTACTGATATCACGATGGCCGGAACCAGCCAGATGCTGGat      >  1:229466/1‑250 (MQ=255)                                                                                                                                                                                                                                        aTGTCCGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGATGCAGGCTATTTCTGGCGTCGGAGCCTTTAGATTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATTGCCGGAACCAGCCAGATGCTGGat      >  1:366214/1‑250 (MQ=255)                                                                                                                                                                                                                                            ccGCGAACAACAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATGGCCGGAACCAGCCAGATGCTGGATATcc  >  1:220774/1‑250 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 ATTGTTGCCCGCGCCTCAACGCCTAATGCCAGCGATATCGCGATGGAAAACAACACCTGGCACCAGTGGTCTTTAGACGCCATTTTGCAACGCTTTGCTGATTGCTGTCGGCAAATCAATAGTGAACTGACTGAATGCCACATCCGCGGTATCGCCGTCACCACCTTTGGTGTGGATG‑GCGCTCTGGTAGATAAGCAAGGCAATCTGCTCTATCCGATTATTAGCTGGAAATGTCCGCGAACAGCAGCGGTTATGGACAATATTGAACGGTTAATCTCCGCACAGCGGTTGCAGGCTATTTCTGGCGTCGGAGCCTTTAGTTTCAATACGTTATATAAGTTGGTGTGGTTGAAAGAAAATCATCCACAACTGCTGGAACGCGCGCACGCCTGGCTCTTTATTTCGTCGCTGATTAACCACCGTTTAACCGGCGAATTCACTACTGATATCACGATGGCCGGAACCAGCCAGATGCTGGATATCC  >  NC_000913/2937546‑2938029

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