Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,036	A→T	6.3%	Q98L (CAG→CTG)	ygiB →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,180,036	0	A	T	6.3%	223.1 / 4.6	63	Q98L (CAG→CTG)	ygiB	DUF1190 family protein
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTC  >  NC_000913/3179891‑3180184                                                                                                                                                  |                                                                                                                                                     gCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCcagca                                                                                                                                                   <  4:344302/149‑1 (MQ=255)       gacTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCAcc                                                                                                                                             >  2:652546/1‑149 (MQ=255)       gacTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCAcc                                                                                                                                             >  2:638275/1‑149 (MQ=255)       gacTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCAcc                                                                                                                                             >  2:440306/1‑149 (MQ=255)               aGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAgg                                                                                                                                     <  1:300982/149‑1 (MQ=255)               aGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAgg                                                                                                                                     <  1:608843/149‑1 (MQ=255)                gCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGc                                                                                                                                    >  2:36926/1‑149 (MQ=255)                  tGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtg                                                                                                                                  >  4:166598/1‑149 (MQ=255)                   gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTCCTGCGCCGAAATACGCGACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                 <  1:47050/149‑1 (MQ=255)                   gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                 >  2:599750/1‑149 (MQ=255)                   gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                 >  4:384909/1‑149 (MQ=255)                     aaaCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggca                                                                                                                               <  3:384909/149‑1 (MQ=255)                     aaaCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggca                                                                                                                               <  4:615558/149‑1 (MQ=255)                       aCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggcatg                                                                                                                             >  4:70671/1‑149 (MQ=255)                         ccAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggcatggc                                                                                                                           >  4:24082/1‑149 (MQ=255)                         ccAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggcatggc                                                                                                                           >  3:332568/1‑149 (MQ=255)                             gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCa                                                                                                                       <  2:267516/149‑1 (MQ=255)                             gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCa                                                                                                                       >  1:543691/1‑149 (MQ=255)                                  aGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGaaaa                                                                                                                  <  1:36926/149‑1 (MQ=255)                                     gcCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCa                                                                                                               >  3:63805/1‑149 (MQ=255)                                      ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGTATGGCACCAGAAAACCAg                                                                                                              <  4:393501/149‑1 (MQ=255)                                      ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAg                                                                                                              <  1:440306/149‑1 (MQ=255)                                      ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAg                                                                                                              <  2:546542/149‑1 (MQ=255)                                        gAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGgc                                                                                                            <  3:166598/149‑1 (MQ=255)                                             tACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAgg                                                                                                       <  4:332568/149‑1 (MQ=255)                                                caccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGccc                                                                                                    >  1:81537/1‑149 (MQ=255)                                                 accGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCa                                                                                                   <  3:643847/149‑1 (MQ=255)                                                 accGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCa                                                                                                   <  1:338819/149‑1 (MQ=255)                                                          aacaaTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAg                                                                                          >  1:247164/1‑149 (MQ=255)                                                             aaTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCgg                                                                                       <  1:691420/149‑1 (MQ=255)                                                               tGCGCTGAAAGAAGCCGAACGTACCGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGa                                                                                     >  3:440175/1‑149 (MQ=255)                                                                    tGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTc                                                                                <  4:160115/149‑1 (MQ=255)                                                                     gAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCt                                                                               <  2:358272/149‑1 (MQ=255)                                                                      aaaGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTg                                                                              <  2:81537/149‑1 (MQ=255)                                                                       aagaagCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGccc                                                                                                                                        <  1:235440/90‑1 (MQ=255)                                                                       aagaagCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGccc                                                                                                                                        >  2:235440/1‑90 (MQ=255)                                                                         gaagCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGAt                                                                           <  2:651545/149‑1 (MQ=255)                                                                           agCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGc                                                                         >  4:301760/1‑149 (MQ=255)                                                                                            cGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGtt                                                          >  2:619755/1‑147 (MQ=255)                                                                                            cGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGtt                                                          <  1:619755/147‑1 (MQ=255)                                                                                             gAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACa                                                       <  1:652546/149‑1 (MQ=255)                                                                                                      cACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGc                                              <  4:543045/149‑1 (MQ=255)                                                                                                                gACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATggg                                    >  3:384166/1‑149 (MQ=255)                                                                                                                 aCTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggc                                   <  1:599750/149‑1 (MQ=255)                                                                                                                      gtTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggcggcgg                              >  4:590986/1‑149 (MQ=255)                                                                                                                       ttGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggcggcggc                             >  3:488977/1‑149 (MQ=255)                                                                                                                        tGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGa                                                                            <  4:102927/101‑1 (MQ=255)                                                                                                                        tGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGa                                                                            >  3:102927/1‑101 (MQ=255)                                                                                                                          cTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACa                                                       <  1:626146/120‑1 (MQ=255)                                                                                                                          cTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACa                                                       >  2:626146/1‑120 (MQ=255)                                                                                                                                ttGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGAttt                    <  3:69/149‑1 (MQ=255)                                                                                                                                  ggTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGc                  <  1:638275/149‑1 (MQ=255)                                                                                                                                  ggTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGc                  <  4:440175/149‑1 (MQ=255)                                                                                                                                   gTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                 >  3:304711/1‑149 (MQ=255)                                                                                                                                       aGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagc             >  2:692112/1‑149 (MQ=255)                                                                                                                                         gTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagcag           >  4:250404/1‑149 (MQ=255)                                                                                                                                         gTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagcag           >  3:305452/1‑149 (MQ=255)                                                                                                                                          ctggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                 <  4:12945/139‑1 (MQ=255)                                                                                                                                          ctggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                 >  3:12945/4‑142 (MQ=255)                                                                                                                                            ggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGa                                        <  3:554169/116‑1 (MQ=255)                                                                                                                                            ggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGa                                        >  4:554169/2‑117 (MQ=255)                                                                                                                                             gTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGc       >  3:409358/1‑149 (MQ=255)                                                                                                                                             gTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGc       >  2:464403/1‑149 (MQ=255)                                                                                                                                              tGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCt      <  1:464403/149‑1 (MQ=255)                                                                                                                                                  agcagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtc  >  3:538940/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTC  >  NC_000913/3179891‑3180184

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