Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,500,389	(T)8→7	coding (61/981 nt)	ydcK ←	putative enzyme YdcK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,500,382	0	T	.	98.4%	268.4 / ‑1.9	64	coding (68/981 nt)	ydcK	putative enzyme YdcK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (30/33);  minor base C (1/0);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACATTTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGA  >  NC_000913/1500244‑1500514                                                                                                                                           |                                                                                                                                     aaCGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtt                                                                                                                                     >  2:451823/1‑140 (MQ=255)          gCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACAtttttttagtg                                                                                                                            <  2:154665/140‑5 (MQ=255)            aTTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGcc                                                                                                                         <  2:1104903/140‑1 (MQ=255)                tGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCt                                                                                                                     >  2:251556/1‑140 (MQ=255)                tGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCt                                                                                                                     <  1:240858/140‑1 (MQ=255)                  gTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTc                                                                                                                   >  1:1030211/1‑140 (MQ=255)                  gTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTc                                                                                                                   >  1:570256/1‑140 (MQ=255)                     gTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTg                                                                                                                >  2:252668/1‑140 (MQ=255)                       aaaTCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGAt                                                                                                              <  1:484098/140‑1 (MQ=255)                               cAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAg                                                                                                      <  2:266837/140‑1 (MQ=255)                                  tttCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCt                                                                                                   <  1:291739/140‑1 (MQ=255)                                          ttGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGttc                                                                                           <  1:996157/140‑1 (MQ=255)                                            gtgCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGttctt                                                                                         >  1:919605/1‑140 (MQ=255)                                                     ccGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGAc                                                                                <  1:410628/140‑1 (MQ=255)                                                     ccGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGAc                                                                                <  1:572168/140‑1 (MQ=255)                                                       gTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGa                                                                              >  2:597374/1‑140 (MQ=255)                                                       gTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGa                                                                              >  1:690471/1‑140 (MQ=255)                                                       gTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGa                                                                              >  2:1107239/1‑140 (MQ=255)                                                        tttCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGat                                                                             >  2:803011/1‑140 (MQ=255)                                                        tttCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGat                                                                             >  2:216816/1‑140 (MQ=255)                                                         ttCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGata                                                                            >  2:1170807/1‑140 (MQ=255)                                                          tCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGatat                                                                           <  2:1169994/140‑1 (MQ=255)                                                           cgcgATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACACTTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGatat                                                                           >  2:1248717/1‑140 (MQ=255)                                                           cgcgATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATAtt                                                                          <  2:542226/140‑1 (MQ=255)                                                                   atcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATa                                                                  <  1:93285/140‑1 (MQ=255)                                                                   atcCAGCCGCCAGCGGTTCCTGCTATCACATCGTCAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTGGCGCATa                                                                  <  1:212780/140‑1 (MQ=255)                                                                    tcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAg                                                                 <  2:919605/140‑1 (MQ=255)                                                                    tcCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAg                                                                 <  2:828643/140‑1 (MQ=255)                                                                             cAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGtt                                                             <  1:857994/135‑1 (MQ=255)                                                                             cAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGtt                                                             >  2:857994/1‑135 (MQ=255)                                                                              aGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCca                                                       >  2:260175/1‑140 (MQ=255)                                                                              aGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCca                                                       >  1:146566/1‑140 (MQ=255)                                                                              aGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCca                                                       >  1:743812/1‑140 (MQ=255)                                                                              aGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCca                                                       >  1:857690/1‑140 (MQ=255)                                                                               gCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCcac                                                      <  1:360837/140‑1 (MQ=255)                                                                               gCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCcac                                                      <  1:451823/140‑1 (MQ=255)                                                                                 ggTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAg                                                    <  2:670140/140‑1 (MQ=255)                                                                                  gTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAgg                                                   <  1:216816/140‑1 (MQ=255)                                                                                     ccTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGa                                                <  1:338847/140‑1 (MQ=255)                                                                                      cTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGAc                                               <  1:840089/140‑1 (MQ=255)                                                                                       tGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACt                                              <  1:1044454/140‑1 (MQ=255)                                                                                             cacaTCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACtat                                        >  2:906816/1‑140 (MQ=255)                                                                                             cacaTCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACtat                                        >  2:1049554/1‑140 (MQ=255)                                                                                             cacaTCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACtat                                        >  2:1038765/1‑140 (MQ=255)                                                                                             cacaTCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACtat                                        >  2:54682/1‑140 (MQ=255)                                                                                                  cGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCg                                   >  2:326671/1‑140 (MQ=255)                                                                                                     tAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGa                                                <  2:816469/124‑1 (MQ=255)                                                                                                     tAAAATCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGa                                                >  1:816469/1‑124 (MQ=255)                                                                                                      aaaaTCGCTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGc                               >  1:1187926/1‑140 (MQ=255)                                                                                                             cTAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCt                        <  1:929921/140‑1 (MQ=255)                                                                                                              tAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTg                       <  2:857690/140‑1 (MQ=255)                                                                                                              tAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTg                       <  1:1038765/140‑1 (MQ=255)                                                                                                              tAATGGCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTg                       >  1:381562/1‑140 (MQ=255)                                                                                                                   gCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGct                  <  2:372122/140‑1 (MQ=255)                                                                                                                   gCGATAATCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGct                  <  1:803011/140‑1 (MQ=255)                                                                                                                         aTCTGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCt            <  1:122910/140‑1 (MQ=255)                                                                                                                            tGCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATg         <  1:54682/140‑1 (MQ=255)                                                                                                                             gCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGc        >  2:286155/1‑140 (MQ=255)                                                                                                                             gCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGc        >  2:1187986/1‑140 (MQ=255)                                                                                                                             gCCGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGc        >  2:834337/1‑140 (MQ=255)                                                                                                                              ccGTAATAACACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCt       >  2:929021/1‑140 (MQ=255)                                                                                                                                   ataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGct                  >  2:187292/1‑124 (MQ=255)                                                                                                                                   ataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGct                  <  1:187292/124‑1 (MQ=255)                                                                                                                                   ataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCgaga  >  1:904704/1‑140 (MQ=255)                                                                                                                                   ataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCgaga  <  2:91852/140‑1 (MQ=255)                                                                                                                                   ataaCACA‑TTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCgaga  <  1:1170807/140‑1 (MQ=255)                                                                                                                                           |                                                                                                                                     AACGCAATAGCATTCTGGTCGTAAATCCAGCAATTTCCTTCTTGTGCTAACACCGTTTCGCGATCGATCCAGCCGCCAGCGGTTCCTGCTATCACATCGTTAAAATCGCTAATGGCGATAATCTGCCGTAATAACACATTTTTTTTAGTGCCATCTTCCTGATAACTAAAGGCTCGCTGTTCTTCACTTAGACGATATTTGCGCATAGCGTTTTCCCACAGGTGACTTACTATAACCGTAGCAAATTCTGCGGCTCTGGCTATGCTCGAGA  >  NC_000913/1500244‑1500514

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