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	97.7%	151.5 / NA	44	P243S (CCG→TCG)	bioP	biotin transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/18);  minor base G (1/0);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACA  >  NC_000913/4011555‑4012038                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             ccTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCgg                                                                                                                                                                                                                                           >  2:280472/1‑251 (MQ=255)  cTGATGGAAACTCGCCCGATGCCACAGCCTAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTGCGGc                                                                                                                                                                                                                                          >  2:133289/1‑251 (MQ=255)    gATGGAAACTCGCCCGCTGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGAACGTTTCGGc                                                                                                                                                                                                                                          >  2:304336/1‑249 (MQ=255)      tGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTCCTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGACTAATATGCACGTTTCGGCAggg                                                                                                                                                                                                                                      >  2:280961/1‑251 (MQ=255)      tGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAggg                                                                                                                                                                                                                                      >  1:106387/1‑251 (MQ=255)              cGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGt                                                                                                                                                                                                                              >  1:177932/1‑251 (MQ=255)                         ccaGCATAACGCCTTTGCGTGGTTCTATCTTGTCGCGTTTCTTGTGGCAGTGATTGCATGGTTCTTGCTGGGCAATGCGCAGAAAATGCCGCAAACCACGCTGCTATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTCATTTATGTGGAACTACGGCGCGACGCCTGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTTGCt                                                                                                                                                                                                                    >  2:187483/2‑250 (MQ=255)                                    ccTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCaa                                                                                                                                                                                                        <  2:263808/251‑1 (MQ=255)                                       ttGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACa                                                                                                                                                                                                      >  2:284929/1‑250 (MQ=255)                                          cGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTCCTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGACCTCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCCGTTCACCCTGGCTATCTGGCACCAACAGCCg                                                                                                                                                                                                  >  2:140646/1‑251 (MQ=255)                                           gTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGc                                                                                                                                                                                                 >  1:196407/1‑251 (MQ=255)                                             ggTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCa                                                                                                                                                                                                >  1:150879/1‑250 (MQ=255)                                              gTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCgg                                                                                                                                                                                                                                           <  1:25088/206‑1 (MQ=255)                                              gTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCgg                                                                                                                                                                                                                                           >  2:25088/1‑206 (MQ=255)                                                       tGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGt                                                                                                                                                                                     >  1:146538/1‑251 (MQ=255)                                                             gTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTttatta                                                                                                                                                                               <  2:301613/251‑1 (MQ=255)                                                                   ggtggCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGgcg                                                                                                                                                                         >  2:97593/1‑251 (MQ=255)                                                                                         ttGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCAc                                                                                                                                                    <  2:263302/250‑1 (MQ=255)                                                                                           gCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACtg                                                                                                                                                  <  1:338450/250‑1 (MQ=255)                                                                                                                gCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCTACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCg                                                                                                                            <  1:164306/251‑1 (MQ=255)                                                                                                                  cGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGGCAGTGGGTCGCt                                                                                                                          >  2:118191/1‑251 (MQ=255)                                                                                                                  cGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTGCAGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGTACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCGGGCACCAACAGCCGCACTGGCCACCGTTTATTACAGGCGCGCTGGGGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCt                                                                                                                          <  1:187483/251‑1 (MQ=255)                                                                                                                     aaaCCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCcg                                                                                                                       >  2:263642/1‑251 (MQ=255)                                                                                                                        ccACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCgcgc                                                                                                                    >  2:358911/1‑251 (MQ=255)                                                                                                                                  aTGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCCGGGCTGCTGGTAAACCTGGCTATCTGGCAACAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACaa                                                                                                           >  1:361419/1‑250 (MQ=255)                                                                                                                                   gggggCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAAc                                                                                                         <  2:272420/250‑1 (MQ=255)                                                                                                                                    ggggCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAAc                                                                                                         <  1:263642/250‑1 (MQ=255)                                                                                                                                         aTTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGCTTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGc                                                                                                      >  2:235929/1‑248 (MQ=255)                                                                                                                                                tgtTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACACTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGc                                                                                            >  2:177348/1‑251 (MQ=255)                                                                                                                                                 ggttCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGc                                                                                            >  1:11719/3‑250 (MQ=255)                                                                                                                                                                 gCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAAGAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATGACAGGCGCGCTGGTGATCCTGTCCGCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCt                                                                           <  1:200336/251‑1 (MQ=255)                                                                                                                                                                     ccGGGATTGGCTACTTTATGTGGAACTACTGCGCGACGCAGGTGGACGCCGGAACGATGGGCATTATGAATCATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAAGCGCGCTGGTGATACTGGCCTAACTGTGGGGGCATCGAAAGTGGGTCGCCCCGCGCTCTTCACAACCGGCAGATGATCGCAGGAGTGATTGCGCGCTGAGc                                                                       >  2:210675/1‑251 (MQ=255)                                                                                                                                                                     ccGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAg                                                                        <  1:118191/250‑1 (MQ=255)                                                                                                                                                                       gggATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGa                                                                     <  1:115531/251‑1 (MQ=255)                                                                                                                                                                          aTTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATa                                                                  <  1:283811/251‑1 (MQ=255)                                                                                                                                                                            tGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATaa                                                                 <  1:304336/250‑1 (MQ=255)                                                                                                                                                                                    ttATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCg                                                        >  1:73805/1‑251 (MQ=255)                                                                                                                                                                                      aTGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTGCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGt                                                       <  1:185423/250‑1 (MQ=255)                                                                                                                                                                                            aaCTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGc                                                <  2:146538/251‑1 (MQ=255)                                                                                                                                                                                             ccTACGGCGAGACGCAGGTGGACGGCGGAGCGCTGGGCATTATGAATAATATGCCCGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCTCCAACAGCGGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCGTCACAGACGGCAGATGATCGCAGGCGTGATTGTGCGCGGTGCGAAGAAACGCTGCCGCAACTGGCt                                               <  2:158937/250‑1 (MQ=255)                                                                                                                                                                                                         ggCAGGTGGACGCCGGAACGCTGGGCATTATGAATAAGAGGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGGGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGc                                   <  2:43804/250‑1 (MQ=255)                                                                                                                                                                                                             ggTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATc                               >  2:136640/1‑251 (MQ=255)                                                                                                                                                                                                                  aCGCCGGAACGCTGGGCATTATGAATAATATGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCgcgcc                          >  2:329837/1‑250 (MQ=255)                                                                                                                                                                                                                                           taataTGCACGTTTCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCcaca  >  1:178863/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             CCTGATGGAAACTCGCCCGATGCCACAGCATAACGCCTTTGCGTGGTTCTATCTTGGCGCGTTTCTGGTGGCAGTGATTGCATGGTTCTTGCTGGGAAATGCGCAGAAAATGCCGCAAACCACGCTGCAATGGGGCATTCTGGTGTTTCTTGGCGTGGTGGCTTCCGGGATTGGCTACTTTATGTGGAACTACGGCGCGACGCAGGTGGACGCCGGAACGCTGGGCATTATGAATAATATGCACGTTCCGGCAGGGCTGCTGGTAAACCTGGCTATCTGGCACCAACAGCCGCACTGGCCAACGTTTATTACAGGCGCGCTGGTGATCCTGGCCTCACTGTGGGTGCATCGTAAGTGGGTCGCTCCGCGCTCTTCACAAACGGCAGATGATCGCAGGCGTGATTGCGCGCTGAGCGAATAAACGCTTCCGTAACTGGCTGACGCTGCTCGCCATCGCGCACGGCGGCGTACAGTCGGCTCCACA  >  NC_000913/4011555‑4012038

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