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	97.9%	148.7 / NA	47	G220D (GGC→GAC)	pitA	metal phosphate:H(+) symporter PitA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/24);  minor base C (1/0);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

GACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAACC  >  NC_000913/3638061‑3638543                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    gACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGAGCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGtttt                                                                                                                                                                                                                                          <  2:7865/251‑1 (MQ=255)    gTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGc                                                                                                                                                                                                                                       <  1:195543/251‑1 (MQ=255)        gtggtgGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACgg                                                                                                                                                                                                                                   <  1:198192/251‑1 (MQ=255)          ggtggATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGGATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCCCGgcg                                                                                                                                                                                                                                 >  1:266918/1‑251 (MQ=255)                   aCTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCTAAAAACGCGCCCGTATCCACCTGCCCCCAGCGGAGCGGAAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGAGCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATg                                                                                                                                                                                                                        <  1:261079/251‑1 (MQ=255)                       aaTATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTTCTGCGTCGCTACTGGAGCGGCCCCAAGAAAAGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAAGCGTCTTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTcc                                                                                                                                                                                                                    >  2:121314/1‑250 (MQ=255)                          atCCCGAAAGTATTAGTTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCGGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGa                                                                                                                                                                                                                  <  2:253431/250‑1 (MQ=255)                            cccGAAAGTATTAATTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAg                                                                                                                                                                                                               <  1:79154/251‑1 (MQ=255)                                        aagtattTGCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGTGCGTGAAAAGAAAGACGGCAAGAAAAATCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTGCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTg                                                                                                                                                                                                    <  1:112281/250‑1 (MQ=255)                                         tgtattTTCGGTTCTTTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACGGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAg                                                                                                                                                                                                               <  2:7711/237‑1 (MQ=255)                                         agtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGAAGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAg                                                                                                                                                                                                               >  1:7711/1‑238 (MQ=255)                                          gtattTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGtt                                                                                                                                                                                                 <  1:249259/251‑1 (MQ=255)                                                             cGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGCGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTgggg                                                                                                                                                                              >  2:153894/1‑249 (MQ=255)                                                                  ccccTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGGCTGGTTATGTTGGTATTGATTGGCGTcgcg                                                                                                                                                                         <  1:247422/251‑1 (MQ=255)                                                                          gTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGGGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAgg                                                                                                                                                                 >  1:277214/1‑251 (MQ=255)                                                                             ggCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGCACGCTGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCtt                                                                                                                                                              >  2:28666/1‑251 (MQ=255)                                                                                cTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGCGCGTGAAAAGAAAGACGGCAAGAAAAACCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGAGTCGCGCCCGCAGGCTTCgt                                                                                                                                                           >  1:49411/1‑251 (MQ=255)                                                                                 tGGTGTTTGCTGGCGGTCTGCTTTTCTTTCTTCGTCGCTACTGGGGCGGCACCAAGAAACGCCCCCGTATCCACCTGACCCCCGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCCGAACGGCATTGGTCTTGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCgtg                                                                                                                                                          >  2:84163/1‑251 (MQ=255)                                                                                       ttGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGCGCGTGAAAAGCAAGACGGCCAGAAAAAGCCGCCGTTCTGGACGCGTCTTGCGCTTATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGGTATGTTGGTATTGATTGGCGTCGCGCCAGCCGGCTTCGTGGTGAAc                                                                                                                                                    >  1:273819/1‑251 (MQ=255)                                                                                              cGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATg                                                                                                                                             <  2:277214/251‑1 (MQ=255)                                                                                               ggTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGc                                                                                                                                            <  1:102118/251‑1 (MQ=255)                                                                                                             gctgcGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCCCCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGaaa                                                                                                                              <  1:153894/251‑1 (MQ=255)                                                                                                                  gTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacc                                                                                                                         >  1:32592/1‑251 (MQ=255)                                                                                                                        ccgggggCGGCACCAATAACCGCGCCCGTATCCACCTTACCCCAGCGGAGCGTGAAAAGATAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCacccgtacc                                                                                                                   <  2:173989/245‑1 (MQ=255)                                                                                                                        aCTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTa                                                                                                                                   >  1:248464/1‑235 (MQ=255)                                                                                                                        aCTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTa                                                                                                                                   <  2:248464/235‑1 (MQ=255)                                                                                                                                ggCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGc                                                                                                           <  2:101239/251‑1 (MQ=255)                                                                                                                                    ccAAGAAACGCGCCCGTACCCACCTGCCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATGGCGCTGATCCTTTCCGCTAGCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAAGGCCACTGGCTACGAAATCACCCGTACCCGTGTTGCCATc                                                                                                       <  2:273819/251‑1 (MQ=255)                                                                                                                                            cgcgCCAGTATCCAACTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGACCCTTTCCGCTATCGACGGGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGGATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACAGGAATGCAACTGGCAACGAAATCGCACGTACACGTGATGCCATCAACAACGt                                                                                               >  1:64555/1‑251 (MQ=255)                                                                                                                                                      tCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGTAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTAc                                                                                     >  2:136727/1‑251 (MQ=255)                                                                                                                                                              aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGTCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCTCCCTTACCCGTTATGCCATCAACAACGTCGAAGCTTACTTTGagca                                                                             <  1:254023/251‑1 (MQ=255)                                                                                                                                                              aCCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGagca                                                                             <  1:136727/251‑1 (MQ=255)                                                                                                                                                                               aaagaaagACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTc                                                             >  1:187831/1‑250 (MQ=255)                                                                                                                                                                                aagaaagaCGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCaa                                                           <  1:79192/251‑1 (MQ=255)                                                                                                                                                                                           cAAGAAACAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCg                                                >  2:153033/1‑251 (MQ=255)                                                                                                                                                                                            aaGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCgg                                               <  1:256008/251‑1 (MQ=255)                                                                                                                                                                                                    gccgccGTTCTGGACGCGTATTGCGCTGATCAATTCCGCTGTAGCCGTAGCGTTTTCGCACGGCGCGAACGATGGTCAGAAACGCATCGGCCTGGTTATGTTGGTATTGATTGGCGTCGCGCCCGCAGACTTCGAGGTGAACATGAATGCCACTGACTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATc                                       >  1:168194/1‑251 (MQ=255)                                                                                                                                                                                                          gTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACa                                                        >  2:93077/1‑228 (MQ=255)                                                                                                                                                                                                          gTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACa                                                        <  1:93077/228‑1 (MQ=255)                                                                                                                                                                                                            tcgtggcgcgTATTGCTCTTATCCTTTCCGCTATCGACGTGGCGGTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTTAGTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACTTTGAAGCTTACGTTGAGCAGCATCCTGCTCTGCTCAAACAGGCTACCGGTGCTGATCATTTAGt                                <  2:134668/244‑1 (MQ=255)                                                                                                                                                                                                                 aCGCGTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGc                          >  1:118738/1‑251 (MQ=255)                                                                                                                                                                                                                   gcgTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTc                        >  1:209897/1‑251 (MQ=255)                                                                                                                                                                                                                     gTATTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGCTGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCAACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCg                      >  1:201340/1‑251 (MQ=255)                                                                                                                                                                                                                       aTTGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCgg                     <  1:28666/250‑1 (MQ=255)                                                                                                                                                                                                                         tGCGCTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGCATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCCGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCCGTTAGTACCGGCCCCGGa                    >  2:129533/1‑249 (MQ=255)                                                                                                                                                                                                                             cTGATCCTTTCCGCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATAAGTTAGTACCGGCTCAGGAAGCTgg              >  1:251025/1‑251 (MQ=255)                                                                                                                                                                                                                                         gCTATCGACGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAAcc  >  1:58279/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    GACGTCAGTGGTGGATGCACTCAATATCCCGAAAGTATTAAGTATTTTCGGTTCTCTGATCGTTTCCCCTATTGTCGGCCTGGTGTTTGCTGGCGGTCTGATTTTCTTGCTGCGTCGCTACTGGAGCGGCACCAAGAAACGCGCCCGTATCCACCTGACCCCAGCGGAGCGTGAAAAGAAAGACGGCAAGAAAAAGCCGCCGTTCTGGACGCGTATTGCGCTGATCCTTTCCGCTATCGGCGTGGCGTTTTCGCACGGCGCGAACGATGGTCAGAAAGGCATTGGTCTGGTTATGTTGGTATTGATTGGCGTCGCGCCAGCAGGCTTCGTGGTGAACATGAATGCCACTGGCTACGAAATCACCCGTACCCGTGATGCCATCAACAACGTCGAAGCTTACTTTGAGCAGCATCCTGCGCTGCTCAAACAGGCTACCGGTGCTGATCAGTTAGTACCGGCTCCGGAAGCTGGCGCAACGCAACC  >  NC_000913/3638061‑3638543

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