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	94.3%	92.3 / NA	37	P243S (CCG→TCG)	bioP	biotin transporter
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base T (20/13);  minor base G (0/2);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAA  >  NC_000913/4011558‑4012039                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               gATGGAAACTAGCCCTATGCCTCAGCATAACGCCTTTGCGTTGTTTTAGCTTGGCGCGTGTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGTAAATGCCGCAAACCACGCTGCAATGTGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACGTTATGTTGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCa                                                                                                                                                                                                                                           <  2:388157/250‑1 (MQ=255)                gATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCCGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGGGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGCATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAAcc                                                                                                                                                                                                                          >  1:350876/1‑251 (MQ=255)                    ccacggccgtgcGCCTTTGCGTGGATGTATCTTGGCGGTTTTCGGGTGGCAGGGATTGCATGGTTCTTTCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTTGGGGTTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGATCGCTGGGCATTATGAATAATATGCACGTTTCGGAAG‑TGCTGCTGGTAAACCGGGc                                                                                                                                                                                                                      <  2:347732/240‑1 (MQ=255)                    ccACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCCGGTAAACCTggg                                                                                                                                                                                                                      >  2:127791/1‑250 (MQ=255)                      acaGAATAACGCCTTTGCGGGGTTCTATCTTGGCGCGTTTATGGTGGAAGTGATTGCCTGGTTCTTGCTGGGAACTGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGAGGCTTCCGGGGTTGGATACTTTATGTGGAACTACGGCGCGACGCAGGAGGGAGCAGGAACGCTGGGCATTATGAATCATATGAACGTTCCGGCAGAGGC‑GCAGGTAAACCTGGca                                                                                                                                                                                                                     >  2:388358/1‑249 (MQ=255)                                 ccTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCa                                                                                                                                                                                                          >  1:106205/1‑250 (MQ=255)                                                ttCTTGGCGCGTTTCTGGGGGCAGTGATTGTATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTGTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGTCATTATGAATAATATGCACGTTGCGGCAG‑GGCTGCTGTTAAACCTGGCTATCTGGCACCAACATCCGCACTgg                                                                                                                                                                                             <  2:170560/247‑1 (MQ=255)                                                                    gCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTATGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCt                                                                                                                                                                      >  1:41284/1‑251 (MQ=255)                                                                          agtGCATGGTTCTTGCTGGGAAAAGCTCAGAAAATGCCGCAAACCACGCTGCAAGGGGGCATTCTGGTGTTGCTGGGCGTGGTGGCTCCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTGCGGCAG‑GGCTGCTGGGAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGAt                                                                                                                                                                <  2:412257/249‑1 (MQ=255)                                                                               aTGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCACAAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTgg                                                                                                                                                           >  1:344769/1‑251 (MQ=255)                                                                               aTGGTTCTTGCTGGGAAATGCGCAGAAAATGCAGCAAAACACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCCGGTGGACGCCGGAACGCTGGGCATTATGAGTAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCCACAGCCGCACTGGCCCACGTTTATTAAAGGCGCGCTGGTGATCCTgg                                                                                                                                                           >  1:157878/1‑251 (MQ=255)                                                                                    tCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGTATTATGAATAATATGCACGTTCCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTTCCGGCGCGCTGGTGATCCTGGCCTcc                                                                                                                                                      >  1:251224/1‑250 (MQ=255)                                                                                      ttGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACt                                                                                                                                                    >  2:379485/1‑251 (MQ=255)                                                                                        gCTGGGAAATGCGTAGAAAAGTCCGCAAACCTCGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTGCTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATCATGAAGAATATGCACGTTTCGGCAG‑GGCTGCTGGTAACCCTGGCTATCTGGCACCAACATCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACtg                                                                                                                                                   <  2:74912/250‑1 (MQ=255)                                                                                                aTGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCAt                                                                                                                                          >  1:54974/1‑251 (MQ=255)                                                                                                     cAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCACGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTa                                                                                                                                      >  1:106185/1‑250 (MQ=255)                                                                                                         aaaTGCCGCAAACCACGCGGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGTCTCCCGGGATTGGCTACTTTAGGTGGAACGACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGTCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTgg                                                                                                                                 <  2:344769/251‑1 (MQ=255)                                                                                                                gCAAACCACGCTGCAATGGGGCATTCTGGTTTTTCTTGGCGTGGTGGCTTCCGGGCTTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGACCGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCCCTAGCCAACGTTTATTACAGGCGCGCTGGTGCTCCTGGCCTCACTGGGGGTGCATCGTAAGGGGGt                                                                                                                               >  2:325627/1‑246 (MQ=255)                                                                                                                gCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCGGGCCTCACTGTGGGTGCATCGTAAGTGGGt                                                                                                                               <  1:325627/246‑1 (MQ=255)                                                                                                                                tGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATc                                                                                                                                                               <  1:229644/198‑1 (MQ=255)                                                                                                                                tGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAAc                                                                                                          >  1:368585/1‑251 (MQ=255)                                                                                                                                tGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGGAAACCTGGCTATCTGGCACCAACAGCCGCCCTGGCCAACGTTTATTACAGGCGCGCTGGTGaac                                                                                                                                                               >  2:229644/1‑196 (MQ=255)                                                                                                                                   ggCATTCTGGTGTTTCTTGGCGAGTTGGCTTCCGGGATTGGCTACTTTATTTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGc                                                                                                       <  2:353806/251‑1 (MQ=255)                                                                                                                                          tGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACCCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAgatgat                                                                                                >  1:124909/1‑251 (MQ=255)                                                                                                                                              gtTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCa                                                                                            <  1:127791/251‑1 (MQ=255)                                                                                                                                                ttCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAgg                                                                                          >  1:146460/1‑251 (MQ=255)                                                                                                                                                  cTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGAATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCg                                                                                        >  1:194388/1‑251 (MQ=255)                                                                                                                                                            tggCTTCCGGGATTGGCTACTTTATGGGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCCCAAACGGCAGGTGATCGCAGGCGTGATTgcgcg                                                                              >  1:83155/1‑251 (MQ=255)                                                                                                                                                                 tCCGGGATTGGCTACTTTATGTGGAACTAAGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGGAAACCTGGATATCTGGAACCAACAGCCGCACTGGCCAACGTTTATTACCGGCGCGCTGGTGAACCGGGCCCCCCTGTCGGCGCATCGTAAGTGGGTCGCTCAGCCCTCTTCACAAAAGGGCGATGGACCCAGGCGTGCTTGCGCGCTg                                                                           >  2:86016/1‑249 (MQ=255)                                                                                                                                                                          gtgtACTTTATGTGGAACTACGGCGCTACGCAGGTGGACGCCGTAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATaaa                                                                 <  1:208098/247‑1 (MQ=255)                                                                                                                                                                                  tATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATAAGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCg                                                         <  1:379485/250‑1 (MQ=255)                                                                                                                                                                                        gAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATACTATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGATTCCGCAACTgg                                                  >  1:342299/1‑251 (MQ=255)                                                                                                                                                                                                  gcgACGCAGGTGGACGCCGGAACGTTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACgccgc                                        >  2:350130/1‑248 (MQ=255)                                                                                                                                                                                                          gttGGACCCCGGAACGCTGGGCATTATGAATAAGAGCCACGTCTCGGCAT‑GGCTGCTGGTAACCCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGGTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTTTGGGTGCATCGTAAGTGGTTCGCGCCGCGCTCTCCACAAACGGCATATGATCGCAGGCGTGATGGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATc                                <  2:176268/249‑1 (MQ=255)                                                                                                                                                                                                                 gCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCAcg                         >  1:242008/1‑251 (MQ=255)                                                                                                                                                                                                                            gggCAGTATGAATAATATGCACGTTTCGGCAG‑GGCTGCTGGGAAACCTGGCTATCGGGCACCAACAGCCGCACTGGCCAACGTTTATGACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGGCGCTCCGCGCTCTTCACAAACGGCAGGTGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCGGACGCTGCTCGCCATCGCGCACGGCGGCGTACa               <  2:194388/250‑1 (MQ=255)                                                                                                                                                                                                                                         aataTGCACGTTTCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACaa  <  1:25747/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               GATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAG‑GGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACAA  >  NC_000913/4011558‑4012039

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