Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,370,266	0	T	G	25.0%	123.9 / 11.3	72	G83G (GGT→GGG)	arnD	putative 4‑deoxy‑4‑formamido‑L‑arabinose‑phosphoundecaprenol deformylase ArnD
Reads supporting (aligned to +/- strand):  ref base T (9/45);  new base G (16/2);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 5.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGCAAAAGAAGCGTTC  >  NC_000913/2370018‑2370512                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        atGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTATGGCGACTGGTGAAGCCACAGTTTTTGTGGCAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCCGGGAC‑AGCCTGTACAGGTat                                                                                                                                                                                                                                                       >  2:330973/1‑250 (MQ=255)      cAAAGTAGGCTTACGCAGTGAGGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTGTGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAgaga                                                                                                                                                                                                                                                  <  2:787475/251‑1 (MQ=255)            aGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTc                                                                                                                                                                                                                                           <  1:34177/251‑1 (MQ=255)              gCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCAt                                                                                                                                                                                                                                         <  1:189301/251‑1 (MQ=255)                 tacgcctggtggttgatcccTTGCGTGGCATCCGTGAAGTCGTGCCTCGTCTGCTGGAAATCTTGAGTAAGAATAATATTCAGGCCAGTATTTTTTTCAGCGTCCGCCTGTAAAATATGGGCCGCCATCTCTTGCGACTGTTGAAGCCACGGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGTATCGCTTTATGGCTGGGATATTTGCCTGGCAGTTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGc                                                                                                                                                                                                                                       <  2:803965/232‑1 (MQ=255)                 tACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATg                                                                                                                                                                                                                                        <  1:223818/249‑1 (MQ=255)                  ccGCATTGATGGCGTGATCTTTCGGTGCACCCGGGAAGGCGTGTCGCTGCGGCTGGAAATATTGAGTAAGCATAATAGTCAGGCCAGCATGTTTGTCAGCGCCGGCCCGGACAATATGGGCCGCCATCTCTGGGGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCGTTATGGCTGGGATATTTTACGGGCAGGTAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGcc                                                                                                                                                                                                                                      <  2:525411/249‑1 (MQ=255)                    caatttaTGTCTATACCTTGCGTGGCCCCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGATTAAGCAAAAGATTCAGGCCAGCATTTTTGTCAGCGTCGTCCCGGACAATATGGGCCGCAATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGTTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGTTCATGCCAAt                                                                                                                                                                                                                                   <  2:53706/245‑1 (MQ=255)                      aTTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTTCTGGCACGGAC‑GGCCGGGGCAGGGAAagaggtttgggattgacagac                                                                                                                                                                                                                                  >  2:197044/1‑234 (MQ=255)                      aTTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGAAAAAGAT‑TGGTCATGCCaaagc                                                                                                                                                                                                                                 >  1:211656/1‑248 (MQ=255)                              cGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCCTGTCAATACCGataacc                                                                                                                                                                                                                         >  1:480489/1‑248 (MQ=255)                              cGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTGCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCAACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCa                                                                                                                                                                                                                         <  2:449250/251‑1 (MQ=255)                                aTACCTTTCGTGGAACCCGTGAAGGCGGGCCGCGTCTGATGGAAATCTTGATTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCAtt                                                                                                                                                                                                                       <  2:309641/251‑1 (MQ=255)                                       tCGTGGCTCCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGaa                                                                                                                                                                                                                 <  2:66913/250‑1 (MQ=255)                                       tCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGaa                                                                                                                                                                                                                 <  1:158741/250‑1 (MQ=255)                                        cGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGc                                                                                                                                                                                                               <  1:439718/251‑1 (MQ=255)                                         gTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATCATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAA‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGGCCATTCCGATATAATTCGTGAAg                                                                                                                                                                                                                >  1:383419/1‑249 (MQ=255)                                         gTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGGAAAAGGA‑TGGTCATGCCAATGGCGATATTATTCGTGAAGCg                                                                                                                                                                                                              >  1:479543/1‑251 (MQ=255)                                                 cGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGTGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGGCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACa                                                                                                                                                                                                      <  2:332594/251‑1 (MQ=255)                                                  gTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAAAGAT‑TTGTCATGCCAATGCCGATATCATTTCTGAAGCGGCTAAACa                                                                                                                                                                                                      >  1:732663/1‑250 (MQ=255)                                                               gcgTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGTCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGcc                                                                                                                                                                                        <  2:751516/251‑1 (MQ=255)                                                                 gTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTg                                                                                                                                                                                      <  2:136127/251‑1 (MQ=255)                                                                 gTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGGAA‑GGGCTGGCCAGGGAAAGAGAT‑TGGTCCTGCCAATGGCGATAGCATTTGTGAAGGGGGTAAAAATAACGaaagcggcatg                                                                                                                                                                                      >  1:100626/1‑242 (MQ=255)                                                                                tGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCa                                                                                                                                                                        <  1:497985/250‑1 (MQ=255)                                                                                 gAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCAcc                                                                                                                                                                      <  1:700996/251‑1 (MQ=255)                                                                                   gTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCAt                                                                                                                                                                    <  2:339949/251‑1 (MQ=255)                                                                                           aaatttCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAAGATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCa                                                                                                                                                            <  2:523091/247‑1 (MQ=255)                                                                                                 cAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCg                                                                                                                                                      <  1:142427/251‑1 (MQ=255)                                                                                                        gCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGGCGGCCTGCTCGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGt                                                                                                                                               <  1:674678/251‑1 (MQ=255)                                                                                                          aTTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGGCCGGCCCGGGAAAGAGAT‑TTGTCATGCCCATGCCGATATCATTTTTGAAGAGGGCAAAAAACACAAAGGCGGCCCGCACGCCTGGGGTCACCaatcccggcaaggccggaggggga                                                                                                                                              >  1:476701/1‑227 (MQ=255)                                                                                                               tttCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGGAC‑GGCCCGGCCAGGGAAAGAAAA‑TTGTCATGCCAATGCCGATATCATTTGTGGAGCGGGTAAAACTCACGAAGACCGCCTGCACACCTGGGATCACCATGCCCGGCCAGCCCGTTGCGGTAAATggg                                                                                                                                        >  1:764895/1‑251 (MQ=255)                                                                                                                  caggtgcGGCCCGGAATATAGGGGCCGCCAGCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCCCGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATc                                                                                                                                     <  2:139232/245‑1 (MQ=255)                                                                                                                  cAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATc                                                                                                                                     <  2:337865/251‑1 (MQ=255)                                                                                                                    gCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGGGGAAGATGTTGCGCTCAAACGCGGCATCGCTTTATGGCTGGTATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGTTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGCATCgg                                                                                                                                   <  1:186584/251‑1 (MQ=255)                                                                                                                        cGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGAGGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCaaa                                                                                                                               <  2:443536/251‑1 (MQ=255)                                                                                                                                aCAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTAGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATc                                                                                                                       <  2:211656/251‑1 (MQ=255)                                                                                                                                   atatGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGCCCAAGCCGATATCCTTCGTGAAGCGGGTAAAAATCACGAAGACGGCCCGCACGCCCGGGATAACCATGCCCAGCAAACCCGTCGCGGGAACTGGGATCGGCAAAAAATGATcgtc                                                                                                                    >  1:662234/1‑249 (MQ=255)                                                                                                                                        ggCCGCCAGCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTGATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGatat                                                                                                               <  2:635125/251‑1 (MQ=255)                                                                                                                                           cgccATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGGAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCCCGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAAATGGGATCGGCAAACAATGATCGACGCTATTGc                                                                                                            >  1:846928/1‑251 (MQ=255)                                                                                                                                                    tGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGCCAATGCCGGTATCATTCTTGAAGCGGCTAAACAACACGAAGGCCGCCTGCACGCCTGGGGTCACCAAGCCTGGCAAGCCCGTAGCGGGAACTGGGGTCGGCAAAAAAAGATCGACGATATTGCACGCGGTCt                                                                                                   >  1:808193/1‑251 (MQ=255)                                                                                                                                                        gACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGGAA‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGGCGGCCTGCAAGCCTGGGATCACCATGCCTGGCAAACCCGTAGCGGTAACTGGGGACGGCAAACAATGATCGGCGATATTGCACACGGTCTTCGc                                                                                               >  1:798876/1‑251 (MQ=255)                                                                                                                                                            ggTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCAct                                                                                            <  1:418350/250‑1 (MQ=255)                                                                                                                                                                           tttGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATc                                                                             <  1:831231/250‑1 (MQ=255)                                                                                                                                                                               tgGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGAAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGGAGAGGCTAAACATCACGAAGGCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGGAACTGGGGTCGGGAAACAATGATCGACGATATTGCACGCGGGCTTCGCACTCTGGAAGAGGATATCgggc                                                                         >  1:539789/1‑248 (MQ=255)                                                                                                                                                                                ggAAGATGCGGCGCTCAAACGCGGCATCGTTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCaa                                                                       <  2:732663/251‑1 (MQ=255)                                                                                                                                                                                     aTGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGGCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCCGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGGTATTGCACGCGGTCTTCGCACTCTGGAAGAAGTTATTGGGCAACCggg                                                                  >  1:632642/1‑250 (MQ=255)                                                                                                                                                                                          gcgcTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑CGCCTGGGCAGGGAAAAAGATATGGT‑ATGCCCATGCCGGTATCTTTCATGTAGCGGCTAAAAATCACGCAGTTGGCCTGCACCACCGGGGTCACCAAGCCCGGCCAGCCCGGAGCGGGAACTTGGATCGGCAAACAAAGGTCGACGAATTTTCACGCGGGGTTCGCCCTCTGGGAGGGGTTATCGGACCACCGGTAACCt                                                             >  1:79302/1‑251 (MQ=255)                                                                                                                                                                                           cgcTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCCGGCCAGGGAAAGAGAT‑TGGGCATGCCAATGCCGATATTATTCGTGAAGCGGCTAAAAATCACGAAGTCGGCCTGCACGCCTGGGATCCCCATGCCCGGCAAGCCCGTAGCGGTAAATGGGATCGGCAAACAATGATAGCCGGTATTGCACGCGGGCTTCGCACTCTGGAAGAGATTATCGGTCAACCggagacctg                                                            >  1:419161/1‑244 (MQ=255)                                                                                                                                                                                               cAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGATATTATCGGTCAACCGGTAACCTGTTCt                                                        <  2:156258/251‑1 (MQ=255)                                                                                                                                                                                                aaaCGCGGAATCGTTTTATGGTTGGGATATTGTTCGGGCAGGTAC‑GGCCGGGGCAGGTAAAGAGAT‑TGTGCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTg                                                       <  2:690152/251‑1 (MQ=255)                                                                                                                                                                                                   cgcgGCATCGCTTTATGGATGGGATATTTTACTGGCAGGTTC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCg                                                    <  2:475857/251‑1 (MQ=255)                                                                                                                                                                                                      ggCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGGAAAGAGAT‑TGGTCATGCCAATGCCGATATAATTCGGGGAGCGGCTAAACATCACGAAGTCGGCCTGCAAGCCCGGGGTCACCATGCCCGGCAAGCCCGGAGCGGTCACTGGGATCGGGAAAAAATGATCGGCGATATTGCACGGGGGTTTCGCACTCTGGGAGGGGTTATCGGGCAAACGGGAACCTGTTCTGCCGCTg                                                 >  1:813974/1‑251 (MQ=255)                                                                                                                                                                                                       gcagcGCTTTATGGCTGGGAGTTTTTACTGGCAGGTAT‑GGCCTGGCCAGGTAAAGAGAT‑TGGTAATGCCAAGGCCGATATCATTCGGGAAGCGGCTAATCATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCTCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGc                                                <  2:480489/247‑1 (MQ=255)                                                                                                                                                                                                          tCGCTTTATGGCTGTGTTATTTGACTGGC‑CGTGCTGGCATGGCCAGGTAAAGAGAT‑TGGTCCTGCCAATGCTGTTTTTATTCCTGAAGCTGCTAAACCTCACGAAGTCGGCCTGCCCGCCTGGGAGCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGTTCGGCAAACAATGATCGTCGATATGTCACGCGGTCTTCGCACTCTGGAAGATATTAGCTGTCACCCGGTAACCTGGTCttcgcctgcgg                                              <  2:781302/250‑11 (MQ=255)                                                                                                                                                                                                            gCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATTCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCCGGGGTCACCATGCCTGGCCAGCCCGTAGCGGGAACTGGGATCGGCAAACAATGATCGACGGTATTGCACGCGGGCTTCCCACTCTGGAAGAGATTTTCGGTCAACCGGTAAACTGGTCTTCCGCTTCggggc                                           >  1:401916/1‑249 (MQ=255)                                                                                                                                                                                                                gatgtcTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCTGCTAAACAGCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTTGAAGAGATTAGCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGc                                        <  2:398766/245‑1 (MQ=255)                                                                                                                                                                                                                  tGGCTTGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTTGGATCGGCAAAAAATGATAGAAGGTATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTTCGCCGCTGCGGGCTGGCGTg                                     >  1:690516/1‑251 (MQ=255)                                                                                                                                                                                                                                 aCTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAATCTGGTCTGCCGCTGCGGGCCGGCGTGCCGACCAGAAGGTgg                      >  1:801536/1‑250 (MQ=255)                                                                                                                                                                                                                                 aCTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTg                       <  2:539789/250‑1 (MQ=255)                                                                                                                                                                                                                                    ggCAGGTCC‑GGCCTGGCCAGGTAAAGAGAT‑TGGCCATGCCAATGCCGATATCATTCGGGAAGCGGATAAACATCACGAAGTCGGCCTGCACGCCTTGGATCCCCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTTCCGACCAGAAGTTGATc                    <  2:674196/250‑1 (MQ=255)                                                                                                                                                                                                                                    ggCAGGTAC‑GGCCTGGCCAGGTAAATAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCg                   <  2:669570/251‑1 (MQ=255)                                                                                                                                                                                                                                    ggCAGGTAC‑GGACTGGCCAGGGAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGTAGCTGCTAAACATCACGAATGCGGACTGCACGCCTGGGAGCACCATGCCTGGCAAGCCCGTAGCGGTAACCGGGATCGGCAAACAATGAACGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTTCCGATGCGGGCTGGCGTGCAGACCAGAAGGTGATCg                   <  2:721755/251‑1 (MQ=255)                                                                                                                                                                                                                                     gCAGGTTC‑GGCCTGGCCAGTTAAAGAGAT‑TGGGCATGCCAATGCCGGTATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCCCCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGa                  <  2:479543/251‑1 (MQ=255)                                                                                                                                                                                                                                     gCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATTTCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTTTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGa                  <  1:605347/251‑1 (MQ=255)                                                                                                                                                                                                                                      cAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGa                  >  1:861515/1‑250 (MQ=255)                                                                                                                                                                                                                                        ggTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGc               <  1:86050/251‑1 (MQ=255)                                                                                                                                                                                                                                         gTAC‑GGCCTGGCCAGGTAACGAGAT‑TGGTCATGCCAATGCCGCTATCATTCGTGAAGCGGCTAAACATCACGACGTCGGCCTGCACGCCTGGGAACACCATACCTGGCAAGCCCGTAGCGGTAAATGGGATCGGCAAACAATGATCGACGATATTGAACGCGGTCTTCGCACTCCGGAAGCGATTTTCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTTGCTGGCCGACCAGACGGGGATCGaacc               >  2:354725/1‑248 (MQ=255)                                                                                                                                                                                                                                         gTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGAAGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGGGCCGACCAGAAGGGGATCGAAGc               >  2:231896/1‑250 (MQ=255)                                                                                                                                                                                                                                           aC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGCaaa            <  1:330973/251‑1 (MQ=255)                                                                                                                                                                                                                                           aC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGCaa             <  1:426943/250‑1 (MQ=255)                                                                                                                                                                                                                                                ccTGGCCAGGTAAAGAGAT‑GGGTCAGGCCAATGCCGTTATCAGTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGCAAaagaa        <  2:686023/251‑1 (MQ=255)                                                                                                                                                                                                                                                      cAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGCAAAAGAAGCGTTc  <  2:199793/251‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        ATGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATGGGCCGCCATCTCTGGCGACTGGTGAAGCCACAGTTTTTGTGGAAGATGCTGCGCTCAAACGCGGCATCGCTTTATGGCTGGGATATTTTACTGGCAGGTAC‑GGCCTGGCCAGGTAAAGAGAT‑TGGTCATGCCAATGCCGATATCATTCGTGAAGCGGCTAAACATCACGAAGTCGGCCTGCACGCCTGGGATCACCATGCCTGGCAAGCCCGTAGCGGTAACTGGGATCGGCAAACAATGATCGACGATATTGCACGCGGTCTTCGCACTCTGGAAGAGATTATCGGTCAACCGGTAACCTGTTCTGCCGCTGCGGGCTGGCGTGCCGACCAGAAGGTGATCGAAGCAAAAGAAGCGTTC  >  NC_000913/2370018‑2370512

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