Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,992,083	G→C	7.4%	noncoding (35/76 nt)	glyW ←	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,992,083	0	G	C	7.4%	233.0 / 9.6	68	noncoding (35/76 nt)	glyW	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (27/36);  new base C (2/3);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.32e-01

AACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTT  >  NC_000913/1991935‑1992228                                                                                                                                                     |                                                                                                                                                  aaCCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                                                                                                                   >  2:216743/1‑149 (MQ=255)        cTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                                                                                                                            <  3:273603/149‑1 (MQ=255)         tAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCg                                                                                                                                           <  3:488655/149‑1 (MQ=255)           gACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTg                                                                                                                                         <  1:156975/149‑1 (MQ=255)           gACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCACGACCTTGCCAAGGTCggg                                                                                                                                         <  1:663545/149‑3 (MQ=255)                aTTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAAATGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTa                                                                                                                                    >  2:233581/1‑149 (MQ=255)                      aaTCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCTGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCACGACCTTGCCAAGGTCgg                                                                                                                                          >  2:575215/1‑136 (MQ=255)                      aaTCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCTGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCACGACCTTGCCAAGGTCgg                                                                                                                                          <  1:575215/137‑2 (MQ=255)                       aTCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTg                                                                                                                             <  2:578547/149‑1 (MQ=255)                        tCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            <  2:633336/149‑1 (MQ=255)                        tCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            <  3:529381/149‑1 (MQ=255)                        tCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            >  2:175115/1‑149 (MQ=255)                        tCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            <  1:236119/149‑1 (MQ=255)                             tACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTTTACCAACTGAGCTAt                                                                                                                       <  2:688792/149‑1 (MQ=255)                               cATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTc                                                                                                                     <  1:12256/149‑1 (MQ=255)                                     cGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCAt                                                                                                               >  1:671286/1‑149 (MQ=255)                                     cGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCAt                                                                                                               <  1:175115/149‑1 (MQ=255)                                        tttGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtca                                                                                                            >  3:621112/1‑149 (MQ=255)                                              aaCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                                                                                                                         >  3:193451/1‑112 (MQ=255)                                              aaCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                                                                                                                         <  4:193451/114‑3 (MQ=255)                                                  cgcCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATc                                                                                                  <  4:177385/149‑1 (MQ=255)                                                   gcCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCa                                                                                                 <  4:549771/149‑1 (MQ=255)                                                     cAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGt                                                                                               >  2:609760/1‑149 (MQ=255)                                                      aaaTTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGtt                                                                                              >  1:296313/1‑149 (MQ=255)                                                      aaaTTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGtt                                                                                              >  4:178595/1‑149 (MQ=255)                                                         ttcttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAAt                                                                                           <  1:543158/149‑1 (MQ=255)                                                         ttcttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAAt                                                                                           >  2:49091/1‑149 (MQ=255)                                                         ttcttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAAt                                                                                           <  2:296313/149‑1 (MQ=255)                                                           cttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCa                                                                                         <  1:216743/149‑1 (MQ=255)                                                            ttcAGGCCTTTCAGCCAGACATGCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCAc                                                                                        >  1:571868/1‑149 (MQ=255)                                                            ttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCAc                                                                                        <  2:154087/149‑1 (MQ=255)                                                                ggCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGa                                                                                    >  2:422336/1‑149 (MQ=255)                                                                ggCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGa                                                                                    <  1:392860/149‑1 (MQ=255)                                                                 gCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAt                                                                                   >  3:652676/1‑149 (MQ=255)                                                                    tttCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAtttt                                                                                >  3:306203/1‑149 (MQ=255)                                                                    tttCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAtttt                                                                                >  3:291131/1‑149 (MQ=255)                                                                      tCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTttat                                                                              <  4:277608/149‑1 (MQ=255)                                                                         gCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTtattat                                                                           <  1:609760/149‑1 (MQ=255)                                                                                 tCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGc                                                                   <  2:671286/149‑1 (MQ=255)                                                                                           gCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgcc                                                         >  3:218776/1‑149 (MQ=255)                                                                                             cGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgccgc                                                       >  2:95439/1‑149 (MQ=255)                                                                                              gATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgccgcc                                                      >  2:470951/1‑149 (MQ=255)                                                                                                    ttttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGa                                                <  3:178595/149‑1 (MQ=255)                                                                                                    ttttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGa                                                >  2:302459/1‑149 (MQ=255)                                                                                                     tttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGAt                                               <  1:233581/149‑1 (MQ=255)                                                                                                          aCTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGtt                                          <  4:306203/149‑1 (MQ=255)                                                                                                           cTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTg                                         <  4:218776/149‑1 (MQ=255)                                                                                                           cTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTg                                         <  4:333002/149‑1 (MQ=255)                                                                                                                 cGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCt                                   <  1:631014/149‑1 (MQ=255)                                                                                                                     aaaCGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACtt                               <  3:627605/149‑1 (MQ=255)                                                                                                                       aCGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCAc                                                                                        >  3:488468/1‑90 (MQ=255)                                                                                                                       aCGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCAc                                                                                        <  4:488468/90‑1 (MQ=255)                                                                                                                          agaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTg                          <  1:49091/149‑1 (MQ=255)                                                                                                                           gaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTgg                         >  2:634789/1‑149 (MQ=255)                                                                                                                              tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGgcgc                      >  4:518765/1‑149 (MQ=255)                                                                                                                               cGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGcgcg                     >  4:490968/1‑149 (MQ=255)                                                                                                                                   cTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGa                 <  2:14929/149‑1 (MQ=255)                                                                                                                                   cTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGa                 <  4:27510/149‑1 (MQ=255)                                                                                                                                    tCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAg                >  1:207079/1‑149 (MQ=255)                                                                                                                                        gACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCaa            <  4:652676/149‑1 (MQ=255)                                                                                                                                         aCCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAAc           >  3:620681/1‑149 (MQ=255)                                                                                                                                           cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAAc           >  4:150042/1‑147 (MQ=255)                                                                                                                                           cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAAc           <  3:150042/147‑1 (MQ=255)                                                                                                                                           cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGa         >  2:311798/1‑149 (MQ=255)                                                                                                                                               aCCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAtt     >  3:552395/1‑149 (MQ=255)                                                                                                                                                 cTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAg                                                                                                       <  2:227272/49‑1 (MQ=255)                                                                                                                                                 cTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAg                                                                                                       >  1:227272/1‑49 (MQ=255)                                                                                                                                                  ttGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTTGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAttttt  <  2:330166/149‑1 (MQ=255)                                                                                                                                                  ttGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAttttt  <  4:621112/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTT  >  NC_000913/1991935‑1992228

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