Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	1,562,470	T→C	19.5%	G58G (GGA→GGG)	Y75_RS07840 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	1,562,470	0	T	C	19.5%	87.9 / 4.9	53	G58G (GGA→GGG)	Y75_RS07840	ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (29/7);  new base C (2/8);  total (38/15)
Fisher's exact test for biased strand distribution p-value = 7.73e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAG  >  NC_007779/1562236‑1562702                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           tattCGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGCAGCGGTAAAACTGTACATACAGGGGCTGGG‑CCCCTccc                                                                                                                                                                                                                                         >  1:936574/1‑238 (MQ=255)     cGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAC‑CAGGTCGCTGCCGTAACGGTAAAACTGTACATACAGTGGCTGGG‑CCAGTCCCAAttt                                                                                                                                                                                                                                    >  1:1994569/1‑238 (MQ=255)     cGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑GCCGGCCCCGtt                                                                                                                                                                                                                                     >  2:233507/1‑238 (MQ=255)     cGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCCGTCCCCGtt                                                                                                                                                                                                                                     >  1:602555/1‑238 (MQ=255)       cACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGCTCCAggt                                                                                                                                                                                                                                     >  1:335535/1‑234 (MQ=255)           aGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑TCCGTCCCatgttccggcg                                                                                                                                                                                                                              >  1:2161096/1‑229 (MQ=255)                 tCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATAtt                                                                                                                                                                                                                        <  1:79396/239‑1 (MQ=255)                 tCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATAtt                                                                                                                                                                                                                        <  1:52059/239‑1 (MQ=255)                 tCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATAtt                                                                                                                                                                                                                        <  1:437291/239‑1 (MQ=255)                 tCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATAtt                                                                                                                                                                                                                        <  1:1741295/239‑1 (MQ=255)                  ccagccagCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGCGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCGGGCCCAACTGCTgggggagagtt                                                                                                                                                                                                                       >  2:932106/1‑230 (MQ=255)                      ccagcGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAATCCCAAGTGGTGGGGGACAtttttc                                                                                                                                                                                                                    >  2:461657/1‑236 (MQ=255)                       cagcGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCAc                                                                                                                                                                                                                  <  1:1588050/239‑1 (MQ=255)                       cagcGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGCACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGGACCAGGTGCTGGCGAATATTTTTCAc                                                                                                                                                                                                                  >  1:316547/1‑239 (MQ=255)                                       cTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGGCCCAGTTTCTGGCGAATATTTTCCACGATAGAATCGcctgc                                                                                                                                                                                                   >  1:1200028/1‑234 (MQ=255)                                               agTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑ACAGTCCCCGGTGCTGGCGGAAATTTTCCACGATAGCATCGCTTGCGCGGGCAc                                                                                                                                                                                           >  1:583146/1‑238 (MQ=255)                                                gTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTTCACGATAGCATCGCTGGCGCGGTCAcc                                                                                                                                                                                          >  1:425702/1‑238 (MQ=255)                                                 tATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGCTCCAGGTGCGGGGGCATATTTTACACGAAAGCATCGCGGGCGCGGTCAccc                                                                                                                                                                                         >  1:2113401/1‑238 (MQ=255)                                                  aTGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTTGCGCGGTCCCCCg                                                                                                                                                                                        >  1:674640/1‑238 (MQ=255)                                                    gCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCCGTTGCTGGCGAATCATTTCCACGGTAGCCACGCTGGCGCGGTCACCCGccc                                                                                                                                                                                     >  1:98259/1‑238 (MQ=255)                                                    gCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGGGAAAATTTTCCACGATAGGATCGCCGGCGCGGGCACCCGCCa                                                                                                                                                                                     >  1:175598/1‑239 (MQ=255)                                                    gCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCCATCCCAATTTCTGGCGAAAATTTTTCAAGATAGCATCGCTGGCGCAGTCACCCGCCa                                                                                                                                                                                     >  1:2228348/1‑239 (MQ=255)                                                    gCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCCTGCGCGGTCAccccccc                                                                                                                                                                                     >  1:1021/1‑238 (MQ=255)                                                    gCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGTCCCAGGTGCTCGGGAATATTTTCCACGAAAGCATCGCTTGGGCGGGCACCCGccg                                                                                                                                                                                     >  1:2082243/1‑238 (MQ=255)                                                         gCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACGGGACGCCCGGAACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGACGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGACCCAGTTGCCGGCGAATATTTTCCACGAAAGCATCCCTGACGCCGGCACCCCCCAGGaa                                                                                                                                                                                 >  1:1811464/1‑238 (MQ=255)                                                         gCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCCGTTCCAGTTGCTGGGGAATATTTTCCACGAAAGAATCGCTGGCGCGGTCCCCCGCCAGtga                                                                                                                                                                                 >  1:117993/1‑236 (MQ=255)                                                         gCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCCGTCCCCGGTGCTGGCGAAAATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGcccggaa                                                                                                                                                                                 >  1:147681/1‑233 (MQ=255)                                                           gggATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGGCCCAGGTGCCGGCGAACATTTTTCACAAAAGCACCCCTGGCGCCGGCCCCCCCCCGTaaac                                                                                                                                                                               >  2:1946506/1‑236 (MQ=255)                                                                            tGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAata                                                                                                                                                                                                                          <  2:40274/178‑1 (MQ=255)                                                                            tGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAata                                                                                                                                                                                                                          >  1:40274/1‑178 (MQ=255)                                                                                  cagcagGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑GCAGTCCCCGTTGCTGGCGAAAATTTTCCACGATAGCCTCGCTGGCGCGGGCACCCGCCCGGAAACGCGCCCGATCACCAGGGATCCGaa                                                                                                                                                       >  1:1164308/1‑238 (MQ=255)                                                                                       gggCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGACCCAGTTGCTGGCGAAAATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCCGGAACCGCCCCGGGTCACCAGGGATCCGAGGCgag                                                                                                                                                  >  1:863638/1‑239 (MQ=255)                                                                                        ggCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGCATATTTTCCACGATCGCATCGCCGGCGCGGGCACCCGCCAGTAAACGCGCCCGATCACCAGGGGACCGATGCgagg                                                                                                                                                 >  1:861310/1‑238 (MQ=255)                                                                                                       cTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGTCCCCGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCCCCAGGAACCGCCCCGGGTCACCCGGGGTCAGATGCGCGATAATAAACGT‑aaac                                                                                                                                  >  1:1503431/1‑236 (MQ=255)                                                                                                                                   aaCCCTCCCACCCCCGGACGCCCGGTAGGAATGGATGTCCCCAGGTCACCATGAAA‑CAGGTGGCTGAGGTAGCGGTAAACCTGTACATACAGTGGCGGGC‑CCAGTCCCGGTTGCGGGAGAATATTTCCCACGATACCATCGCTGGCGGGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCgag                                                                                                                                                  <  1:561410/195‑1 (MQ=255)                                                                                                                                                  ggACGCCCGGTACGAATGGATGTTCCCAGGTCACCTTGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCATTCCCCGTTGCTGGCCGATATTTTTCACGCAAGCCTCGCTTGCGCGGGCACCCCCCAGTAAACGCGCCGGGTCACCCAGGAACACAAGCGCGATAATAAACGTAAACCACGC‑GACGCCCCCCACCACCCGCACACAAcccccgacaacttc                                                                                      >  2:549293/1‑230 (MQ=255)                                                                                                                                                                 atgtatttccccAGGCCCCCATGAAA‑CAGGTCGCTGAGGTGGCGGAAAACCTGTCCATACAGGGGCGGGC‑CCGGCCCCGGTTGCTGGCAAATATTTCCCACAATAGCATCGGGGGGGCGGCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCa                                                                      <  1:259954/230‑1 (MQ=255)                                                                                                                                                                         ccccAGGTC‑CCAAGAAACCGGGTCGCTCACGTGGGGATAAAACTGTACATACAGGGGCTGCC‑CCACCCCCAGTTGCGGGCAAATTTTTCCCAGAATAGCATCGCGGGCGCGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCa                                                              <  1:233507/238‑1 (MQ=255)                                                                                                                                                                           ccggaccACCATAAAA‑GTGGTCGCTGCGGTGGGGAAAAAACTGTACATACGGTGG‑TGGTCCCACCCCCGGTTGTGGGCAATTTTTTCCCACGATAGCCTCGCGGGCGCGGTCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACg                                                            <  1:21146/233‑1 (MQ=255)                                                                                                                                                                                               gTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAg                                                                                                                            >  1:2000443/1‑157 (MQ=255)                                                                                                                                                                                               gTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGG‑CCAGTCCCCGGTGCTGGCGGATATTTTCCACCATAGCATCGCTGGCGCCGTCCCCCCCCCGGGACCGCGCCCGATCACCAGGGATCCACTGCGGGATAAAAAAAGT‑AATACCCCCGACCCCCCCCACCACCAGCAA‑CAACCCCCCGAAACCGTGGCGCAAAACACCCCCGAACGACACACTTTTTACTCACggag                                            >  2:1620193/1‑234 (MQ=255)                                                                                                                                                                                                 cGCTGACGTAGCGGTAAAACTCTACATACGGTGGCGGCC‑CCAGCCCCAGTTGCGGGAAATTATTTTCCACGATCGCTCCGGGGGCGCGCCCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGtattga                                        <  1:461657/238‑1 (MQ=255)                                                                                                                                                                                                         tgggggaaaacCTGTACACACGGGGGCGGTC‑CCGGCCCCGGTTGGGGGCAAATTTTCCCCACGATAGCAGCGCGGGCGGGGCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACCCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAg                                <  1:1946506/231‑1 (MQ=255)                                                                                                                                                                                                              gTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAAACGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCACGGtattgat                                       >  1:682344/1‑226 (MQ=255)                                                                                                                                                                                                                     tttacataggGTGGGTGGC‑CCGGCCCCGGGGGCGGGCAAATATTGCCCAGTAGAGCGGCGGTGGCGCGACCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCa                   <  1:1208743/237‑1 (MQ=255)                                                                                                                                                                                                                       gacataggggggCTGCT‑CCGGCCCCAGTGGCGGGCGAATATTTCCCGCGATAACATCGGTGGCGGGGACACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATg                 <  1:1620193/229‑1 (MQ=255)                                                                                                                                                                                                                               gTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGa                               >  1:258681/1‑217 (MQ=255)                                                                                                                                                                                                                                tGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCCCCCGCCCGGAACCGCGCCGGATCACCACGGATCAGATGCGGGATAATAAACGT‑AATCACCCAGACGCCCGCCCCCCCGAGCAC‑CAGCCCCCACCAGCGTTGGCGTAAAAAACCCCACAAGGGCACGATATTTACTCATGGTATTGTTATTTAAGACCTGTTTCAGCATGGGATTGAac        >  2:2075121/1‑239 (MQ=255)                                                                                                                                                                                                                                tGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGGCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAac        >  1:1688194/1‑239 (MQ=255)                                                                                                                                                                                                                                 ggCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGaa                                 >  1:1062672/1‑213 (MQ=255)                                                                                                                                                                                                                                    tGGT‑CCAGCCCCTTTCGCGGCCAAATTTTTCCCAGGAACGCATCGCTGCCGGGCCCCCCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACaa    <  1:778351/239‑1 (MQ=255)                                                                                                                                                                                                                                      gT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CCGCCCCCAGCAGCGTTGGCGCAAAATACTCCAGAAGGACACggt                                                          >  1:389851/1‑181 (MQ=255)                                                                                                                                                                                                                                      gT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTCGGGCGGGCACCCGCCAGTAACCGCGCCGGAACACCAGGGCTCCGATGCGAGAAAATAAACGT‑AATCCCCCAGACCCCCGCCCCCACCAGCAA‑CAACCCCCACCAGGCCTTGCGTAACATACTCCCAAAGGGCACGGTATTTACTCATGGTATTGATATTTAAAACCCGTTCCCGCATGGGATTTAaccaaaaa  >  2:861528/1‑233 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           TATTCGCACCAGATGATCCAGCCAGCGATTTCGCCAGACTGCAGAGAGTATGCCCAGCGGGATGCCAATCAGGAGTGCCAGCAGCAGGGCACCAAAAGCCAGCTCCAGCGTCGCCGGGAAAAATATACGCAACTCTTCCAGCACCGGACGCCCGGTACGAATGGATGTTCCCAGGTCACCATGAAA‑CAGGTCGCTGACGTAGCGGTAAAACTGTACATACAGTGGCTGGT‑CCAGTCCCAGTTGCTGGCGAATATTTTCCACGATAGCATCGCTGGCGCGGTCACCCGCCAGTAACCGCGCCGGATCACCAGGGATCAGATGCGAGATAATAAACGT‑AATCACGCAGACGCCCGCCACCACGAGCAC‑CAGCCCCCAGCAGCGTTGGCGTAAAATACTCCAGAAGGTCACGTTATTTACTCATGGTATTGATATTGAAGACCTGTTCCAGCATGGGATTGAACACAAAG  >  NC_007779/1562236‑1562702

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