Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,483	T→A	S7S (TCT→TCA)	cysZ →	sulfate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,483	0	T	A	97.4%	131.7 / NA	39	S7S (TCT→TCA)	cysZ	sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (24/14);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGC  >  NC_000913/2531240‑2531729                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       ctctgttgatctctaattttcccaacacGGCAAACATTTAATCACTAAGAGTATTTACCGATTACCTCAAGTGCAAGTGCACTATTAACTTTAACAGCACAAAGGTAGATGACATCGAGCTTTTTGCTGGCTTTTCGAACATATCCTAACTGTCCATTGCGAAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATACTGATTATGGTTTCATCCTTCACATCAGCCCCa                                                                                                                                                                                                                                                  >  2:403455/26‑250 (MQ=255) ctctGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGCTTACCTCAAGTGCAAGTGCCCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGTGATTGGTCACAAGGAGCTCATCCTGACTATGGTTTCATCATTCACATCAGCCCCAc                                                                                                                                                                                                                                                 >  2:19476/1‑251 (MQ=255)                         ccacGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAAGCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAATTATCTCATCCTGATTAGGTTTTTATCTTTCACATCAGCCCCACGTAGCTGGTTTTACTATTTTgcgc                                                                                                                                                                                                                         <  2:270039/251‑1 (MQ=255)                                    aTTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGAACTATTAACTTTCACAGCCCAAAGATAGATGACATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCAGGTACAATCTGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGTCCCACGCAGCGGTTTTTACTATTTTGCGCAAgtatgga                                                                                                                                                                                                                >  2:59925/1‑243 (MQ=255)                                        aCTCACTAAGAGTATGTGCCGATTACCGCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTGTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGTTTTCAGCATTCACATCAGCCCCACGGAGCGGTTTTTACTAGATGGCGCAAGGCTGGAAGCTc                                                                                                                                                                                                           <  2:409855/250‑1 (MQ=255)                                            aCTAAGAGTATTTGCCGAGTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCCCAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCGAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGTTCATCCGGATTATGGTTTCATCATTCACAGCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGtctc                                                                                                                                                                                                      <  1:366302/251‑1 (MQ=255)                                                 gagTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGAAGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCacc                                                                                                                                                                                                 >  2:167104/1‑249 (MQ=255)                                                          ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTc                                                                                                                                                                                        <  1:301737/251‑1 (MQ=255)                                                                  cTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTc                                                                                                                                                                                >  1:307523/1‑251 (MQ=255)                                                                     aagtgcaagtgcACTATTAACTTTCACAGCACAAAGATAGATGAAATAGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGtt                                                                                                                                                                             >  2:361753/1‑251 (MQ=255)                                                                       gtgcaagtgcaCTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGATTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCCTTGCGCAATTACCCGGTAACATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTAt                                                                                                                                                                           >  1:221363/1‑251 (MQ=255)                                                                                    aTTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCgctgct                                                                                                                                                               <  1:340511/250‑1 (MQ=255)                                                                                         cTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCa                                                                                                                                                          >  1:87863/1‑250 (MQ=255)                                                                                         cTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCa                                                                                                                                                          >  1:75327/1‑250 (MQ=255)                                                                                              acaGCAAAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTTGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTCTGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCAAGTCTCGCAACCTGGGATTCGGCGTTTCGTTCTTTTACCGCTGCTGGTCAATAttt                                                                                                                                                    >  1:298287/1‑251 (MQ=255)                                                                                                                 aaTCGTGCTTTTTGCTGTTTTTTAGATCATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCACAAGGAGCCCATCCTGATTATGGTTTCATCCTTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCACGTCTCGCAACCTGGGATTCGGCGTTTAGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCAtt                                                                                                                                 >  2:154336/1‑251 (MQ=255)                                                                                                                  aTCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCACAAGGAGCTCATCCTGCATATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGACAATATTTTGTTGATGGGGGGCGCATTc                                                                                                                                >  1:272240/1‑251 (MQ=255)                                                                                                                  aTCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTGTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTc                                                                                                                                <  1:361753/251‑1 (MQ=255)                                                                                                                              gCTGTTTTTTCGAACATATCCTAACTGTCAATTGCGCCATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGCTGGGGGTCGCATTCTGGTGGCTCttt                                                                                                                    >  1:140389/1‑251 (MQ=255)                                                                                                                                tgtattctCGAACATATCCTAACTGTCCATTGCGCAAATACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGATAATCCTGATTATGGTTTAATCCTTCACATCAGCCCCACGAAGCGGTTTTTACTATTTTGCGCACGGCTGGAAGCTCGTCTCACAACCTGGGATTCGGCGTTTCGTTATTTTACCGCGGCTGTTCAATATCTTGTTGATGGGGGGCGCATTCTGTGGGCTCttaac                                                                                                                  >  2:284487/8‑248 (MQ=255)                                                                                                                                               atCCTAACTGTCCCTTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                                                                   >  2:91992/1‑251 (MQ=255)                                                                                                                                               atCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                                                                   >  2:268376/1‑251 (MQ=255)                                                                                                                                                         tCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGCTGTCTGGATCCCGAct                                                                                         >  2:142308/1‑251 (MQ=255)                                                                                                                                                           cATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTAGGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTAtttt                                                                                                                                                                                                                             <  1:72153/117‑1 (MQ=255)                                                                                                                                                           cATTGCGCAATTACCAGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTAtttt                                                                                                                                                                                                                             >  2:72153/1‑117 (MQ=255)                                                                                                                                                                   aaTTACCCGGTAAAATACGCAGAATTTTCCTGGGCTTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCCTGAGtt                                                                               >  2:85896/1‑251 (MQ=255)                                                                                                                                                                          cGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGGTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTcc                                                                        <  1:268376/251‑1 (MQ=255)                                                                                                                                                                              aaaaTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGAc                                                                    >  1:15896/1‑251 (MQ=255)                                                                                                                                                                               aaaTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACt                                                                   <  1:167104/251‑1 (MQ=255)                                                                                                                                                                               aaaTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACt                                                                   <  1:82494/251‑1 (MQ=255)                                                                                                                                                                                   aCGCAGAATTTTCCTGGGATTGGTCAAAAGGAGATCATCCTTATTATGGTTTCATCAGTCACATCAGCCCCACGCAGCTGTTTTTGCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCGCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCt                                                               <  2:372884/251‑1 (MQ=255)                                                                                                                                                                                       aGAATTTTCCTGCGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACa                                                            >  1:357784/1‑250 (MQ=255)                                                                                                                                                                                          aTTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTGCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGACTCGGCTTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCACATTCTGGTGGCGCTTTACACAGCCCGAGGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATgg                                                        <  1:397005/251‑1 (MQ=255)                                                                                                                                                                                                                   gCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCCGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGTCGGAGTTATCTGTTGTCGCGTCTgg                               >  2:166393/1‑251 (MQ=255)                                                                                                                                                                                                                      cATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTg                            >  1:333311/1‑251 (MQ=255)                                                                                                                                                                                                                            gATTATGTTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCTGACTGGCTACAATGGCTGAGTTATCTGGTGTGGCCTCTGGCTGTCAtc                       <  1:303157/250‑1 (MQ=255)                                                                                                                                                                                                                                              cacaTCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGATTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGttt     <  1:278476/250‑1 (MQ=255)                                                                                                                                                                                                                                              cacaTCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTg    <  2:15896/251‑1 (MQ=255)                                                                                                                                                                                                                                                caTCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTAGTCTCGCAACCTGGGATTAGGCGTTTCGTTATTTTCCCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGAGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCCCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGc  >  1:37882/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       CTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGC  >  NC_000913/2531240‑2531729

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