Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,312,434	C→T	G153S (GGC→AGC)	yciB ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,434	0	C	T	98.0%	169.8 / NA	50	G153S (GGC→AGC)	yciB	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/26);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCCAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCATAAATGACAGTAACCTTCCATTTAATAAACTC  >  NC_000913/1312189‑1312677                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      gCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGGCAGGCTAAAGa                                                                                                                                                                                                                                                  <  2:536897/251‑1 (MQ=255)          aaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGGGGCGGTAGATAAAGATACCGCTTAACCATGTGAAGAATAAAGGAAAGGCGGTCCGGCTAAAGACTTTaaaa                                                                                                                                                                                                                                          >  2:513516/1‑249 (MQ=255)              aCGTAAAACAACATCACCATGAGGGACGTTATGTGTTGTAGACATGGTAAAAAAGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGTTTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCATGCGGCATGTGGCGGGAGATATAGATACCGCTTAACAAAGTGAAGATTAAGGTGAGGGCGGGCAGGCTAAAGCATTTAAAggtggc                                                                                                                                                                                                                                     >  2:71284/1‑246 (MQ=255)                   aaaCAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCaa                                                                                                                                                                                                                                 >  1:122994/1‑250 (MQ=255)                        aGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATAttt                                                                                                                                                                                                                           >  2:445222/1‑251 (MQ=255)                                  gAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTTCGGCAg                                                                                                                                                                                                                   <  1:512066/249‑1 (MQ=255)                                     ggACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGcc                                                                                                                                                                                                                 >  1:333399/1‑248 (MQ=255)                                        cGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCGTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTGGCGGCAGCCAGAATg                                                                                                                                                                                                           <  1:75264/251‑1 (MQ=255)                                          ttATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCg                                                                                                                                                                                                         >  1:550806/1‑251 (MQ=255)                                           tATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGGCAGGCTAAAGACTTTAAAGGTGACCCAAATATTTTGCGGCAACCAGAATGCg                                                                                                                                                                                                         >  1:454952/1‑250 (MQ=255)                                                   gtAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATg                                                                                                                                                                                                <  1:445222/251‑1 (MQ=255)                                                   gtAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCAAAATATTTTGCGGCAGCCAGAATGCGATGTAGAt                                                                                                                                                                                                 <  1:39227/250‑1 (MQ=255)                                                    tAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAAATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGGAGATGt                                                                                                                                                                                               >  2:296810/1‑251 (MQ=255)                                                     aGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTTAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCgatgttgatgtt                                                                                                                                                                                              >  2:438468/1‑251 (MQ=255)                                                          tGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGGAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGcc                                                                                                                                                                                          >  2:64528/1‑250 (MQ=255)                                                           ggTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAg                                                                                                                                                                                        <  2:517883/251‑1 (MQ=255)                                                              aaaaCCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACACTGTGAAGATTAAAGGCAGGGCGGTCAGGCTAAAGACTTTAAAGTTTTCCCAAATATTTTTCCGCAGCCCGAATGCGATGTAAATTTTTTCCAGaac                                                                                                                                                                                     >  2:270874/1‑249 (MQ=255)                                                              aaaaCCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAc‑c                                                                                                                                                                                    <  2:230681/251‑1 (MQ=255)                                                                    gACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCGTAACAATGTGAAGATTAAGGTAAGGGCGGGCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTGGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAg                                                                                                                                                                              <  1:350558/251‑1 (MQ=255)                                                                              gAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGaaa                                                                                                                                                                    <  1:87543/251‑1 (MQ=255)                                                                                              aTGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGAt                                                                                                                                                   <  1:58729/251‑1 (MQ=255)                                                                                                      aTGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTc                                                                                                                                           <  2:180259/251‑1 (MQ=255)                                                                                                      aTGACTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCtt                                                                                                                                            <  2:211195/250‑1 (MQ=255)                                                                                                            ttCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCa                                                                                                                                      <  2:167756/250‑1 (MQ=255)                                                                                                                         cAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCa                                                                                                                        <  1:487696/251‑1 (MQ=255)                                                                                                                          agagGAGATACACCGACAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGCGGCGGTCGATATAGAAACCGCTTAACAATGAGAAGAGTGAGGGAAGGGCGGTCAGGCTAAAGACTTTAAAATTGACCCAAATAATTTGCAGCAGCCAGAATGAGATGTCGATGTTTGCCAGACACG‑AAAGGATAAAGAAAACAGACCAAGCCCAAATT‑AGCTTCCACCATACCTGTTGCGGccg                                                                                                                       >  2:466300/1‑249 (MQ=255)                                                                                                                                       cgcCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTa                                                                                                          <  1:90770/251‑1 (MQ=255)                                                                                                                                         ccAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTAcc                                                                                                        <  1:433132/251‑1 (MQ=255)                                                                                                                                                                          ttATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGCAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGAACATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCCTCAc                                                                       >  1:377842/1‑251 (MQ=255)                                                                                                                                                                           tATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGGAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATGGAGCGGCAGCCAGAATGCGATGTAGGTGGTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCTAGTCCAGATTCAGCTTCGACCAGACCGGTTGCGGCAGCGTGATTTGCTTACCCAGGATCCGCTGAATTAGTGGcaatttgatcac                                                                       <  2:375920/250‑12 (MQ=255)                                                                                                                                                                             tCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTTAAGATTAAGGGAAGGGCGGTAAGGCCAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGAAGTTTGCCAGAA‑CGCAAAGGATAAAGAAAACAA‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAACGTGCGGTCTTTACCCCGCATCCGCTGAATTAGCGGCGTTTTCATCAcccc                                                                    >  2:122556/1‑250 (MQ=255)                                                                                                                                                                                    gCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTg                                                             <  1:102833/251‑1 (MQ=255)                                                                                                                                                                                      ggCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCCTTGGCTGAc                                                           >  2:352309/1‑251 (MQ=255)                                                                                                                                                                                        cATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTTGCttaat                                                          >  2:300519/1‑246 (MQ=255)                                                                                                                                                                                           gtgGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAAc                                                       <  2:214131/250‑1 (MQ=255)                                                                                                                                                                                                  tagatatagataCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGATcc                                               >  1:327886/1‑251 (MQ=255)                                                                                                                                                                                                     atatagataCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGc                                            >  1:478114/1‑251 (MQ=255)                                                                                                                                                                                                        tagataCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGCCCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCaa                                          >  2:380681/1‑250 (MQ=255)                                                                                                                                                                                                         agataCCGCTTCACAATGTGGAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACACAAATATTTTGCGGCAGCCAGAATGCGATGTAGCTGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCCGCTTCGACCATACCGGATGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCCTTGGCTGACTAACAGGGCTCCCGCAAAt                                        >  2:80186/1‑251 (MQ=255)                                                                                                                                                                                                          gataCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGTCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTGTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATa                                       <  2:282787/251‑1 (MQ=255)                                                                                                                                                                                                          gataCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAAt                                        >  2:184518/1‑250 (MQ=255)                                                                                                                                                                                                          gataCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATa                                       <  1:296810/251‑1 (MQ=255)                                                                                                                                                                                                             aCCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATagag                                    <  1:530975/251‑1 (MQ=255)                                                                                                                                                                                                                gCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCAt                                 <  1:423772/251‑1 (MQ=255)                                                                                                                                                                                                                       aaTGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAg                                                                                                                                                                              >  1:350513/1‑104 (MQ=255)                                                                                                                                                                                                                       aaTGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAg                                                                                                                                                                              <  2:350513/104‑1 (MQ=255)                                                                                                                                                                                                                        aTGTGAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCATAAATGACa                         <  2:522886/251‑1 (MQ=255)                                                                                                                                                                                                                            gAAGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCATAAATGACAGTa                      >  1:403293/1‑250 (MQ=255)                                                                                                                                                                                                                              aGATTAAGGTAAGGGCGGTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CTCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCATAAATGACAGTAAcc                   <  1:80186/251‑1 (MQ=255)                                                                                                                                                                                                                                               gTCAGGCTAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCATAAATGACAGTAACCTTCCATTTAATAAACtc  >  2:568382/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      GCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCTTCCTGCGGCATGTGGCGGTAGATATAGATACCGCTTAACAATGTGAAGATTAAGGTAAGGGCGGTCAGGCCAAAGACTTTAAAGTTGACCCAAATATTTTGCGGCAGCCAGAATGCGATGTAGATGTTTGCCAGAC‑CGCAAAGGATAAAGAAAACAG‑CCCAGGCCAGATTCAGCTTCGACCATACCGGTTGCGGCAGCGTGAGTTCTTTACCCAGCATCCGCTGAATTAGCGGCTTTTTCATCACCCATTGGCTGACTAACAGGGCTCCCGCAAATAGAGCATAAATGACAGTAACCTTCCATTTAATAAACTC  >  NC_000913/1312189‑1312677

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