Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,878,973	A→C	16.1%	R382R (CGT→CGG)	gyrB ←	DNA gyrase, subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,878,973	0	A	C	16.1%	133.9 / 2.7	57	R382R (CGT→CGG)	gyrB	DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (10/34);  new base C (7/4);  total (17/40)
Fisher's exact test for biased strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTT  >  NC_000913/3878745‑3879209                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 gAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTcacgggcacggg                                                                                                                                                                                                                                        <  2:1306775/239‑1 (MQ=255)   aGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCAcgggcacgggca                                                                                                                                                                                                                                      <  1:949030/239‑1 (MQ=255)       cgcTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTTACGCGCGCGGCGCGCCG‑CCTCCCGGGGCCGGgcagcc                                                                                                                                                                                                                                   >  1:249868/1‑237 (MQ=255)         cTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTTTCGCGCGCGCCGCGCCG‑CCTCCCGGGCACGCGCGGCACCg                                                                                                                                                                                                                                >  1:1251898/1‑239 (MQ=255)            tctcGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTcacgggc                                                                                                                                                                                                                                             <  1:312726/223‑1 (MQ=255)                   gTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGc                                                                                                                                                                                                                      <  1:1159707/239‑1 (MQ=255)                      gAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCGTTTTTCCCGCCCCCCGCCCCC‑CTCCCCGCGCACCGgcagccaccaccccctttcca                                                                                                                                                                                                                    >  2:312726/1‑222 (MQ=255)                            tttACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCCTTTTACGCGCGCGGCGCGCCG‑CCTCCCGGGGACGGGGCGCCCCGCGAATTTTTCCAACCACg                                                                                                                                                                                                              >  1:758664/1‑238 (MQ=255)                                   ttCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcgc                                                                                                                                                                                                      <  1:41455/239‑1 (MQ=255)                                    tCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTGCTTTCACCCCCGCGACGCGCCG‑CCTCAACGGCACGGGGAGGATCGATAATTTTTCCAACCACCATTTTcgcc                                                                                                                                                                                                     >  1:345746/1‑238 (MQ=255)                                     cAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGt                                                                                                                                                                                                    <  1:1278762/239‑1 (MQ=255)                                     cAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCCTTTCACCCGCGCGGCGCGCCG‑CCTCCCGGGCCCGGGGCGCATCGATAATTTTTCCCACCACCATTTTcgcg                                                                                                                                                                                                     >  1:1468455/1‑238 (MQ=255)                                            aGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCTACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑Tg                                                                                                                                                                                             <  1:963419/238‑1 (MQ=255)                                                  gCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTT‑CACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTcc                                                                                                                                                                                     <  1:276116/239‑1 (MQ=255)                                                      ggTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCcagc                                                                                                                                                                                  <  1:247846/239‑1 (MQ=255)                                                      ggTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCcagc                                                                                                                                                                                  <  1:49552/239‑1 (MQ=255)                                                        ttCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCcagca                                                                                                                                                                                 <  1:825128/238‑1 (MQ=255)                                                          cTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGt                                                                                                                                                                              <  1:501869/239‑1 (MQ=255)                                                            tGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTAt                                                                                                                                                                            <  1:799202/239‑1 (MQ=255)                                                                 ttACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGGCTTTTAAGCGCCGGCCCCGCCGCCTTCCCGGG‑GCGGGGCGCAACGGCAATTTTTTCAACCACCACTTTTTCCGCGGTGTGGGTTTTCcacccgggttttt                                                                                                                                                                         >  2:1446006/1‑227 (MQ=255)                                                                  tACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGcc                                                                                                                                                                      <  2:100628/239‑1 (MQ=255)                                                                  tACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGcc                                                                                                                                                                      <  2:1198107/239‑1 (MQ=255)                                                                  tACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGc                                                                                                                                                                       <  2:50841/238‑1 (MQ=255)                                                                     gCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCcagca                                                                                                                                                                                 <  1:1176423/225‑1 (MQ=255)                                                                      agagccagtCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCcagca                                                                                                                                                                  <  1:161030/231‑1 (MQ=255)                                                                        ccTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAgt                                                                                                                                                                <  1:1052727/239‑1 (MQ=255)                                                                         cTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAgtt                                                                                                                                                               <  1:1537835/239‑1 (MQ=255)                                                                              tCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgttc                                                                                                                                                          <  1:159149/239‑1 (MQ=255)                                                                                  agagCTGCGCGCGGAGTCCCCTTCCCACAGGTACAGTTCGGCAGTCGCCGGATCGCGGTCCTGGCGTGCTGCCAGTTTGCCGCGC‑AGGCCTGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGGGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATc                                                                                                                                                       <  2:746181/238‑1 (MQ=255)                                                                                  agagCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTTGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATc                                                                                                                                                       <  1:1446006/238‑1 (MQ=255)                                                                                  agagCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATc                                                                                                                                                       <  1:1369052/238‑1 (MQ=255)                                                                                  agagCCGCCCGCGGAGTCCCCGTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATc                                                                                                                                                       <  2:19807/238‑1 (MQ=255)                                                                                                          ccccggttccAGTGCGGAAGCCGCCGGACGGCGTTCCTGGCAGTCTGCCCGTTTGCCCGGC‑GGGCCGGCTAGGCCGGGCCCCCCTTTCCGGGGGTTCTTTTCCCGCGCGCGACGCGCCG‑CTTCCGGGGCACGGGCAGCATCGATAATTTTCCCACCCCCGATTTCCGCGTGGTT‑TGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCt                                                                                                                            <  2:684450/230‑1 (MQ=255)                                                                                                                    tgagcggaaggctggggcgcgcgTCCCTGGATGTCGGCCAGTTTCACCGAG‑AGGCCCGGTAAGTAGAGGGCCTCTTGAGGGGGGTTCATTTCACGCGCGCGACGCGCCG‑TTGCACGGTCCCGGGCAGCTACGATAATTTTGCCAACCACGTTTTTCGCCGCGGT‑TGGTTTTCCCGGCAGGTATTCTGACAGCAGTTCGTTCTTCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAAc                                                                                                                  <  2:686551/220‑1 (MQ=255)                                                                                                                       tCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGaaa                                                                                                                   <  1:1265995/235‑1 (MQ=255)                                                                                                                         ggAAAGCGCGGGATCGCGTCCCTGGCAGTCGGCCAG‑TTGCCCGGCGGGGCCTGCAAGGGCGAGCGCCCCTGTACGGCGGTCCATTTCGCGCGCGCGACGCGCCG‑CTACCGGGGCGGGGGCAGAGTCAATAATTTTGCCACCCCCGATTTTCGGGTCGTT‑TGGGTTTTCAAGCGGGTTTCCTGCCAGAGTTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGt                                                                                                              <  2:714939/238‑1 (MQ=255)                                                                                                                               cgcggtgtCGCGTCCCTGGCGGGCTGCCGGTTTCCCCGGC‑AGGCCCGCTAAGGCGAGCGCCCCTTTACGGCGGGTCATTTCACGCGCGCGCCGCCCCG‑CTTCCCGGGCACGGGCAGCATCGATAATTTTGCCACCCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCt                                                                                                        <  2:1332024/231‑1 (MQ=255)                                                                                                                               cgcgGGGTCGCGTTCCTGGCCTGCGGCCAGTTTCCCCGGC‑AGGCCCGCTAAGTGGGGCGCCCCTTTGCGGCGGGTCATTGCACGCCCGCGCCGCGCCG‑TTTCAGGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTGGTT‑TGGGTTTTCCAGCAGGTTTTCTGCCGGCAGTTCTTTCTTCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCt                                                                                                        <  2:1239168/238‑1 (MQ=255)                                                                                                                                    gATCGCGTTCCGGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCCCCCCTTTACGGCGGGTCATTTCCCGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGt                                                                                                   <  1:990554/238‑1 (MQ=255)                                                                                                                                    cgcgcgtttCCGGGGAGCCGGCGAGTTCGCCGGGC‑GGCCCCGTAAGGCCGGGCCCCCCTTTACGGGGGGTCTTTCCCCGCGCGCGACGCGCCG‑CTTCACGGGCAGGGGCAGCATCGATAATTTTGCCACCCACGTTTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCGGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGt                                                                                                   <  2:1359115/231‑1 (MQ=255)                                                                                                                                                              ttccccGGC‑AGGCCCGCAAAGGCGGGCGCCCCTTTCCGGCGGTTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCgg                                                                          <  1:1287758/233‑1 (MQ=255)                                                                                                                                                               tGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCAATTCCAGCGCGCGGCGCGCCG‑CCTCACCGGCACGGGGCGCACCGGTAAAATTTCCAACCACGATTTTTTCCTCGGG‑TGGGGTTTTCTTCCGGGATTCTTCCAGCCGTTCCTTCCTCTGA‑TGTTTCAACTCCCGTCTTTCCCCCGAACGAAACC‑CTTTGTCTTTTGTTTGGGTGGGGGAATTTTGGGTcggcg                                                                         >  2:1135960/1‑233 (MQ=255)                                                                                                                                                               tGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCATTTTACGCGCGCGGCGCGCCG‑CCTCCCCGGCGCGGGCGGCGACGGTAATTTTTCCCCACACCCCTTTTTGGTCGGG‑TGGGGTTTTCCGCCGGGATTTTGCCCCCCGTTCGTTCATTCGCTTT‑TTTAACCCCCG‑TTTTGCCCCACAAGGAAACCGCTTGTCTTTGTTGTGTG‑GGGGGAGATTCGGGTCCGGcg                                                                        >  2:683340/1‑236 (MQ=255)                                                                                                                                                                                aaaGTCGAGCCCCCCTTTCGGGCGGGTCTTTTCACGCGCGCGCCGCGCCG‑CTTCGCGGGCACGGGCACCATCGATATTTTTGCCAACCACGTTTTTCGCGTCGGT‑GGGGTTTTCAAGCAGGTATTCTGCCAGAAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGAcc                                                       <  2:345746/237‑1 (MQ=255)                                                                                                                                                                                    tCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTAAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCttt                                                                                                      <  1:131906/188‑1 (MQ=255)                                                                                                                                                                                     cctgcgccctttgccggcggGGCATTCCCCGGGCGCGCCGCCCCG‑CTGCACGGGCACGGGCAGCAGCGATAATTTTGCCAACCGCGTTTTGCGCGTGGGT‑GGGTTTTTCCAGCGGGTATTCTGCGAGCAGTTCGTTCATTTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATc                                                 <  2:758664/225‑1 (MQ=255)                                                                                                                                                                                         cgcCCCTTTACGGCGGGTCATTTCCCGCGCGCGCCGCGCCG‑CTTCCCGGCCACGGGCAGCATCGATAATTTTGCCACCCCCGATTTTCGCGTCGGT‑TGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCTTTCATCTGT‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAgg                                              <  1:351277/238‑1 (MQ=255)                                                                                                                                                                                                  ccggcggGTCATTCCCCGCGCGCGCCGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCAcgcg                                    <  1:360631/238‑1 (MQ=255)                                                                                                                                                                                                              tccACGCGCGCGACGCGCCG‑CTGCGCGGGCGCGGGCAGCATCGATAATTTTCCCAACCACGATTTGCGCGGCGGT‑TGGTTTTTCCAGCAGGTATCCTGCCAGAAGTTCGTTCATCTGC‑TG‑TCAAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCgg                        <  1:67309/237‑1 (MQ=255)                                                                                                                                                                                                                  cgcgcgcgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgc                    >  1:1060035/1‑239 (MQ=255)                                                                                                                                                                                                                  cgcgcgcgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGGCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgc                    >  1:1177275/1‑239 (MQ=255)                                                                                                                                                                                                                   gcgcgcgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgc                    >  1:345696/1‑238 (MQ=255)                                                                                                                                                                                                                   gcgcgcgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGCTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTTACGCGCGTCGTCAACGgggga                    >  2:731530/1‑237 (MQ=255)                                                                                                                                                                                                                    gcggccgcgggtcgcctgcacgggcacgggCCGCATCGAAAATTTTGCAAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgcg                   <  2:1251898/221‑1 (MQ=255)                                                                                                                                                                                                                    cgcgcgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgc                    >  1:1233005/1‑237 (MQ=255)                                                                                                                                                                                                                        cgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGCTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGa               >  1:84226/1‑238 (MQ=255)                                                                                                                                                                                                                        cgACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTAGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTg                >  1:805222/1‑237 (MQ=255)                                                                                                                                                                                                                           ggccccg‑cTTCAGGGGCAGGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACtt            <  1:1354102/234‑1 (MQ=255)                                                                                                                                                                                                                           cgcgCCG‑CTTCACGGCCCCGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGCGGGT‑TGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCGCC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACttt           <  1:831507/239‑1 (MQ=255)                                                                                                                                                                                                                                    ttCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCttttt   >  1:1337972/1‑239 (MQ=255)                                                                                                                                                                                                                                      cacgggcacgggCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCtttttt  >  1:1340511/1‑238 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 GAAGCGCGCTTTCTCGACGTTGAGGATTTTACCCTTCAGCGGCAGAATCGCCTGGTTCTTGCGGTTACGCCCCTGCTTCGCAGAGCCGCCCGCGGAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGC‑AGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCG‑CTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGT‑TGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGC‑TG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGG‑AGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTT  >  NC_000913/3878745‑3879209

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