Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	663,503	C→T	S154N (AGC→AAC)	dacA ←	D‑alanyl‑D‑alanine carboxypeptidase DacA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	663,503	0	C	T	95.6%	122.0 / ‑2.3	45	S154N (AGC→AAC)	dacA	D‑alanyl‑D‑alanine carboxypeptidase DacA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (20/23);  total (21/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 = 6.28e-01

GGCTGTTATCCCATAACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGCTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTG  >  NC_000913/663258‑663739                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ggCTGTTATCCCAAAACTGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGAGCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGTCCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCTTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGGCCTGGTTACCa                                                                                                                                                                                                                                          <  1:93909/251‑1 (MQ=255)      ttATCCCATAACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGc                                                                                                                                                                                                                                     >  1:362716/1‑251 (MQ=255)                    gCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATgg                                                                                                                                                                                                                       >  1:519779/1‑251 (MQ=255)                         tACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGCTACCAGCGGCAAAATCGGCCATGGCGac                                                                                                                                                                                                                   <  1:102497/250‑1 (MQ=255)                           cGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGacac                                                                                                                                                                                                                 >  1:462242/1‑250 (MQ=255)                             gTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTTGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAg                                                                                                                                                                                                              <  1:558484/251‑1 (MQ=255)                                        gAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGtt                                                                                                                                                                                                       <  1:490401/247‑1 (MQ=255)                                        gAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGtt                                                                                                                                                                                                       >  2:490401/1‑247 (MQ=255)                                         aaTACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAAGCGGCCAGGGCGACACAAGCATCGTTACCCg                                                                                                                                                                                                  <  2:100332/251‑1 (MQ=255)                                                    aCGTCAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCATTGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCCTCAAGCCAACCAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTGCCCGATTTCAgggt                                                                                                                                                                                        >  2:408728/1‑248 (MQ=255)                                                             ctttttctttATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTAGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGCAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCgg                                                                                                                                                                              >  1:620381/1‑251 (MQ=255)                                                                tttctttATCGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCAGATAAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGCTACCAGCGGCAACATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAAGTTAATACCGCGGATc                                                                                                                                                                           >  2:174737/1‑251 (MQ=255)                                                                           aTCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAAc                                                                                                                                                                <  1:239429/251‑1 (MQ=255)                                                                             cGAGTATTCATTCGATACGTCACGGATCAATGTCTGACCGATCAGAGCCATAACTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTCCCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGCTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTTCCAGCGGCAAAATCGGCCTTGGCTACACAATCATCGTTCCCCGATTGCAGGTTAAGGCAGAGGATCAGCTgaggaaacg                                                                                                                                                              >  2:106762/1‑245 (MQ=255)                                                                                 tATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAAc                                                                                                                                                          >  2:199234/1‑251 (MQ=255)                                                                                     cattcGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACGGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAAGCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGTATCAGCTGAGAAACCGGAACCTGc                                                                                                                                                      <  2:293015/251‑1 (MQ=255)                                                                                                  ccGGATCAATGCCTGGCCGATCAGCGCCATAGCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTTCGTTAACGTAGCGGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGa                                                                                                                                         <  1:364719/250‑1 (MQ=255)                                                                                                        cAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAAc                                                                                                                                    <  1:606637/250‑1 (MQ=255)                                                                                                            gCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAg                                                                                                                               <  1:199234/251‑1 (MQ=255)                                                                                                             gctcgCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGc                                                                                                                              >  2:26240/5‑251 (MQ=255)                                                                                                             ccTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGGATGTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTCGAGGAACATCAGc                                                                                                                              <  2:46631/251‑1 (MQ=255)                                                                                                                  ccGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAACAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGGGGATCAGCTGCGCAACCGGAACCTGCATGCCCGGTTTGAGGCAAATCAGCgaagg                                                                                                                         >  2:502895/1‑250 (MQ=255)                                                                                                                          cgcCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTaa                                                                                                                 <  2:362716/251‑1 (MQ=255)                                                                                                                           gcCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTaaa                                                                                                                <  1:253394/251‑1 (MQ=255)                                                                                                                                                  cTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATg                                                                                        >  2:155048/1‑251 (MQ=255)                                                                                                                                                    gACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGcg                                                                                      >  1:533638/1‑251 (MQ=255)                                                                                                                                                    caCCATCAGCCTCCAGACCATGTACCGTCTGGAAGTGGGTATGTAGCATGCCCAGTGCGGTAACGTAGCTGGTCATCAAGCCAACAAAAGCGTCCTGTTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGGTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGTAACAGCAGCGAAGAACCTTTAAACGTCCGGGG‑ACCGGTGGCCCATgtg                                                                                      <  2:59540/250‑3 (MQ=255)                                                                                                                                                      ccATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGcgtc                                                                                    <  2:106524/251‑1 (MQ=255)                                                                                                                                                          cAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGGTACCGGTGGCCCATGCGTCGtt                                                                                 <  1:444757/250‑1 (MQ=255)                                                                                                                                                                ccAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCACCAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCCTGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGGTCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGAt                                                                           >  2:248713/1‑250 (MQ=255)                                                                                                                                                                ccAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCCTGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTAACCGGTGGCCCAGGcgtcgtgcccgaa                                                                           >  2:510600/1‑243 (MQ=255)                                                                                                                                                                 cAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCa                                                                                                                                                                                            <  1:179274/137‑1 (MQ=255)                                                                                                                                                                 cAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCa                                                                                                                                                                                            >  2:179274/1‑137 (MQ=255)                                                                                                                                                                               tCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCCGtt                                                           >  2:140647/1‑251 (MQ=255)                                                                                                                                                                                        gggTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTaaa                                                   <  2:272374/250‑1 (MQ=255)                                                                                                                                                                                          gTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCTGGTTACCGGTTGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAAttt                                                <  2:386024/251‑1 (MQ=255)                                                                                                                                                                                                  cAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGcgt                                        <  2:68348/251‑3 (MQ=255)                                                                                                                                                                                                    ggCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑TCCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTc                                      <  1:393256/251‑1 (MQ=255)                                                                                                                                                                                                                 aaCGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCg                          >  2:32377/1‑250 (MQ=255)                                                                                                                                                                                                                  aCGTAGCTGTTCATCAAGCCAACAAAAGCGGCCGGTTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCAGCGTTACCCGATTGCAGGGTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTGACCTGTGGCCCATGCGTCGTGGCCGATAGTTGCTAGATCATTTTCTTTAAATTTACCGGCTTTCACTTCCTGGCCGAt                        <  2:302184/251‑1 (MQ=255)                                                                                                                                                                                                                      aGCTGTTCATCAAGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAAc                     <  1:511452/250‑1 (MQ=255)                                                                                                                                                                                                                                  aGCCAACAAAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGtcatc        >  1:440621/1‑251 (MQ=255)                                                                                                                                                                                                                                   gcaaacaaAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACAAAAGCAGCGTTACCCGAGTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGTAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCAGTGCCTTGCCGATAACGTAACTGTtcatca       <  1:276455/248‑1 (MQ=255)                                                                                                                                                                                                                                   gCCAACAAACGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACTTAACTGGtcatc        >  2:144397/1‑250 (MQ=255)                                                                                                                                                                                                                                        caaAAGCGTCCTGGTTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAACTTTACCGGCTTTCATTGCCTGGCCGATCACGTAACTGGTCATCAttttt  >  2:506691/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              GGCTGTTATCCCATAACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCATTCGGTACGTCACGGATCAATGCCTGGCCGATCAGCGCCATATCTCGCGCGGAGCTGTACTGACCATCAGCATCCAGACCATGTACCGTCTGGAAGTGGGTATTTTTCAGGCCCAGTGCGTTAACGTAGCTGTTCATCAAGCCAACAAAAGCGTCCTGGCTACCAGCGGCAAAATCGGCCATGGCGACACAAGCATCGTTACCCGATTGCAGGTTAATACCGCGGATCAGCTGAGAAACCGGAACCTGCATGCCCGGTTTGAGGAACATCAGCGAAGAACCTTTAAAC‑ACCGGGTTACCGGTGGCCCATGCGTCGTTGCCGATAGTGACTAAATCAGTTTCTTTAAATTTACCGGCTTTCATTGCCTGGCCGATAACGTAACTGGTCATCATTTTG  >  NC_000913/663258‑663739

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