Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,625	G→A	S70F (TCC→TTC)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,625	0	G	A	98.7%	251.0 / NA	76	S70F (TCC→TTC)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/44);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCG  >  NC_000913/3176379‑3176868                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    gCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTg                                                                                                                                                                                                                                                 <  2:195744/251‑1 (MQ=255)   cattttACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACTCCTGGCTATCCAGCAACATGATTTGCCATTGCTCACCAATAAACGCGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAAGCGTGGTGGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGAt                                                                                                                                                                                                                                               <  1:262810/248‑1 (MQ=255)    aGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTAATc                                                                                                                                                                                                                                              >  1:485887/1‑251 (MQ=255)        tACGTTACAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCCAACACCTGGATATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAACACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTgc                                                                                                                                                                                                                                          >  2:57228/1‑251 (MQ=255)         aCGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTgc                                                                                                                                                                                                                                          <  1:33626/250‑1 (MQ=255)         aCGTTCCAGCCACTCAAGCTGAAACTCGCGCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTgcg                                                                                                                                                                                                                                         <  1:541549/251‑1 (MQ=255)          cGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTgcgc                                                                                                                                                                                                                                        >  1:75059/1‑251 (MQ=255)           gTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTTCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCt                                                                                                                                                                                                                                       <  1:134916/251‑1 (MQ=255)              ccAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTaaa                                                                                                                                                                                                                                    <  1:478740/251‑1 (MQ=255)               cAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTa                                                                                                                                                                                                                                      >  1:286760/1‑248 (MQ=255)               cAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACGTGCGATGCCTTCAGCGAAATGCTGAGAGGCCGCAGAGAATTGATCCTGCGCTa                                                                                                                                                                                                                                      <  2:286760/248‑1 (MQ=255)               cAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGCCCCGCCAACCACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACAACGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAAt                                                                                                                                                                                                                                   >  2:500583/1‑251 (MQ=255)                    aCTCAAGCGGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCt                                                                                                                                                                                                                              <  2:433301/251‑1 (MQ=255)                        aaGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAACCACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATAACCTGTcg                                                                                                                                                                                                                          >  2:398729/1‑251 (MQ=255)                         aGCTGAAAATCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACCGGATTTGCCATTGCTCACAAAAAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCCGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGCTCCTGCGCTAACTCACCTGTcgc                                                                                                                                                                                                                         >  1:511947/1‑251 (MQ=255)                             gAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGAc                                                                                                                                                                                                                      >  1:284112/1‑250 (MQ=255)                                    gctcagctcaCCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAAGCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGt                                                                                                                                                                                                              <  2:445002/251‑1 (MQ=255)                                                  gCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCAGTGCTTACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCAGCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGGTGCCCGGCAGCCAGACGCAGGGCGCACTAAAACTTGCGATGCCTGCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGTGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTg                                                                                                                                                                                                 <  2:569229/250‑1 (MQ=255)                                                              acacCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATc                                                                                                                                                                                    <  2:201136/251‑1 (MQ=255)                                                               cacCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCg                                                                                                                                                                                   <  2:149396/251‑1 (MQ=255)                                                                     gCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTcc                                                                                                                                                                              >  2:453751/1‑250 (MQ=255)                                                                      cTATCCAGCAAAAGGATTTGCCATTGCTCACCAATAAACAAGCGCTTGACCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGACATCGTGGTTGCCCGGCAGCCAGAAGAAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTcc                                                                                                                                                                              >  1:332602/1‑249 (MQ=255)                                                                              cAACAGGAGTTGCCATTGCGCGCCAATAAGCACGCGCTTCGCCGGGGAGATACCCGCATCCTGGAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGCCGCAGGGCGCACGAAAACTTGCGGTGCCTTCAGCGAAATGCTGAGAGGCCGCAGAGAATTGATCCTGCGCTAGAGCACCTGTCGCGACAATCAGGTCGAATTCGGGCTGGTGCGGCCGAAGCGCCTCCAGGACCGCCt                                                                                                                                                                    <  2:88468/251‑1 (MQ=255)                                                                               aaCAGGATTTGCCATTGCTCACCAATAAACCCGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTg                                                                                                                                                                   <  2:417110/251‑1 (MQ=255)                                                                                         gCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTcc                                                                                                                                                         >  1:478151/1‑251 (MQ=255)                                                                                          ccATTGCTCAACAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAAGCGCCTCCAGCACCGCCTGGTAACTCTccc                                                                                                                                                        <  2:547629/251‑1 (MQ=255)                                                                                           cATTGCTCACCAAGAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCGCCAGCACCGCCTGGTAACTCTCCCa                                                                                                                                                       <  1:53672/251‑1 (MQ=255)                                                                                               gCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGACATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGCGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGAGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGgtg                                                                                                                                                   >  1:509426/1‑251 (MQ=255)                                                                                                    ccAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTAc                                                                                                                                              <  1:469546/251‑1 (MQ=255)                                                                                                     cAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGACTTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTAcc                                                                                                                                             >  2:48154/1‑251 (MQ=255)                                                                                                     cAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTAcc                                                                                                                                             <  2:451016/251‑1 (MQ=255)                                                                                                      aataagaacGCGCTTCGCCGGGGATATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGCAATCGTGGTTGCCAGGCAGCCAGACGCAGGGAGCACGAACACTTGCGATGCCTTCAGCGACATGCTGGTAGGCCGCCGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGCCGACTTCGTGCTGGTGTGGCCGAATCGCATCCAGCCCCGGCTGGTAACTCTCCCATGAGTTTAcac                                                                                                                                            >  2:548897/8‑249 (MQ=255)                                                                                                         aaccacGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCTAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTa                                                                                                                                         <  1:234359/248‑1 (MQ=255)                                                                                                         aaaCACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTa                                                                                                                                         <  2:409997/251‑1 (MQ=255)                                                                                                            cacGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAAc                                                                                                                                       >  1:488518/1‑250 (MQ=255)                                                                                                             ccgcgCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACAGCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTTGCCGAATCGCCTCCAGCACCGCCTGGGAACTCTCCCAGGTGTTTACCCCTAACAg                                                                                                                                     <  2:196283/250‑1 (MQ=255)                                                                                                             acGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAAc                                                                                                                                       >  1:381337/1‑249 (MQ=255)                                                                                                             acGCGCGTCGCCGGGGAGATACCCGCAGCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGAGCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACAGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAg                                                                                                                                     <  1:500583/251‑1 (MQ=255)                                                                                                                cgcTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCCCCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGc                                                                                                                                  >  2:424565/1‑251 (MQ=255)                                                                                                                    tCGCCGGAGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCCGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTGTACCCCTAACAGGGCTTc                                                                                                                               <  1:441614/250‑1 (MQ=255)                                                                                                                     cGCCGGGGAGATACCCGAATCCTGTAACGCGCTGGACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAACCTTGCGATGCCTTCAGCGAAATGCTTAGAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATTAGGTCGAATTCTGGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCGCCCAGGTGGTTACCCCTAACAGGGCTTCGt                                                                                                                             <  2:259003/251‑1 (MQ=255)                                                                                                                             agaTACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTtgt                                                                                                                     <  1:453751/251‑1 (MQ=255)                                                                                                                                 aCCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGCAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCa                                                                                                                  >  1:380314/1‑250 (MQ=255)                                                                                                                                  cccGCATCCTGTAACGCGATGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCaa                                                                                                                 >  1:105369/1‑250 (MQ=255)                                                                                                                                     gCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAAc                                                                                                               >  2:128018/1‑249 (MQ=255)                                                                                                                                      cATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCGAACAgg                                                                                                            >  1:369264/1‑251 (MQ=255)                                                                                                                                       aTCCTGTACCGCGCGGGACATCGCGGGCTGGAAAGCGTGGGTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACGTGCGATGCCTTCAGCGAAATGCTGATAGGCGGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGGCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGTTAACTCTCCCAGGTGTTTTCCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGt                                                                                                           <  2:511947/251‑1 (MQ=255)                                                                                                                                         ccTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCCGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAACTGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATAGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACACCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTgg                                                                                                         >  1:362911/1‑250 (MQ=255)                                                                                                                                         ccTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTg                                                                                                          >  1:390949/1‑250 (MQ=255)                                                                                                                                         ccTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGa                                                                                                         <  1:1942/251‑1 (MQ=255)                                                                                                                                            gTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGGAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAgtg                                                                                                      <  2:251408/251‑1 (MQ=255)                                                                                                                                            gTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAgt                                                                                                       <  1:521332/250‑1 (MQ=255)                                                                                                                                              aaCGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCTCGAAACCTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTc                                                                                                    <  2:247729/251‑1 (MQ=255)                                                                                                                                                 gcgcTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGt                                                                                                 <  1:374150/251‑1 (MQ=255)                                                                                                                                                   gcTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTaa                                                                                               <  1:498886/251‑1 (MQ=255)                                                                                                                                                       tACATCGCGAGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGAGCGCACGAAAACTTGCGATGCCTTCAGCGAACAGCTGATACGCCGCAGAGCATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCCGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGAGAGTGTCGGTAATTTg                                                                                           >  1:268544/1‑251 (MQ=255)                                                                                                                                                        aCATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGt                                                                                          <  1:57228/251‑1 (MQ=255)                                                                                                                                                            cgcgGGCTGGAAATCGTGGTTGCCCGGCAGCCAGCCGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCAAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTaaa                                                                                       >  2:93979/1‑250 (MQ=255)                                                                                                                                                                ggCTGGAAATCGGGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCt                                                                                  <  2:284112/251‑1 (MQ=255)                                                                                                                                                                    ggAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGAc                                                                               <  2:68896/250‑1 (MQ=255)                                                                                                                                                                          cGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACGTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAACCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCGTGTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGc                                                                         <  2:380314/250‑1 (MQ=255)                                                                                                                                                                           gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTTATCCTGCGCTAAATCACCAGTCGCGACAATCAGGTCGAATTCGGGCTGGTGTGTCCGAAGCGCCTCCAGCACCGCCTGGTAACTCTCCCTGGTGTTTACCCCTAACATGTCTTCGTGCGTTTGTGCAAAGAGGTGAGTTTCGGTAAGTTTTAAAGTCCTGACTCTGGCCt                                                                       <  2:227829/251‑1 (MQ=255)                                                                                                                                                                           gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCt                                                                       >  1:513422/1‑251 (MQ=255)                                                                                                                                                                           gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCt                                                                       <  1:563460/251‑1 (MQ=255)                                                                                                                                                                           gTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGACTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCt                                                                       >  1:356627/1‑251 (MQ=255)                                                                                                                                                                                 gCCCGGCAGCCAGACGCAGGGCGCACGAAAACTGGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCAccag                                                                 <  2:486530/251‑1 (MQ=255)                                                                                                                                                                                        aGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAgg                                                          <  2:182697/251‑1 (MQ=255)                                                                                                                                                                                        aGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAgg                                                          <  1:325489/251‑1 (MQ=255)                                                                                                                                                                                               cGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGtt                                                   >  1:444104/1‑251 (MQ=255)                                                                                                                                                                                                              aaCTTGCGATGCTGTCAGCGAAATGCTGATAGGCCGCAGAGAATTGATCCGGCGCTAAATCACCTGGCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCTAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCaaa                                    <  2:335177/251‑1 (MQ=255)                                                                                                                                                                                                                                    aTGCTGATAGGCCGCAGAGAATTGATCCTGCGGTAAATCCACTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGTCCGAATCGCCTCCAGCACCTGCTGGTAACTCTTCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTTTGCAAACAGGTTAGTGTCGGTAATTTTGAAAATCCTGACTCTGGCCGCACCATCCAGAGTAATGTTTAACATTCTTTCTAAAGGGTTTCCTTAGGTTTCacgac              <  2:295992/251‑1 (MQ=255)                                                                                                                                                                                                                                     tGCTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGCCAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCacgacg             >  1:510925/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCt           >  2:38461/1‑251 (MQ=255)                                                                                                                                                                                                                                       cTGATAGGCCGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCAAAATGGTGTCCTTAGGTTTCACGACGCt           >  1:446268/1‑251 (MQ=255)                                                                                                                                                                                                                                               ccGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAAcc   >  1:346346/1‑251 (MQ=255)                                                                                                                                                                                                                                                cGCAGAGAATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACAACGCTAATAAACCg  <  1:211876/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GCCAGTTTACGTTCCAGCCACTCAAGCTGAAACTCGCTCAGCTCACCGTGCGGCACGCCAAACACCTGGCTATCCAGCAACAGGATTTGCCATTGCTCACCAATAAACACGCGCTTCGCCGGGGAGATACCCGCATCCTGTAACGCGCTGTACATCGCGGGCTGGAAATCGTGGTTGCCCGGCAGCCAGACGCAGGGCGCACGAAAACTTGCGATGCCTTCAGCGAAATGCTGATAGGCCGCAGAGGATTGATCCTGCGCTAAATCACCTGTCGCGACAATCAGGTCGAATTCGTGCTGGTGTGGCCGAATCGCCTCCAGCACCGCCTGGTAACTCTCCCAGGTGTTTACCCCTAACAGGGCTTCGTGCTTTTGTGCAAACAGGTGAGTGTCGGTAATTTGTAAAATCCTGACTCTGGCCTCACCAGCCAGAGGAAGGGTTAACAGGCTTTCCAAATGGTGTCCTTAGGTTTCACGACGCTAATAAACCG  >  NC_000913/3176379‑3176868

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