Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,256,301	G→A	S117S (TCC→TCT)	cyuA ←	putative L‑cysteine desulfidase CyuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,256,301	0	G	A	97.6%	118.8 / ‑6.7	41	S117S (TCC→TCT)	cyuA	putative L‑cysteine desulfidase CyuA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (22/18);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGC  >  NC_000913/3256057‑3256528                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    cAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACt                                                                                                                                                                                                                               >  1:480044/1‑251 (MQ=255)       ggACTTCATTGACGAACTTCAGGATCTCAGCCCGCGTCGTGCGGGAAAGCACCGTCAGCGTAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTcccc                                                                                                                                                                                                                         <  1:410120/251‑1 (MQ=255)                   cGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGt                                                                                                                                                                                                             >  2:134671/1‑251 (MQ=255)                   agaatttCAGGATCACACCCAGCGGCGTTCGGGAAAGCCCCGGCAGCGGAGACCCTTGCTCGCCATCTGCCCCACACGCCTGCGGGGAAAACCCCCCCCCAGCGGGCGTCTCGATATGCAC‑AATGGTGGGAGGCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGt                                                                                                                                                                                                             <  2:412483/246‑1 (MQ=255)                      aCTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCCGCAGTGCt                                                                                                                                                                                                          >  2:43449/1‑251 (MQ=255)                         tcgtgatcacccccagcTTCGTCCTTGAAACCACCGACAGGGGCGCCTCTTGCCCGCCCCCTTCCACCCCCCCCTTCGGGGTAAACCCCACCCCACGGTGCGTCTCGATACGCAC‑AATGTTGGTATGCCCCCCGACGATGGGGACCCACGCCCTCTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGACCGTAACAGATACTTTCCCCGCCGCCAGCAGTg                                                                                                                                                                                                            <  2:285994/234‑1 (MQ=255)                         caaggatgtcagcagccgtcgtTCTGGAAAGCACCGTCAGGGGAGCCTCTTGCCCGCCCTCTCCCACCCACCCCGGCTGGTTAAACACCACACCATCGTGCGTCTCGATATGCAA‑AATTTTGGTATGCCCGCCGACGAGGGTGACACACGCCCACTTCCCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCATGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCtt                                                                                                                                                                                                         <  1:197249/233‑1 (MQ=255)                             gATCTCAGCCAGCGTCGTTCTGGACAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCCAGGTTCCTGGATCTTAACAGAGACTTTCCCCCCCGCCAGCCGTGCTTTGGCAt                                                                                                                                                                                                   >  2:484315/1‑251 (MQ=255)                              agctcagccggctgcgtTCTGGAAAGCCCCGTCAGCGGATCCTCTTGCTCGCCCCCTGCCGCTCACGCCTGCTGGTTAACCACCACCCCATCTTGCGTCGCGATATGCTC‑AATGGTGTTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATc                                                                                                                                                                                                  <  1:450194/237‑1 (MQ=255)                                  cagccagcGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCCCACCACCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCa                                                                                                                                                                                              <  1:134671/251‑1 (MQ=255)                                         cgtcgtTCTGGAAAGCACCGTCAACGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCCTCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCAGAAATTGCCt                                                                                                                                                                                       >  2:87508/1‑251 (MQ=255)                                                          ccGTCAGCGGAGCCTTTTGCCCCCCCTCGGCCCAACACGCCGGCGGGGAAAACCCCACCCCACGGGGCGTCTCGATATGCAA‑AATGTGGGTATGCCCGCCGACGATGGTGACACACGCCCACTGCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTc                                                                                                                                                                      <  2:164329/251‑1 (MQ=255)                                                           cGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACCCACGCCTGCTGGGTAAACACCACCCCATCGTGCGTCTCGATAGGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCa                                                                                                                                                                     <  1:362233/251‑1 (MQ=255)                                                             ccAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAA‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGc                                                                                                                                                                   <  1:131893/250‑1 (MQ=255)                                                                   gagaCTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTcc                                                                                                                                                             >  2:62855/1‑251 (MQ=255)                                                                            gCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCACTTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCg                                                                                                                                                    >  1:278040/1‑251 (MQ=255)                                                                                      cgcccccctcgctgggtgGGTAACCACCACCCCATCGTGCATCCCATTATGCAC‑AATGTGGGTATGCCCGCCGACGAGGGTGACACACGCCCACTTCTCACCGTTCCAGACTTGGGCGCGGGAAAAGAGGATTTCCACGCAAGGTTCCTGGATCTTAACAGAGCCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTc                                                                                                                                          <  2:113668/235‑1 (MQ=255)                                                                                           cccccGCCTGCTGGGTAACCCCCACACCACCGTGCTTCTCGAT‑TGCCCCAAGTTTTGTATCCCCGCCGACGATGGTGACACACGCCCACTTCTCCCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCACCGCAAGGTTCCTGGATCTTAACAGATCCTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCt                                                                                                                                     <  1:215096/247‑1 (MQ=255)                                                                                              acGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAAc                                                                                                                                  >  2:150592/1‑251 (MQ=255)                                                                                                      tgGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGccccc                                                                                                                            >  2:353414/1‑249 (MQ=255)                                                                                                         gTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAgcgc                                                                                                                       >  2:158017/1‑251 (MQ=255)                                                                                                            aaCACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGcgccgc                                                                                                                    >  1:94983/1‑251 (MQ=255)                                                                                                                 cacaccaTCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACCCACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCt                                                                                                                                                                         >  1:170385/1‑193 (MQ=255)                                                                                                                      aaTCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGATGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCCTTTCCACCTAACGCCCCCAGCGCCGCCGCCATCGGc                                                                                                          >  2:170805/2‑251 (MQ=255)                                                                                                                             cGTCTCGATATGCAC‑AATGTTGGAATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCATCACTTCCAGCCCGGCGTTGGCATTTCCACCTCACGCCCACAGCGCCGCCGCAATCGGCAGcccc                                                                                                    >  2:407888/1‑250 (MQ=255)                                                                                                                              gTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATc                                                                                                                                                                                                  >  2:226553/1‑155 (MQ=255)                                                                                                                              gTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATc                                                                                                                                                                                                  <  1:226553/155‑1 (MQ=255)                                                                                                                                ctcGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCccacc                                                                                                 >  1:199584/1‑250 (MQ=255)                                                                                                                                           c‑aATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCGCCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCg                                                                                     >  1:375548/1‑251 (MQ=255)                                                                                                                                                        ccccGCCTACGATGGTGACACCCCCCCGTTTCTCCCCGTTCCAGACTTTGGCGCGTTAGAAGAGGTTTTCATCGCAAGGTGCTTGGATTGTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGTTGTTGCGGCTTTCAGCACTTCCAGCCCGGCGTTTGCATTTCCACCTAGCGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATGCCCGTGCCGGGAACGGTGACg                                                                         <  1:457966/250‑1 (MQ=255)                                                                                                                                                                acggtgtggacccacGCCCCCTCCCCACCGTCCCAGCCTTTGGCGCGTGAGAAAAGGTGTTCATCCCAAGGTTCTGGGATCTTAACGAGGTCTTCCCCCGCCGCCAGCAGTGCTTTGGCACCGGCAATTGCCTGCGCTGTTGCGGCTTTCAGCACTTCCAGCCCGGCGTTGGCATTGCCCCCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGAcc                                                                 <  2:99852/239‑1 (MQ=255)                                                                                                                                                                   atggggacccccccccACTTCCCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGAAAGGTTCCTGGATCTTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGtt                                                              <  1:317603/238‑1 (MQ=255)                                                                                                                                                                        gACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtca                                                         <  1:150592/251‑1 (MQ=255)                                                                                                                                                                           acacGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTtcatc                                                       >  1:475728/1‑250 (MQ=255)                                                                                                                                                                              cGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGGGACGCCCAGACCGTTCTTCATCAgct                                                   >  1:176900/1‑249 (MQ=255)                                                                                                                                                                                    cTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCg                                             >  1:299797/1‑251 (MQ=255)                                                                                                                                                                                      tctcACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACAGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGa                                            >  2:35006/1‑250 (MQ=255)                                                                                                                                                                                                           cgcgTGAGAAGAGGATTTCATCGCAAGGTTACTGGATCTTCACAGAGACTTACCCCGCCGCCAGCAGTGCTTTGGCATCGGCACTTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCCTTTCAACCTAACGCCCCCAGCGCCGCCGCCATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGCCCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTc                       >  2:301210/1‑250 (MQ=255)                                                                                                                                                                                                                   aagaGGATTTCACCGCAAGGTTCCTGGTTCTTAACAGAGTCTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGAcc              <  1:158017/251‑1 (MQ=255)                                                                                                                                                                                                                     ggggattccaccgaagggTTCCGGGATTGTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCtt            <  2:69974/235‑1 (MQ=255)                                                                                                                                                                                                                               ccGCAAGGTTCCTGGATCTTAACAGAGACTTCCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGc  <  1:234062/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                    CAAACGGGACTTCATTGACGAACTTCAGGATCTCAGCCAGCGTCGTTCTGGAAAGCACCGTCAGCGGAGACTCTTGCTCGCCCTCTGCCACACACGCCTGCTGGGTAAACACCACACCATCGTGCGTCTCGATATGCAC‑AATGTTGGTATGCCCGCCGACGATGGTGACACACGCCCACTTCTCACCGTTCCAGACTTTGGCGCGTGAGAAGAGGATTTCATCGCAAGGTTCCTGGATCTTAACGGAGACTTTCCCCGCCGCCAGCAGTGCTTTGGCATCGGCAATTGCCTGCGCTGTTGCGTCTTTCAGCACTTCCAGCCCGGCGTTGGCATTTCCACCTAACGCCCCCAGCGCCGCCGCAATCGGCAGCCCCACCATTCCCGTGCCGGGAACGGTGACGCCCAGACCGTTCTTCATCAGATTTGGCGAAACCCAGGCTTCTACACGTTCAACCGGACCTTCCAGTTCTGC  >  NC_000913/3256057‑3256528

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