Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,180,039	A→T	6.3%	Q99L (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,039	0	A	T	6.3%	223.1 / 4.3	63	Q99L (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.98e-01

GCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTC  >  NC_000913/3179891‑3180186                                                                                                                                                     |                                                                                                                                                    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                                                            <  1:619755/147‑1 (MQ=255)                                                                                            cGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGtt                                                            >  2:619755/1‑147 (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                    <  4:440175/149‑1 (MQ=255)                                                                                                                                  ggTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGc                    <  1:638275/149‑1 (MQ=255)                                                                                                                                   gTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                   >  3:304711/1‑149 (MQ=255)                                                                                                                                       aGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagc               >  2:692112/1‑149 (MQ=255)                                                                                                                                         gTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagcag             >  3:305452/1‑149 (MQ=255)                                                                                                                                         gTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagcag             >  4:250404/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)                                                                                                                                                    cagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtctc  <  4:304711/149‑1 (MQ=255)                                                                                                                                                    cagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtctc  <  4:488977/149‑1 (MQ=255)                                                                                                                                                    cagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtctc  <  3:118252/149‑1 (MQ=255)                                                                                                                                                    cagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtctc  <  3:590986/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTCTC  >  NC_000913/3179891‑3180186

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