Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,814,459	Δ1 bp	coding (275/516 nt)	luxS ←	S‑ribosylhomocysteine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,814,459	0	A	.	95.7%	187.9 / 5.0	47	coding (275/516 nt)	luxS	S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (21/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

GATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTAT  >  NC_000913/2814221‑2814702                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    gATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCatat                                                                                                                                                                                                                                         >  1:30739/1‑250 (MQ=255)      gccgTTCCTGCCACTTCTCGTTCGGCTGTGCCAGTTCCTCGTTGCTGATGATGCGTACGTCCTGTTCCAGAATGCGACGCGCAATATCCTGCGCTTCCTGCAACGTGTGCAGCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCGTCCATTGCCGCTTTCCAGGCATCAGCAACACGCGGCTCATCTGGCGTAC‑AATTAGACTCATATAAAAAc                                                                                                                                                                                                                                   <  2:41991/248‑1 (MQ=255)                atacttattttgcGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTAACGTTACAGAATGCTACGCGCAATAGCCTGCGCTTCCTGAACCGAGTGTATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTGCAGCACGTCTTCTATTGCCGCTTTCCAGGCATCAGCAACACGCGGCTCATCTGGAGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCa                                                                                                                                                                                                                         <  1:318146/239‑1 (MQ=255)                 aaCTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATTCGTACGTCACGTTCCAGAATGCTACGCGCAATAGCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAg                                                                                                                                                                                                                        <  2:480559/251‑1 (MQ=255)                          gtctgCAGTGCCGGTTCTTCGTTGCGGTTGATGCGTACGTCACGTTCCTGAATGCTACGCGCAATAGCCTGCTCTTCCTGCACCGAGTGCATCGTGTAAGTGCCACACTGGTCGACTTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGGCTTCCATTGCCGCTTTCCTGGCATCAGCAGCTCGCTGCTCTTCTGGCGTTCCA‑TTAGACTCATATAAAAGCCAGTGCGGCAGCCCATTgg                                                                                                                                                                                                                <  2:6904/246‑1 (MQ=255)                                aGTGCAAGTTATTCGTTGCTGTTGATGAGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCCCCGAGTGCGTCTGGTAAGGGCCACACTGGTAGACGTTCAGTTCCGAGATCTGATTATGATCCTACACATGAAGCACGTATTCAATTGCCGATTTCCAGACTTCAGCAAAACACTGCTCATCTGGAGTCCCATGCA‑CCTCATATAAAAAcatacgcggcagcacattcgcgcgaa                                                                                                                                                                                                          >  2:336425/1‑225 (MQ=255)                                     cAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGAt                                                                                                                                                                                                    >  1:451165/1‑251 (MQ=255)                                      aGTTCTGCGTTGCCTTTGATGCGTACGTCACGGTCCCGATGGCTACGCGCAATATCCTGCGCGTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCAGGTGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCTTTGGCGAGATATCGATa                                                                                                                                                                                                   <  2:459858/251‑1 (MQ=255)                                             cGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATctct                                                                                                                                                                                             <  1:26373/250‑1 (MQ=255)                                             cGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGACTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGCCAGCCCATTGGCGAGATATCGATAATctct                                                                                                                                                                                             >  2:469025/1‑250 (MQ=255)                                                    gTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCAtt                                                                                                                                                                                     >  1:19309/1‑251 (MQ=255)                                                    gTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGATCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCCTTGGCGAGATATAGATAATCTCTACACCAtt                                                                                                                                                                                     >  1:42214/1‑251 (MQ=255)                                                     ttGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCAtt                                                                                                                                                                                     >  2:148052/1‑250 (MQ=255)                                                                                    gcgcAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTa                                                                                                                                                                                            >  1:64413/1‑212 (MQ=255)                                                                                       cAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAAc                                                                                                                                                   >  1:19859/1‑250 (MQ=255)                                                                                                  gcTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGATTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAgg                                                                                                                                       <  1:86870/251‑1 (MQ=255)                                                                                                      ccTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTAt                                                                                                                                   >  2:431714/1‑251 (MQ=255)                                                                                                           aaCGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATc                                                                                                                              >  1:232050/1‑251 (MQ=255)                                                                                                              gAGTCCATCTGGTACGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATAACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCGGGGTATGGATCCct                                                                                                                           >  2:305464/1‑251 (MQ=255)                                                                                                                   cATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGtt                                                                                                                                                                      <  1:115955/203‑1 (MQ=255)                                                                                                                   cATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGtt                                                                                                                                                                      >  2:115955/1‑203 (MQ=255)                                                                                                                     tcgggTAAGTGCCACACTGGAAGACGTTACGTTCCCGGAACTGCTTCTGCTCCTGCACTTTCAGAACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCGGCTCATCTGGCGTACCA‑TTAGACTCATATACAAACCGGTGCGGCAGCCCAGAGGCGAGATATCGCTAATCTCTAAACCCATACCGTTAAGGTGGCTCCCCATCTAAACAGCAAACAGGAGCTCCTGGGAATGGATCCCTCTTTCTg                                                                                                                    >  2:554655/4‑251 (MQ=255)                                                                                                                     tCTGGTAAGGGCCCCACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCt                                                                                                                     <  2:30739/250‑1 (MQ=255)                                                                                                                         gTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCAt                                                                                                                >  1:368349/1‑251 (MQ=255)                                                                                                                            aGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATTATCTCTACACAATTACCGTTAAGATGGTTACGCATAAAACCAGCAAAAAGGTGCACAAGGGTATGGATCCCTCATTCTGGAATCa                                                                                                              >  2:397795/1‑250 (MQ=255)                                                                                                                             gTGCCACACTGGTCGACGTTCAGTTCCGGGATCTGATTCTGCTCCTGCCCTTTCAGCACGTCTTCAATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGAGCGGCAGCCCATTGGCGAGCTATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAAACAGCAAAAAGGTCCTCCAGGGTATGGATACCTCTTTCTGGCATCaat                                                                                                            >  2:34435/1‑249 (MQ=255)                                                                                                                                  acacTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCttt                                                                                                       <  1:431714/251‑1 (MQ=255)                                                                                                                                      gggTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTc                                                                                                   <  1:260794/250‑1 (MQ=255)                                                                                                                                            aCGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCAGCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACg                                                                                             <  2:19309/251‑1 (MQ=255)                                                                                                                                                tCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAg                                                                                          <  1:453044/250‑1 (MQ=255)                                                                                                                                                   gTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGATCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAgc                                                                                      >  1:509206/1‑251 (MQ=255)                                                                                                                                                        gggATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCCCGCAGAAgcgc                                                                                    >  2:111250/1‑248 (MQ=255)                                                                                                                                                                                    gCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGc                                                      >  2:89987/1‑251 (MQ=255)                                                                                                                                                                                       cGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCGGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCgggg                                                  <  1:305464/251‑1 (MQ=255)                                                                                                                                                                                                 tGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGTTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCTTCGCCATGCGTGGTGTTCATTGt                                        >  2:89992/1‑251 (MQ=255)                                                                                                                                                                                                     gatttCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACAATTACCGTTAAGCTGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGAACGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTc                                    >  2:343565/3‑251 (MQ=255)                                                                                                                                                                                                        ttCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCCGAAGCGCAGATCGAACACGATGATTGCGTCGCCATGCGGGGGGTTCATTGTTTTCGcc                                 >  2:356456/1‑251 (MQ=255)                                                                                                                                                                                                         tCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCAGTGGCGAGAGATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTTTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGcc                                 <  1:554655/250‑1 (MQ=255)                                                                                                                                                                                                                aTCAGCAACACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCATCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGGTACGCATAAATCCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCGTTGTGCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCTCCCGAAc                         <  1:353294/251‑1 (MQ=255)                                                                                                                                                                                                                      aaCACGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAgg                   <  1:469025/251‑1 (MQ=255)                                                                                                                                                                                                                       acacGCTGCTCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGgc                  <  1:356456/251‑1 (MQ=255)                                                                                                                                                                                                                           gctgctCATCTTGCGGTCCA‑TTAGACCCATATTAAAACCGGTGCTTCAGCCCATTGGCGAGATATCGATAATTTCTACACTATTACCGTTAAGGTGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCGGGCATCACATCTTTGTTCGGCACGCAGAAGCGCAGATCGAGCACGGTGATTGCGTCTCCATGCGGGGGGGTCATTGGTTTCGCCGCCCGTACTGCAGGCGCtt              <  2:232050/251‑1 (MQ=255)                                                                                                                                                                                                                           gctgctCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCtt              <  1:148052/251‑1 (MQ=255)                                                                                                                                                                                                                              gctCATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTcc            <  1:100969/250‑1 (MQ=255)                                                                                                                                                                                                                              gcacATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTc             <  2:451165/246‑1 (MQ=255)                                                                                                                                                                                                                                 cATCTGGCGTACCA‑TTAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACTCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCAt          <  1:39874/249‑1 (MQ=255)                                                                                                                                                                                                                                       gCGTACCA‑TTAGACTCAAATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCACCAGGGTATGGATCCCTCTGTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTAt  <  2:147623/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    GATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATCCTGCGCTTCCTGCAACGAGTGCATCTGGTAAGTGCCACACTGGTAGACGTTCAGTTCCGGGATCTGATTCTGATCCTGCACTTTCAGCACGTCTTCCATTGCCGCTTTCCAGGCATCAGCAACACGCTGCTCATCTGGCGTACCAATCAGACTCATATAAAAACCGGTGCGGCAGCCCATTGGCGAGATATCGATAATCTCTACACCATTACCGTTAAGATGGTTACGCATAAAACCAGCAAACAGGTGCTCCAGGGTATGGATCCCTCTTTCTGGCATCACTTCTTTGTTCGGCACGCAGAAGCGCAGATCGAACACGGTGATTGCGTCGCCATGCGGGGTGTTCATTGTTTTCGCCACCCGAACTGCAGGCGCTTCCATCCGGGTAT  >  NC_000913/2814221‑2814702

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