Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,011,802	C→T	P243S (CCG→TCG)	bioP →	biotin transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,011,802	0	C	T	90.9%	181.1 / 4.0	67	P243S (CCG→TCG)	bioP	biotin transporter
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base T (33/27);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.99e-01

CCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGC  >  NC_000913/4011555‑4012047                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ccTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCgg                                                                                                                                                                                                                                                     >  1:405732/1‑251 (MQ=255)   tGATGGAAACTCGCCCGATGCCACCGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGCAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGGGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGGGGAACAACGGCGCGACGCAGGAGGGCGCCGGAACGCTGGGCATTATGAATAATCTGCACGTTTCGGc                                                                                                                                                                                                                                                    >  1:241669/1‑250 (MQ=255)   ggATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGGGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCa                                                                                                                                                                                                                                                   <  2:233277/250‑1 (MQ=255)     aatggaacctggccGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAgg                                                                                                                                                                                                                                                 <  1:440021/239‑1 (MQ=255)     aTGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAgg                                                                                                                                                                                                                                                 >  1:259632/1‑251 (MQ=255)       ggAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGgc                                                                                                                                                                                                                                               <  2:162668/251‑1 (MQ=255)        gAAACTCGCCCGCTGCCACAGCATAACGCCTTTGCGGGGTTCTATCTGGGAGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAAGGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCGTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGTCGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATTCATGTTTCGGCAGGgct                                                                                                                                                                                                                                              <  2:279891/251‑1 (MQ=255)               gCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGt                                                                                                                                                                                                                                        >  1:150062/1‑250 (MQ=255)                 ccGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGCTTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTaa                                                                                                                                                                                                                                      >  1:66707/1‑250 (MQ=255)                 ccGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTaa                                                                                                                                                                                                                                      >  2:123659/1‑250 (MQ=255)                             cATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTg                                                                                                                                                                                                                         >  1:96109/1‑251 (MQ=255)                               tAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGc                                                                                                                                                                                                                       >  1:285998/1‑251 (MQ=255)                                                tCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTGCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCGTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCGGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGAACTg                                                                                                                                                                                                       <  2:69853/250‑1 (MQ=255)                                                            cgTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttatt                                                                                                                                                                                          <  1:43831/251‑1 (MQ=255)                                                            cgTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGACTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttatt                                                                                                                                                                                          >  1:140081/1‑251 (MQ=255)                                                                  tggtggCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCGATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCATCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCCGGGGGACGCCGGAAAGCTGGGCATTGTGAATAATATGCACGTTTCGGCAGGGCGGCTGGTACACCTGGCTATCTGGCCCCAACAGGCGCACGGGCCAACGGTTATTAAAggg                                                                                                                                                                                    >  2:781/1‑250 (MQ=255)                                                                            gATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTg                                                                                                                                                                           <  2:259632/250‑1 (MQ=255)                                                                            gATTGCATGGTGCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATGGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGa                                                                                                                                                                          <  2:89559/251‑1 (MQ=255)                                                                             aTTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGAt                                                                                                                                                                         >  1:97824/1‑251 (MQ=255)                                                                              atGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACAACGCTGCAATGGGGCATTCTGGCGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTAAGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCGGGCGGTAAAACTGGCTATCTTGAACCAAAAACCGGACTGGCCAACGATTATTCCAGgcacgacgggatac                                                                                                                                                                        >  2:354366/2‑239 (MQ=255)                                                                                    ggTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGcc                                                                                                                                                                  <  1:123659/251‑1 (MQ=255)                                                                                    ggTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGc                                                                                                                                                                   <  1:338935/250‑1 (MQ=255)                                                                                     gTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAAAAGCCGCACTGGCCAACGTTTATTACCGGCGCGCTGGTGATCCTGGccc                                                                                                                                                                 >  2:438482/1‑250 (MQ=255)                                                                                             tGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTgg                                                                                                                                                         >  1:277293/1‑251 (MQ=255)                                                                                                gAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGTGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGGGGAACTACGGCGCGTCGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCGATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTg                                                                                                                                                      <  2:349623/251‑1 (MQ=255)                                                                                                   aTGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCAt                                                                                                                                                   <  1:318523/251‑1 (MQ=255)                                                                                                    tGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATc                                                                                                                                                  >  1:369729/1‑251 (MQ=255)                                                                                                      cgcAGAAAATGCCGCAAACCACGCTGCAATAGGGCATTCTGGTGATTCTTGGCGTAGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAACCGAACTGGCCAACGATTATTACAGGCGCGCTGGTTATCCTGGCCTCACTGTGGGTGCATCGt                                                                                                                                                >  1:195669/1‑251 (MQ=255)                                                                                                                          aCGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCCACAGCCGCACTGGCCACCGTTTATTACAGGCGCGGTGGAGAACCTGGCCTCACTGTGGGGGCATCGGAAGGGGGTCGCTCCGCGctc                                                                                                                            >  2:381392/1‑251 (MQ=255)                                                                                                                            gCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCtt                                                                                                                          <  1:270249/251‑1 (MQ=255)                                                                                                                                           tcgggtttgtCTTGGCGTGGTGTCTTCCGGGATTTGCGACTTGATGGGGAACTCCGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCCCAAACGGCAgat                                                                                                            <  2:286314/241‑1 (MQ=255)                                                                                                                                           tCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAgat                                                                                                            <  1:456359/250‑1 (MQ=255)                                                                                                                                                   ttCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAg                                                                                                    <  1:142559/250‑1 (MQ=255)                                                                                                                                                    tCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGc                                                                                                  >  1:428476/1‑251 (MQ=255)                                                                                                                                                     cTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGACGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGccgggg                                                                                                 >  2:68401/1‑246 (MQ=255)                                                                                                                                                            tggtggCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTg‑c                                                                                         >  1:229054/1‑251 (MQ=255)                                                                                                                                                             ggtggCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTg‑cg                                                                                        >  1:65775/1‑251 (MQ=255)                                                                                                                                                                    tcagggATTGGCTCCTTTATGTGGAACTACGGCGCGACGCAGGGGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGGGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCCCAAACGGCAGCTGTTCTCCGGCGGGGTTg‑cgcgccg                                                                                   >  2:185960/4‑247 (MQ=255)                                                                                                                                                                    tCCGGGATTGGCTACGTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTg                                                                                   <  1:436119/249‑1 (MQ=255)                                                                                                                                                                             ggCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAAc                                                                        >  1:242750/1‑251 (MQ=255)                                                                                                                                                                             ggCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAAc                                                                        <  1:148528/251‑1 (MQ=255)                                                                                                                                                                              ggttcTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAAAAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACg                                                                       <  2:150323/247‑1 (MQ=255)                                                                                                                                                                              gcttcTTTCTGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACg                                                                       <  2:405732/247‑1 (MQ=255)                                                                                                                                                                                  cTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACATGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCACATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTc                                                                   >  1:222856/1‑251 (MQ=255)                                                                                                                                                                                      aTGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTa                                                               <  1:781/251‑1 (MQ=255)                                                                                                                                                                                          ggAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACGGGCCAACGTTTATTACAAGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCAGCGCTCTTCACAAACGGCAGCTGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAGAGCTTCCGTAACTg                                                           >  1:273285/1‑251 (MQ=255)                                                                                                                                                                                           gAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTgg                                                          <  1:68401/251‑1 (MQ=255)                                                                                                                                                                                                     gggACGCAGGTGGACGCCGGAACGTTGTGCATTATGAAGAATATGCCCGTTCCGGAAGGGCTGCTGGTAAGCCTGGCTATCTGGCACCAACAGTCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGATGAGCGAATAAACGCTTCCGTAACTGGCTGACgctgc                                                <  2:323894/249‑1 (MQ=255)                                                                                                                                                                                                       gACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCCCTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCCCTGTGGGGGCATCGTAAGTGGGTCGCTCCGCGCTCTTCCCAAACGGCCGATGATCGCAGGCGTGATTGACGC‑CTGAGCGAATAAACGCTCCCGTACCTGGCTGACgctcccc                                              >  2:105795/1‑247 (MQ=255)                                                                                                                                                                                                       gACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTc                                              >  1:32938/1‑251 (MQ=255)                                                                                                                                                                                                        ccgagggggggcGCCGGAACGCTGGGCTTTAGGAAGAATGTGCACGTTTCGGCATGGCTGCTGGTAAACCTGGCTATCTGGCCCCAACAGCCGCACTGGTCAACGTTTATTACAGGCGCGCTGGTGAGCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGAGTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCg                                             <  2:407245/240‑1 (MQ=255)                                                                                                                                                                                                                ggACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATcgcg                                     >  1:342087/1‑251 (MQ=255)                                                                                                                                                                                                                     acGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCCCTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGCTTG‑CGCGCTGGGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCCTCGCGCAcgg                                >  2:84617/2‑251 (MQ=255)                                                                                                                                                                                                                       ggAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCAcggcg                              <  1:438482/251‑1 (MQ=255)                                                                                                                                                                                                                         aaCGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACggcggc                            <  1:406268/251‑1 (MQ=255)                                                                                                                                                                                                                          aCGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGgcggcg                           >  1:237896/1‑251 (MQ=255)                                                                                                                                                                                                                              tGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGCGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTAc                        >  1:186354/1‑250 (MQ=255)                                                                                                                                                                                                                              tGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGAGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGAGGATCGCAGGCGTGGTTG‑CGCGCTGAGCGAAAAAACGCTTCCGTAACGTGCTGACGATGCTCGCCCTCGCGCACGGCGGCGTAcc                       >  2:130346/1‑250 (MQ=255)                                                                                                                                                                                                                                    ttATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATATGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGGGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGCGCGACTAAACGCTTCCGTAACAGGCTGAAGCTGCTCGCCATCGCGCACCGCGGCGTACAGTccgc                 >  2:265254/1‑248 (MQ=255)                                                                                                                                                                                                                                     tATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGACTAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGcc                >  1:193007/1‑250 (MQ=255)                                                                                                                                                                                                                                         aataatATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGTCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGTAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACg                                                                       <  1:288558/192‑1 (MQ=255)                                                                                                                                                                                                                                            aataTGCTCGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGAt                                                                                             <  2:99803/168‑1 (MQ=255)                                                                                                                                                                                                                                            aataTGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGAt                                                                                             >  1:99803/1‑168 (MQ=255)                                                                                                                                                                                                                                             atatGCACGTTGCGGCAGGGCGGGTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCGTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCCCAAGc        <  2:66707/251‑1 (MQ=255)                                                                                                                                                                                                                                               aggCACGTTTCGGCAGGGCTGCTGGTAACCCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCt      <  2:366582/249‑1 (MQ=255)                                                                                                                                                                                                                                                tGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCtt     >  1:277878/1‑251 (MQ=255)                                                                                                                                                                                                                                                   aCGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGc  >  1:454089/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTG‑CGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAAGCCTTCGC  >  NC_000913/4011555‑4012047

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