Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,814,904	C→T	R137R (CGC→CGT)	slmA →	nucleoid occlusion factor SlmA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,814,904	0	C	T	98.2%	192.3 / NA	56	R137R (CGC→CGT)	slmA	nucleoid occlusion factor SlmA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/27);  minor base G (0/1);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CCGCATGTTCGATAGCCTGATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAACCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGCCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCG  >  NC_000913/3814661‑3815145                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   tcGCATGTTCGATAGCCTGGTTGCGTTTATCGAAGCTAGCCTGATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGACTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCCTGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTAt                                                                                                                                                                                                                                             >  2:192398/2‑251 (MQ=255)  cGCATGTTCGATAGCCTGATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCGGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTGTGAACAGGATCGCCGGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGGCAGGTa                                                                                                                                                                                                                                              <  2:288610/249‑1 (MQ=255)       gTTCGATAGCCTGATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTg                                                                                                                                                                                                                                       <  2:12530/251‑1 (MQ=255)              aGCCTGATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGTCCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTTCGTGAAAagag                                                                                                                                                                                                                                <  2:292292/251‑1 (MQ=255)                   gATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAgg                                                                                                                                                                                                                                                >  2:89492/1‑230 (MQ=255)                   gATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAgg                                                                                                                                                                                                                                                <  1:89492/230‑1 (MQ=255)                                           gATTACTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCg                                                                                                                                                                                                   <  2:338627/251‑1 (MQ=255)                                               aCTCGCATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTAGTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATg                                                                                                                                                                                                <  1:192398/250‑1 (MQ=255)                                                    caccaatcTGATGCTGAAAGAAGAGAAAGACTCCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCGTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCTTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAAccc                                                                                                                                                                                          <  2:114341/244‑1 (MQ=255)                                                    cATCAATCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAAccc                                                                                                                                                                                          <  2:32945/251‑1 (MQ=255)                                                                   gAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGcc                                                                                                                                                                              >  2:148000/1‑248 (MQ=255)                                                                   gAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGcc                                                                                                                                                                              >  2:131724/1‑248 (MQ=255)                                                                     aaGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATc                                                                                                                                                                         <  2:126017/251‑1 (MQ=255)                                                                            gaAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCt                                                                                                                                                                  <  1:103109/251‑1 (MQ=255)                                                                            gaAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCt                                                                                                                                                                  >  1:256655/1‑251 (MQ=255)                                                                             aaaGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCtt                                                                                                                                                                 >  2:237028/1‑251 (MQ=255)                                                                                 acacCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGCCCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGCAAAGAGAATGCGTGAGGGTGAAGGTTACACCAACGCTACAACCCTGCTGGCAAGCCCGATCCTGGCCTTCtgt                                                                                                                                                             >  2:243065/1‑251 (MQ=255)                                                                                            cgcCTGCGTCTGATTGTGTTGCGGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCt                                                                                                                                                  <  2:181563/251‑1 (MQ=255)                                                                                             gcCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTg                                                                                                                                                 <  1:423000/251‑1 (MQ=255)                                                                                                      tGATTGTGTTGCTGCATCACGATTTTGGTGACCGTAATCCTCGCCTAACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAAACAGCTGTTCGAGCGTCTTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGGGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGtttt                                                                                                                                        >  2:285501/1‑251 (MQ=255)                                                                                                       gATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGCCCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGtt                                                                                                                                                                                                            >  2:17868/1‑182 (MQ=255)                                                                                                       gATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGtt                                                                                                                                                                                                            <  1:17868/182‑1 (MQ=255)                                                                                                                gctgctTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCg                                                                                                                              <  1:82365/251‑1 (MQ=255)                                                                                                                          ggTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGAGCTAATGCTTGAACAGGATCGCCTGCTAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCATGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCCCCGATGAAACCATGCTGGCAAGCCAGATCCTGGCCTTCTGTGACGGTATGCTGTCAAGGCTTGTCGGCAGCGAGTTTAAAt                                                                                                                     >  2:349945/1‑250 (MQ=255)                                                                                                                          ggTTTTAGCGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGCAACGTTTTGTCCGCAGCGAGTTTAAAt                                                                                                                     >  2:89646/1‑250 (MQ=255)                                                                                                                                ggttaGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGcc                                                                                                               <  2:73202/246‑1 (MQ=255)                                                                                                                                                gCCTGACCCGCCTCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGAc                                                                                              <  2:210157/251‑1 (MQ=255)                                                                                                                                                gCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGAc                                                                                              >  2:136684/1‑251 (MQ=255)                                                                                                                                                   tGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGTTCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGcc                                                                                           >  2:37471/1‑251 (MQ=255)                                                                                                                                                    gACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGccc                                                                                          >  1:60905/1‑251 (MQ=255)                                                                                                                                                       ccGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCt                                                                                       <  1:131724/251‑1 (MQ=255)                                                                                                                                                              atCACTGGGCATGCTCTCATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACAAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGTGTGAAGGTTACACCACCCATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATGTTGACGCTCGCTgggctct                                                                                <  2:8480/250‑6 (MQ=255)                                                                                                                                                                        aTGCGCTAAGGTGTGAACAGGAGCGCCTGCAAGGGCGCACCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGca                                                                      <  1:192907/251‑1 (MQ=255)                                                                                                                                                                          gcgcTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAACCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGcat                                                                     <  2:129774/250‑2 (MQ=255)                                                                                                                                                                                   tttGAACAGGATAGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGCAAAGAGACTGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGt                                                           >  2:438994/1‑250 (MQ=255)                                                                                                                                                                                        aCAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGGGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATg                                                     <  2:141238/251‑1 (MQ=255)                                                                                                                                                                                          aGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGAc                                                   <  1:237028/251‑1 (MQ=255)                                                                                                                                                                                            gATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGACACCCTGCTGGCAAGCCAGATCATGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCAGACGGATGATTTTGACGCCCGCTGGCCGCTACTTGCGGCA‑CAGTTGCAGTAATATGACGc                                                 >  1:167443/1‑251 (MQ=255)                                                                                                                                                                                             aTCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGcc                                                <  1:136684/251‑1 (MQ=255)                                                                                                                                                                                               cGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGGAGGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGAGTTTGTCGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCg                                               <  1:438583/250‑1 (MQ=255)                                                                                                                                                                                                   tgcaaggcgacattaaacAGCTGGTCGTGCGAATTGAAGCGCAGCAGCGTCAGCGACTGCGTGCTAAGAGAATGCGTGAGGGAGAAGGTTAAACCACCGATGAAAACCTGCTGGCAAGCCAGATCATGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATCCCGCCCGACGGGTGATTTTAACGCCAGCGGGCCGCTAATTGCGGCA‑CAGTTGCAGGAATCTGACGACGGATg                                           >  2:258773/18‑250 (MQ=255)                                                                                                                                                                                                      aaGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACtt                                       >  2:342714/1‑251 (MQ=255)                                                                                                                                                                                                          gcgcATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTc                                    >  1:324606/1‑250 (MQ=255)                                                                                                                                                                                                            gcATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTCAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGCTGCAGTAATATGACGCCGGATGACTTTTCATc                                 >  2:173224/1‑251 (MQ=255)                                                                                                                                                                                                                cAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGc                             >  2:343719/1‑251 (MQ=255)                                                                                                                                                                                                                     aGCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGtt                        <  1:199840/251‑1 (MQ=255)                                                                                                                                                                                                                      gCTGTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGttt                       <  1:28531/251‑1 (MQ=255)                                                                                                                                                                                                                         gTTCGAGCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCtt                    >  1:104468/1‑251 (MQ=255)                                                                                                                                                                                                                               gCGTATTGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGTGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACg              <  2:110170/251‑1 (MQ=255)                                                                                                                                                                                                                                cGTATTGAAGCACAGCTGCGTCAGCAATTGCGTAAAAAGAGAATGCGTGAGGGCGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGTCATTCTGTGAAGGTCTGCTGTCCCGTTTTGTCCGCAGCGAGTTTAAATACCGCACGACGGATGATTTTGACGCCGGCTGGCCGCTAATT‑CAGCAACAGTTACGGTATTATGACGCCGGCTGACTTTACATCCGGCGAGATTCTTTAAACGc             >  2:309218/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttGAAGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCaa          >  1:149585/1‑250 (MQ=255)                                                                                                                                                                                                                                      gAAGCGCAGCTGCGTCAGGTATTGCGTGAGAAGAGTATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGAGCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTAGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCAGGCGAGTTTCTTTAACCGCCAAAcc       >  2:359182/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaGCGCAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAActc      >  1:117114/1‑251 (MQ=255)                                                                                                                                                                                                                                       aaGCGCAGCAGCGTCAGGTATTGCGTGAAAAGAGCATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCAGGCGCGTTTCTTTAAACGCCAAActc      >  1:82002/1‑251 (MQ=255)                                                                                                                                                                                                                                         gcgcAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCtt    <  2:93076/251‑1 (MQ=255)                                                                                                                                                                                                                                           gcAGCTGCGTCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAGTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTcg  >  1:341093/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   CCGCATGTTCGATAGCCTGATTGAGTTTATCGAAGATAGCCTGATTACTCGCATCAACCTGATTCTGAAAGATGAGAAAGACACCACAGCGCGCCTGCGTCTGATTGTGTTGCTGCTTCTCGGTTTTGGTGAGCGTAATCCTGGCCTGACCCGCATCCTCACTGGTCATGCGCTAATGTTTGAACAGGATCGCCTGCAAGGGCGCATCAACCAGCTGTTCGAGCGTATTGAAGCGCAGCTGCGCCAGGTATTGCGTGAAAAGAGAATGCGTGAGGGTGAAGGTTACACCACCGATGAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCA‑CAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCG  >  NC_000913/3814661‑3815145

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