Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,638,300	G→A	G220D (GGC→GAC)	pitA →	metal phosphate:H(+) symporter PitA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,638,300	0	G	A	96.8%	186.2 / NA	64	G220D (GGC→GAC)	pitA	metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (31/30);  minor base C (2/0);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAACC  >  NC_000913/3638058‑3638543                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    cGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCCGCGGAGCGTGAAAAAAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGCTCCATTCCGCTATCgacgtggcgg                                                                                                                                                                                                                                             >  2:20072/1‑250 (MQ=255)           gtggtgGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACgg                                                                                                                                                                                                                                   >  1:288843/1‑251 (MQ=255)              gtggATGCCCTCAATATCCCGACAGAATTAAGTATTTTCGGTTCTCTGATCGATTCCCATATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTCGGTCCCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGATCGTGAAAAGAAAGCCGGCAAGCAAAAGCCGCCGTTCCGGCCGCGGATTGCGATGATACTTTCCGCTATTGACGTGGCGTTTTCGCACGgcgc                                                                                                                                                                                                                                >  2:33618/1‑251 (MQ=255)                         cAATAACCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGCGTTTGCTGTCGGTCTGATTTTCTTGCTGCTTCGCGCCTGGAGCGGAACAAAGAAACGCGCCCGTATCCACCTGAGCCCAGCTGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCGATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTc                                                                                                                                                                                                                     <  1:370665/251‑1 (MQ=255)                          aaTATCCCGAAAGTATTAAGTATTGTCGGTTCTCTGACCGTTTCCCCTTTTGCCGGCCTGGTGTTTGCTGGTGGTCTGATTTTCTTGCTGCGTCGCTCGTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCTTTTTCGCACGGCGCGAACGGTgggca                                                                                                                                                                                                                    <  2:162488/251‑4 (MQ=255)                             atCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGGTGGTccga                                                                                                                                                                                                                  >  2:111114/1‑247 (MQ=255)                             accccGAAAGTTTTAAGTATTTTCGGTTCTCTGTGCGTTTCCCCTATTGTCGGCGTGTTGTTTGCTGGCGTTCTGATTTTCTTGCTGCGTCGCTCCGGGAGCGGCACCAAGAAACGCGCCCGTATCCACCGGACCTCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGa                                                                                                                                                                                                                  <  1:282226/248‑1 (MQ=255)                              tCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCCAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGCCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGaa                                                                                                                                                                                                                 >  1:337100/1‑250 (MQ=255)                               cccGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGaa                                                                                                                                                                                                                 >  1:313268/1‑249 (MQ=255)                               cccGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAg                                                                                                                                                                                                               <  2:107319/251‑1 (MQ=255)                                          taagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGGCGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTg                                                                                                                                                                                                    >  2:412057/1‑251 (MQ=255)                                           aagtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTTCTGGAGCGGCCCCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGCAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGACAGGCATTGGTCTgg                                                                                                                                                                                                   >  2:384739/1‑251 (MQ=255)                                             gtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCTGGGTGTTTGCTGGCGGTCTGATTTTCGTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCAGTGGTCTGGtt                                                                                                                                                                                                 <  2:362866/251‑1 (MQ=255)                                                   tCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGGCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTg                                                                                                                                                                                           >  1:221770/1‑251 (MQ=255)                                                        tctcTGATCGTTTCCCCTATTGTCCGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTTGCTACCGGGGCGGCAACAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGCGCGGGAAAAGAAAGACGGCCAGAAAAAGCCGCCGTTCTGGACGCGTAATGCGCTGATACTTTCCGCTATCGGCGTGGCGCTTTCGCACCGCGCGAAAGGTGGTCTGAAAGGCCTTGGTCTGGTTATGTTGGTatt                                                                                                                                                                                      >  2:252390/1‑251 (MQ=255)                                                           ctGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTATGATTTTCTTGCCGAGTCACTACTGGAGCGGAACCAAGAACCGCGCCCGTATCCACCTGAACCCAGCGGAGCGTGGAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGAGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGGTATGTTGGTcttgat                                                                                                                                                                                   >  2:283067/1‑251 (MQ=255)                                                              ttCGTTTCCCCGTTTGTCGGCGTGGTGTTTGTTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTattgattg                                                                                                                                                                                 <  2:233324/249‑1 (MQ=255)                                                               tCGTTTCCCCTATGGGCGGCCTGGGGTTTGCTGGCGGTCTGATTTTCTTGCTGGGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCACAGCGGAGCGTGAAAAGAAAGACGGCACGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTGATGTTGGTATTGATTgg                                                                                                                                                                                <  2:131459/250‑1 (MQ=255)                                                                 gTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGt                                                                                                                                                                             >  1:395884/1‑251 (MQ=255)                                                                     caccTATTGTCTGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTTCGTCGCTACTTGAGCGGAACCAAGAAACGCGCCCGGATCCACCTGACCCCAGGGGAGCGTGCAAAGAAAGACGGCCAGAAAAAGCCGCCGTTCTGGGCGACAATTGCGCTGATCCTTTCCGCTATCGACGGGGCGTTTTCGCACGGGCCGCACGATAGTCAAAAAAGCAATGGTCTGGGTGTGTTGGTATTGATTGGCGTcgcg                                                                                                                                                                         >  2:328024/3‑251 (MQ=255)                                                                       ccTATTGTCGGACTGGTGTTTGCTGGCGGTCTGAATTTCTTGCTGCCTCGCTACTGAAGCGGCACCAAGAAACGCGCCCGTATACACCTGACCCCAGCGGCGCGTGAAAAGAAAGACGGCAAGAAAAAACCGCCGTTCTGGACGCGTATTGCCCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGCCCCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTccagcc                                                                                                                                                                       >  2:380550/1‑246 (MQ=255)                                                                             gTCGGCCTGGTTTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCcagcag                                                                                                                                                                  <  1:52098/250‑1 (MQ=255)                                                                               cGGCCGGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGGGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGc                                                                                                                                                                <  2:384060/250‑1 (MQ=255)                                                                                 gcctgttttttgctgtctgacTGATTTTTTTGCTGCGTCGCGACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCGGCGGAGCGTTAAAAGAAAGACGGCAAGAAAAACCCGCCGTTCTGGACGCGTATTGCGCTTATCCTTTCCGCGATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTc                                                                                                                                                             <  1:427858/231‑1 (MQ=255)                                                                                                          tttCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGAGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCt                                                                                                                                    >  1:387300/1‑251 (MQ=255)                                                                                                                gctgcGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCGGTCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGaaa                                                                                                                              <  2:395884/251‑1 (MQ=255)                                                                                                                      tCGCTACTGGAGCGGCACCAAGAAACGCGCCCGAATCCACCTGACACCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGATCTGGACGCGTATTGCGCTGATCCTTTCAGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGATCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacac                                                                                                                        >  1:371735/1‑249 (MQ=255)                                                                                                                       cGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAATGCATTGGTCTGGTTATGTTGTTATTGATTGGCGTCTCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCaccc                                                                                                                        <  1:267122/250‑1 (MQ=255)                                                                                                                        gCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTAGCGACGTGGCGTTTTCGCACGGCGCGAACGATGGGCATAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGGGAACATGAATGCCACTGGCTACGAAATCacccgt                                                                                                                      <  1:38276/251‑1 (MQ=255)                                                                                                                           aCTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgtacc                                                                                                                   >  2:105663/1‑251 (MQ=255)                                                                                                                             gggAGCGGCACAAAGAAACTCGTCCGTATCCCACTGACCTCAGGGGTGGGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGAGTATTGCGCTTATCCTTTCCGCGATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGGATTGATTGGCGGCGCGCCAGCAGGTTTCGTGGTGAACATGAATGCCACTGTCTACGAAATCacccgtacccg                                                                                                                 <  2:250114/250‑1 (MQ=255)                                                                                                                               gAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTg                                                                                                               >  1:389120/1‑251 (MQ=255)                                                                                                                                  cGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGAGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAAAATGAATGCCACAGGCTACGAAATCACCCGTACCCGTgctg                                                                                                            >  2:7437/1‑248 (MQ=255)                                                                                                                                  cGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGAGAAAAGAAAGCCGGCAACAACAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGCCGTGGCGTTTTCGCACGGCGCGAATGATGGTCAGAAAGGCATTGGTCTGGTTCTGTTGGTATTGCTTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCccccgtacccgttgtg                                                                                                            >  2:328285/1‑247 (MQ=255)                                                                                                                                  cGGCACAAAGAAACCCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGCCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTAGCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGTTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATg                                                                                                            <  2:387300/251‑1 (MQ=255)                                                                                                                                         aaGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATca                                                                                                      >  1:24470/1‑250 (MQ=255)                                                                                                                                                gcgcCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTc                                                                                              >  1:177839/1‑251 (MQ=255)                                                                                                                                                 cgcCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCg                                                                                             <  1:13131/251‑1 (MQ=255)                                                                                                                                                      gTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCt                                                                                        <  1:33618/251‑1 (MQ=255)                                                                                                                                                      gTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCt                                                                                        <  1:289993/251‑1 (MQ=255)                                                                                                                                                                 gcccacgCGGGGCGTGAACAGGAAGACGGCAAGAAAAAGCCGCCGTGCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGGTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGagca                                                                             <  2:96639/245‑1 (MQ=255)                                                                                                                                                                 aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGagca                                                                             <  1:195479/251‑1 (MQ=255)                                                                                                                                                                                    tgaaagaCGGAAAGGAACAGCCGCCGGTCTGGACGCGTATAGCGATCATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGCTGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGGGAACATGAATGCCTCTGGCTACGAAATCACCCGTTCCCGTGACGCCATCAACAACGTCGAAGCTTACTTTGAGCCGCATCCTGCGCTGATCaa                                                           >  2:135312/2‑250 (MQ=255)                                                                                                                                                                                      aaagaCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGGACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCGGCATCCTGCGCTGCTCAAAc                                                         >  2:14350/1‑250 (MQ=255)                                                                                                                                                                                       aagaCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAAc                                                         >  2:387761/1‑249 (MQ=255)                                                                                                                                                                                        agaCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCAGTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCGCAAACAgg                                                      <  2:337100/251‑1 (MQ=255)                                                                                                                                                                                         gaCGGCAAGAAAAAGCCGCCGTTCTGGACGAGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTAGCACGGCGCGACCGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGACTGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGACGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGc                                                     >  1:129773/1‑251 (MQ=255)                                                                                                                                                                                             gCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTAcc                                                 <  1:405222/251‑1 (MQ=255)                                                                                                                                                                                                aGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGt                                              >  1:15604/1‑251 (MQ=255)                                                                                                                                                                                                 aaaaaaGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTg                                             <  1:385545/250‑1 (MQ=255)                                                                                                                                                                                                    aaaGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGATGCTg                                          >  1:3616/1‑251 (MQ=255)                                                                                                                                                                                                         cgccgTTCTGGACGCGTATTGCGCTGATCCGTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAg                                     <  2:217609/251‑1 (MQ=255)                                                                                                                                                                                                            cgTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGtt                                   >  1:200612/1‑250 (MQ=255)                                                                                                                                                                                                                cTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTAc                              >  1:102504/1‑251 (MQ=255)                                                                                                                                                                                                                  ggACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCg                            <  1:20072/251‑1 (MQ=255)                                                                                                                                                                                                                   gACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCgg                           <  1:22821/251‑1 (MQ=255)                                                                                                                                                                                                                       cgTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATAAGTTAGTACCGGCtac                       >  1:139008/1‑249 (MQ=255)                                                                                                                                                                                                                              cgcTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCt                <  1:384739/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGgcgc           >  1:167404/1‑251 (MQ=255)                                                                                                                                                                                                                                      cTTTCAGCTATCGACGTGGCGTTTTCGCACGGCACGAACGAAGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCGACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGCTCAATTAGTACCGGCTCCGGAAGCTGGcgcaa         >  1:362105/1‑250 (MQ=255)                                                                                                                                                                                                                                       tttCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGcgcaacg       >  2:349824/1‑251 (MQ=255)                                                                                                                                                                                                                                        ttCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGcgcaac        <  1:111114/249‑1 (MQ=255)                                                                                                                                                                                                                                          ccGCTATCGACGTGGCGTTTTCGTACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGTTATGTTGGTATTGATTGGCGGCGCGCCAGCAGGCTTCGTTGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGTCTCCGGAAGCTGGcgcaacgca     <  2:225243/250‑1 (MQ=255)                                                                                                                                                                                                                                           cGCTATCAACGTGGCGTTTGCGCACGGCGAGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATCGATTGGCGTCGCGACAGCAGGCTTCGTGGCGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGATTACTTTAAGCAGCATCCTGCGCTGCTCAAACAGACTCCCGGTGCTGATCAGCTAGTACCGGCTCCGGAAGCTGGCgcaacgcaac   >  2:374747/1‑251 (MQ=255)                                                                                                                                                                                                                                            gCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAAcc  <  2:313268/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    CGGGACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAACC  >  NC_000913/3638058‑3638543

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