Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	81,102	G→A	S286L (TCG→TTG)	leuB ←	3‑isopropylmalate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	81,102	0	G	A	94.6%	99.4 / NA	38	S286L (TCG→TTG)	leuB	3‑isopropylmalate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (21/14);  minor base C (2/0);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTAACGA  >  NC_000913/80861‑81336                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           ccATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGccgcccg                                                                                                                                                                                                                                   >  1:310664/1‑251 (MQ=255) ccATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCCATATATGGTGCCCAGccccccc                                                                                                                                                                                                                                   >  2:487476/1‑250 (MQ=255)            ccccTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGACACGGGCTACATCCCCGGTGAGAATACCTTCTTCTAATGCGCGGGTAATGGCGCGGTCAAAGGCGCAAGACGCATTCTCGGCATTCAGGCTGTAACGCAGC‑AGCACGGCCAGCGCAAGGGGTTGTGCAATCGTGTTGGCGATGTTTTTGCCTGGGATATCTGGGGCCAAGACGCACGCCGGGTCAtcc                                                                                                                                                                                                                        >  2:53975/1‑249 (MQ=255)                ttCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATACCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCCTCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGAATTGTGCAATCGGGTTTGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTc                                                                                                                                                                                                                    >  1:20086/1‑251 (MQ=255)                 tCTGCTACATAGCGGGCAATGATATCGCCCTTTTCATCGTTACTAACGGCGGCAGCGCCACGGGCTAAATCCCGGGTGCGAATGCCTTCTTCTAAGGCGCGGTTAATGGCGCGTTAAATGGCGCAAGCCGCATCAGCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGt                                                                                                                                                                                                                     <  1:210989/249‑1 (MQ=255)                    gCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAACCCCCGGTGCGAATGCCTTCTTCTAAGGCGCGTTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTcc                                                                                                                                                                                                                   <  1:412626/248‑1 (MQ=255)                              gggCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCg                                                                                                                                                                                                      >  1:83900/1‑251 (MQ=255)                                          cGCCAATTTCAGCGGTACTACCGGCGGCAGCGCCCCGGGCTAAATCCCCGGTGCGAAGGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCg                                                                                                                                                                                          <  1:347954/251‑1 (MQ=255)                                           gCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCgg                                                                                                                                                                                         >  1:300124/1‑251 (MQ=255)                                           ccccATTTCATGGGTACTAACGGGGGCAGCGCCACGGGCTAAATCCCGGGGGCGAATGCTTTCTTCTAAGGCGCGGTAAAGGGCGCGTACAATGGCGAAAGCCGCACCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCGGCGATATCTGGTGCCAGGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCgg                                                                                                                                                                                         <  2:107667/250‑1 (MQ=255)                                                                   ggcAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGcc                                                                                                                                                                   >  2:133644/1‑249 (MQ=255)                                                                           cACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTc                                                                                                                                                                                                                    >  2:429348/1‑192 (MQ=255)                                                                              gggctaatccccgggggcgatggCCTGCTTCTAAGGCGCGGTTAATGGCGCGTCAAAGGGCGAAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGAAGTGCCAGCGAAAGGATTTGTGAAATCGGGTGGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTgcg                                                                                                                                                      <  1:184218/229‑1 (MQ=255)                                                                              gggCTACATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCCTCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAACGGATTTGTGCAATCTGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCACAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGGTCGTTCAGGCTGGCGGGAGGCCACATcccactcggcacagtggtcatttcg                                                                                                                                                      >  2:43061/1‑229 (MQ=255)                                                                              gggCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAAGGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCACCATCGAGCCAGTGATCATTgcg                                                                                                                                                      <  1:51198/251‑1 (MQ=255)                                                                                          cGGTGCGAATGCCTTCTTCCAATGCGAGGTCAATGGCGCGTTCAAAGGCGCAAGCAGCATCATCGGCATCCAGACTGTAACGCAGC‑ACCCGTGCCAGCGAACGGATTTGTGCAATCGGGATGGCGCTGTTTTTGCATGCGATATATGGTGCCAAGTCGCACGCCGGTTCCTACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGGTCATTgcgccatagtcaga                                                                                                                                           >  2:457330/1‑240 (MQ=255)                                                                                               cGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAAt                                                                                                                                     >  1:273164/1‑251 (MQ=255)                                                                                                                  cgcgGTTAGTGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCCTCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACATTGCTCGTTCCGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGCCAGAATGTGGCCAAACAGGTTgggg                                                                                                                  >  2:496693/1‑249 (MQ=255)                                                                                                                  cgcgGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAACCAGGTTgg                                                                                                                    >  1:437518/1‑249 (MQ=255)                                                                                                                     ggTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCAATGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAgaac                                                                                                               >  1:483972/1‑248 (MQ=255)                                                                                                                         actGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGAAGTGCCAGCGAAAGGATATGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGACGCCCGCCGGTTCATACAGTCCACAACCTTGCTCGTTCAGGCTGGCGGAAGGCACCATCACCATCGAGCCAGAGATCATTGCGCACTCGTCAGACAGCATGTCGCCAAAAAGGCTGGGGCAcaga                                                                                                            >  2:103925/3‑249 (MQ=255)                                                                                                                                 gTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTc                                                                                                    >  1:196758/1‑250 (MQ=255)                                                                                                                                     aaTGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCCAACt                                                                                               >  1:325319/1‑251 (MQ=255)                                                                                                                                         gcgaaaggcgctgcatcGGCATCCAGGGTGAAACGCAGGGAG‑AGGGCCGGGGAAAGGATTGGTGAAAGCGGGTTGGCGATTTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGa                                                                                           <  2:300124/238‑1 (MQ=255)                                                                                                                                                    atcatcGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCCGTGATCATTGCGCACTCGTCAGACCGAATGTCGCCAAACAGGTTGGAGCACAGCTGAACGTCAAACTTTGATGGTTCTTTa                                                                                 >  2:346030/1‑250 (MQ=255)                                                                                                                                                      catcGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATc                                                                              <  1:189174/251‑1 (MQ=255)                                                                                                                                                         cGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGc                                                                           <  1:171430/251‑1 (MQ=255)                                                                                                                                                            cATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGc                                                                        <  1:257460/251‑1 (MQ=255)                                                                                                                                                                              gc‑agcagTGCCAGAGCAAGGATTTGTGCAATCGGGATGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGAGGAGTTGTCGATg                                                      >  1:356503/1‑251 (MQ=255)                                                                                                                                                                               c‑agcagTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGt                                                     >  2:18300/1‑251 (MQ=255)                                                                                                                                                                                      ggacAGCGAAAGGATTGGTGCAACCGGGTGGGCGATGTTTTGGCCTGCGATATCTGGGGCCAAGCCGCCCGCAGGTTCAGACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACatat                                               <  1:260160/248‑1 (MQ=255)                                                                                                                                                                                         cAGCGAAAGGATTTGTGTAAGGGGGTTGGCTATGTTTTTGCCTGCGTTATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATgcg                                            <  1:450274/251‑1 (MQ=255)                                                                                                                                                                                                     ttGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGCGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACAt                                >  2:421818/1‑251 (MQ=255)                                                                                                                                                                                                          cAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGCTGTACATATGCGCCAGTTCGACATCCggg                           >  2:106699/1‑251 (MQ=255)                                                                                                                                                                                                                   gggCGATTTTTTGGCCTGCGATTTCTGTTGCCAAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCTTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTg                  <  1:276702/250‑1 (MQ=255)                                                                                                                                                                                                                    ggCGATGTTTTTGCCTGCGATATCTGGTGCCAAGACGCCCGCCGGTTCATACAGTCCAAAACATTGCTCTTTCGGGCTGGCGAAAGGCAACATCCACATAGCGACCGTGATCATTGCGCACTCGTCAGACAGAATTTCGCCAAACAGATTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCGACATGGTGGCGTAGTCGGTGAACATATGCCATAGTTCAACATCCGGGtctacagt                   >  2:462928/1‑242 (MQ=255)                                                                                                                                                                                                                                gcatGCGATACCTGGTGCCAGGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTaa     <  1:487476/248‑1 (MQ=255)                                                                                                                                                                                                                                    gCGATAGCTGGTGCCACGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTAACGa  <  1:130362/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           CCATGATTACACCCCTTCTGCTACATAGCGGGCAATGATATCGCCCATTTCATCGGTACTAACGGCGGCAGCGCCACGGGCTAAATCCCCGGTGCGAATGCCTTCTTCTAATGCGCGGTTAATGGCGCGTTCAATGGCGCAAGCCGCATCATCGGCATCCAGGCTGTAACGCAGC‑AGCAGTGCCAGCGAAAGGATTTGTGCAATCGGGTTGGCGATGTTTTTGCCTGCGATATCTGGTGCCGAGCCGCCCGCCGGTTCATACAGTCCAAAACCTTGCTCGTTCAGGCTGGCGGAAGGCAACATCCCCATCGAGCCAGTGATCATTGCGCACTCGTCAGACAGAATGTCGCCAAACAGGTTGGAGCACAGCAGAACGTCAAACTGTGATGGATCTTTAATCAGCTGCATGGTGGCGTTGTCGATGTACATATGCGCCAGTTCGACATCCGGGTATTCCGTGGCGATCTCGTTAACGA  >  NC_000913/80861‑81336

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