Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,791	A→T	D516V (GAC→GTC)	rpoB →	RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,182,791	0	A	T	95.7%	146.2 / NA	47	D516V (GAC→GTC)	rpoB	RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (29/16);  minor base G (0/2);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGT  >  NC_000913/4182547‑4183035                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     atCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCATCAGGGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGaa                                                                                                                                                                                                                                                <  1:298636/251‑1 (MQ=255)          ggTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCt                                                                                                                                                                                                                                       >  2:66292/1‑251 (MQ=255)          ggTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCt                                                                                                                                                                                                                                       >  1:47244/1‑251 (MQ=255)           gTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTg                                                                                                                                                                                                                                      >  1:192044/1‑251 (MQ=255)               aGGCGAAGTCGATTATATCGACCCCCTCGGCAACCGTCTTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCTCGTTGGCCTGGTACGTGTAGAGCGTGCGGTTAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACGT‑GAGATGATCAACTCCAAGCCGATTTCCGCAGCAGTTAAAGTGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTGATGTGCCAGAACAACCCGCTGTCTg                                                                                                                                                                                                                                  <  2:289472/251‑1 (MQ=255)                     aGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTCGGCGAAATGGCGGAAAACCCGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCCACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTACGTTTATGGTCCAGAACAACCCGCTGTCTGAGAtt                                                                                                                                                                                                                             >  1:128520/1‑250 (MQ=255)                         gatgatATCGACCACCTAGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGAGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGcc                                                                                                                                                                                                                        >  2:105625/1‑250 (MQ=255)                                      accTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAAAACGCTGTCTGCGATTACGCACAAACGACGTAtc                                                                                                                                                                                                           >  2:82019/1‑251 (MQ=255)                                       ccTCGGCCACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAATTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCTGTTCCAGCCAGCTGTCTCAGTTTATTGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAtct                                                                                                                                                                                                          <  1:268616/251‑1 (MQ=255)                                          cGGCAACCGACGTATCCGATCCGTTGGCGACATGGCGGAAAACCAGTTACGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGACGATCTGGATACCCTGATGCCACAG‑GCTATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTGCTTCGGTTCCAGCCAGCTGTCTCCGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGGATCTCCg                                                                                                                                                                                                       >  1:60686/1‑251 (MQ=255)                                                      tATCCGATCAGTTGGCGCAATGGCGGAAAACCAGTTCCGCGTTGGGCTGGCACGTGTAGCGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGCTACCCTGATGCCACAG‑GATATGATCAACGCCCAGCCGATTTCCGCAGCAGTGCAAGCGTTCTTCGGTTCCAGCCAGCTGTCTCCGTTTATGGTCCAGAACAAACCGCTGTCTGAGATTACGCCCCACCGTCGAATCTCCGCACTCGcccccg                                                                                                                                                                                           >  2:279344/1‑249 (MQ=255)                                                        tccgttccgtTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGc                                                                                                                                                                                         >  2:162973/1‑251 (MQ=255)                                                                gttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                                                                                                                                                                                 <  2:192044/251‑1 (MQ=255)                                                                gttGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGACTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCCGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTCCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                                                                                                                                                                                 >  1:752/1‑251 (MQ=255)                                                                         aTGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGCATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                                                                                                                                                                        >  2:146512/1‑251 (MQ=255)                                                                                          ccGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc                                                                                                                                                       <  1:279957/251‑1 (MQ=255)                                                                                                  gCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                                                                                                                               <  1:146512/251‑1 (MQ=255)                                                                                                   ccTGGTACGTGTAGAGCGTTCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTGCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                                                                                                                              <  2:147979/251‑1 (MQ=255)                                                                                                   ccTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCGCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                                                                                                                              <  2:47244/251‑1 (MQ=255)                                                                                                    cTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                                                                                                                              >  2:233868/1‑250 (MQ=255)                                                                                                     tGGTACGTGTAAAGCGTGCGGTGAAAGAGCGTCTGTCTATGGGCGATCTGGGTACCCTCATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTATTCGGTTCCAGCCCGCTGTCTCCGTTTATGGTCCAGAACAACCCGCTGTCTGCGATTCCGCATATACGTCATATCTCCGCACTCGGCACAGGCGACCCGACCCGCGAACGAGCAGGCTTCGAAGTTCgagaagaccac                                                                                                                                            >  2:27875/1‑244 (MQ=255)                                                                                                             ggatagCGTGCGGTGAAAGGGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGAGA‑ATGATACACGCCAAGCCGCTTTCCGCAGCAGTGCAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCCGTTTATGGTCCAGAACAACCCGCTGTATGAGATTACGCACAACCTTCGTAGCTACGCAATCGGCACAGTCGGTCTGACCCGTGAACGTGAAGGTTTCGAAGTACGAGACGTAACCCCGACTc                                                                                                                                     >  2:173345/5‑250 (MQ=255)                                                                                                                  gCGTGCGGTGAAAGAGCGGCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTCCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGACTCACTACg                                                                                                                               <  2:97107/251‑1 (MQ=255)                                                                                                                           gAAAGAGCGTCTGACTCTGGGCGAGCGGGATACCCTGATGCCACAG‑GATATGATCACCGCCAAGCCGATTTCCTCAGCAGTGAAAGAGTCCTTCGGTGCCAGCCAGCGGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGGATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAt                                                                                                                      <  2:92691/251‑1 (MQ=255)                                                                                                                                                cGATCTGGCTACCCTGCTGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGCGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                                                                                   >  2:57510/1‑249 (MQ=255)                                                                                                                                                      ggttACCCTTAGGCCACAG‑GATATGATCCACGCCAAGCCGATTTCCGCAGCAGTGAAAGCGTTCGTCGGTTCTAGCCAGCTGTCTCAGTTTAAGGTCCAGAACAACCCGCCGTCTGAGATTACGCACAAACGTCGTATCTCCGCCCTCTGCCCAGGCGGTCTGACCCGTGAACGTGCTGGCTTCGAAGTTCGAGACGTACACCCTACTCAATACGGTTGCGTATGTCCAATCGAAACCCCTGAAGGt                                                                                               <  2:56729/244‑1 (MQ=255)                                                                                                                                                             cTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCgg                                                                                    >  1:76587/1‑251 (MQ=255)                                                                                                                                                                aTGCCACAG‑GATATGATCACCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                                                                                 >  1:271867/1‑251 (MQ=255)                                                                                                                                                                   ccACAG‑GATATGATCAACGCCACGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGACGCGTATGTCCAATCGCAACCCCTGAAGGTCCGAACATCGGTCTGAt                                                                              >  2:37855/1‑251 (MQ=255)                                                                                                                                                                          gATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActc                                                                        >  1:141724/1‑251 (MQ=255)                                                                                                                                                                             atGATCAACGCAAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCCGAAGGTCCGAACATGGGTCTGCTCAActctct                                                                     >  2:249630/1‑251 (MQ=255)                                                                                                                                                                               gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                                     >  2:84782/1‑249 (MQ=255)                                                                                                                                                                               gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctct                                                                     <  1:84782/249‑1 (MQ=255)                                                                                                                                                                                           gCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                       >  1:226173/1‑251 (MQ=255)                                                                                                                                                                                             cGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTCCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                       >  2:82561/1‑249 (MQ=255)                                                                                                                                                                                               aTTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                                                   >  1:265968/1‑251 (MQ=255)                                                                                                                                                                                                 tttcGCAGCAGTGAAAGAGTGCTCCGTTGCCAGCCAGCTGTCTCAGGTTATAGTCCAGAACAACCCGCTGTGTGAGATTACGCACAAACGTCGTATCTCCGCGCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACTGTCGCGTATTTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                                  <  1:206078/247‑1 (MQ=255)                                                                                                                                                                                                         cagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAATTTCGAGACGTACTCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACAGCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                         <  1:23004/251‑1 (MQ=255)                                                                                                                                                                                                         cagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                         >  1:131711/1‑251 (MQ=255)                                                                                                                                                                                                                   gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGTCCAGCACAACACGCTGTCTGAGCTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAAAGTGCAGGCTTCGAAGTTAGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGACGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTCCGCCCAGCCTAACGAATACGGCTTCCTTg                               >  2:288905/1‑251 (MQ=255)                                                                                                                                                                                                                              ccagccagCTGTCTCAGTTTCTGGTCCAGAACAACCCGCTGTCTGAGAATACACACAAACGTCGTATCTCCGCACTCGCCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCAAGACGCACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAACTCGGTCTGAACACCTCTCTGTCCGTGTCCGCACAAACTAACGACTCCGGCTTCCTTGAGACTACGTac                    >  1:148455/1‑250 (MQ=255)                                                                                                                                                                                                                              ccagccagCTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa                     <  2:79341/250‑1 (MQ=255)                                                                                                                                                                                                                               cagccagcTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATc                   >  1:94343/1‑251 (MQ=255)                                                                                                                                                                                                                               cagccagcTGTCTCAGTTTATGGTCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGATTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGCATACGGCTTCCTTGAGACTCCGTATc                   >  1:109083/1‑251 (MQ=255)                                                                                                                                                                                                                                aggcagcTGTCGCAGTTTATGGTGCAGAACAACGCGGTGTCTGAGATTACGGACAAGCGGCGTATCTCCGCGCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTGCGAGACGTACGCCCGACTCACTACGGGCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCGTGAGACGCCGTATCg                  <  2:118172/248‑1 (MQ=255)                                                                                                                                                                                                                                             cgTTTATGGTCCAGAACATCCCGCTGTCTGAGATTACGCCCAAACGTCGTATCTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGGATCGTAAAGTGACCGAc     <  1:219715/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tATGGGCCAGAACAACCCGCTGTCTGATATTACGCACAAACGCCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAGCGAAACCCCTGAAGGGCCGAACATCGGGCTGATCAACTCTCTGGCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGgt  <  1:288905/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     ATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAG‑GATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGT  >  NC_000913/4182547‑4183035

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