Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,347,745	C→T	5.4%	E390K (GAA→AAA)	rnb ←	ribonuclease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,347,745	0	C	T	5.4%	241.7 / 2.4	74	E390K (GAA→AAA)	rnb	ribonuclease II
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base T (2/2);  total (37/37)
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.46e-01

ACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAATTTGCG  >  NC_000913/1347600‑1347891                                                                                                                                                  |                                                                                                                                                   aCATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTcgac                                                                                                                                                 <  3:271495/149‑1 (MQ=255)     gCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTTGACGATa                                                                                                                                             >  6:398636/1‑149 (MQ=255)     gCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGAAGCCAACACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTTGACGATa                                                                                                                                             >  6:404717/1‑149 (MQ=255)      cAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGCTCAAAGCCCATATGCACGTTATAGATGCCAAACCCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATAc                                                                                                                                            >  6:365857/1‑149 (MQ=255)      cAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATa                                                                                                                                             <  7:150021/148‑1 (MQ=255)       aGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACg                                                                                                                                           >  8:376570/1‑149 (MQ=255)           cGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTg                                                                                                                                       >  7:399140/1‑149 (MQ=255)            gTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGCGCTTATCGCGAAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTgg                                                                                                                                      >  8:190779/1‑149 (MQ=255)             tGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTgg                                                                                                                                      >  5:251789/1‑148 (MQ=255)             tGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGc                                                                                                                                     <  7:376570/149‑1 (MQ=255)               cGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCaa                                                                                                                                   <  6:93063/149‑1 (MQ=255)                     cAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGcga                                                                                                                              <  2:174626/148‑1 (MQ=255)                          aCGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGc                                                                                                                        >  2:159790/1‑149 (MQ=255)                           cGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCt                                                                                                                       >  8:411105/1‑149 (MQ=255)                                  aGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGcgac                                                                                                                <  6:211035/149‑1 (MQ=255)                                  aGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGcgac                                                                                                                >  5:211035/1‑149 (MQ=255)                                       gTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGata                                                                                                            <  4:59785/148‑1 (MQ=255)                                       gTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGata                                                                                                            >  3:59785/1‑148 (MQ=255)                                        tCGGCATTCGCCGGATCAAAGCCCATATGCCGGTTATAGATGCCACAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATc                                                                                                          >  6:116306/1‑149 (MQ=255)                                            cATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGc                                                                                                      <  1:159790/149‑1 (MQ=255)                                             aTTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTTGCAATGCGACGAGGCTCGGCGACGATATCCAGCa                                                                                                     <  5:17610/149‑1 (MQ=255)                                              ttCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCAc                                                                                                    >  5:343808/1‑149 (MQ=255)                                                 gCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACtt                                                                                                  >  7:27406/1‑148 (MQ=255)                                                      aTCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCtt                                                                                            >  4:76900/1‑149 (MQ=255)                                                        cAAAGCCCATATGAACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCtttt                                                                                          >  3:98652/1‑149 (MQ=255)                                                         aaaGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGcgac                                                                                                                             <  1:7204/113‑1 (MQ=255)                                                         aaaGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGcgac                                                                                                                             >  2:7204/1‑113 (MQ=255)                                                         aaaGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGTGACGATATCCAGCACTTCACCttttt                                                                                         >  4:47058/1‑149 (MQ=255)                                                           aGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCa                                                                                       >  3:136297/1‑149 (MQ=255)                                                                atatGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATg                                                                                                                                 >  7:215701/1‑102 (MQ=255)                                                                atatGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATg                                                                                                                                 <  8:215701/102‑1 (MQ=255)                                                                  atGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCgaga                                                                                 >  8:385804/1‑148 (MQ=255)                                                                       cGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAg                                                                           <  7:318168/149‑1 (MQ=255)                                                                       cGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAg                                                                           <  7:411105/149‑1 (MQ=255)                                                                         ttATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCg                                                                         <  2:155071/149‑1 (MQ=255)                                                                         ttATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCg                                                                         <  2:114556/149‑1 (MQ=255)                                                                         ttATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCg                                                                         <  3:76900/149‑1 (MQ=255)                                                                            taGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGt                                                                       <  5:117987/148‑1 (MQ=255)                                                                                tGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATcc                                                                  <  7:385804/149‑1 (MQ=255)                                                                                  ccAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCgg                                                                >  4:223709/1‑149 (MQ=255)                                                                                   ccaaaCCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAAt                                                                                                                                  >  6:405944/3‑82 (MQ=255)                                                                                   ccaaaCCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAAt                                                                                                                                  <  5:405944/80‑1 (MQ=255)                                                                                   cAACACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGCCGAGGCTCGGCGACGATATCCAGCAGTTCACCTTTGTCACCGAGATTAAAGCGGTAATCCggg                                                               >  4:103686/1‑149 (MQ=255)                                                                                   cAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCggg                                                               >  7:109775/1‑149 (MQ=255)                                                                                    aaaaCCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGc                                                                                                                                >  5:61944/1‑83 (MQ=255)                                                                                    aaaaCCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGc                                                                                                                                <  6:61944/83‑1 (MQ=255)                                                                                    aaaaCCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGGTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCCCCTTTTTCACCGAGAATAAAGCGGTAATCCGGGc                                                              >  7:321013/1‑149 (MQ=255)                                                                                           aGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGTGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCttt                                                       <  3:47058/149‑1 (MQ=255)                                                                                           aGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCAc                                                                                               >  5:314386/1‑109 (MQ=255)                                                                                           aGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCAc                                                                                               <  6:314386/109‑1 (MQ=255)                                                                                                aTCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGTCGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCATCACTTCACCTTTTTCAGCGAGAATAAAGCGGTAATCCGGGCGAGCTTTaaa                                                    <  4:98652/147‑1 (MQ=255)                                                                                                aTCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGTCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTaaa                                                    <  8:200887/147‑1 (MQ=255)                                                                                                aTCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTaaa                                                    <  8:27406/147‑1 (MQ=255)                                                                                                 tCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCACACATCGCTTCGTCGACGATACGGTTGGCAATGCGACGAGGCTCGGAGACGATATCCAGCCCTTCACCTTATTCACCGAGAATAAAGCGGTAATCCGCGCGATCTTTAAacac                                                 >  4:58975/1‑149 (MQ=255)                                                                                                  cgcgCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAacac                                                 >  7:37953/1‑148 (MQ=255)                                                                                                     gcAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGt                                             <  4:136297/149‑1 (MQ=255)                                                                                                        gTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTaaa                                                    <  1:207387/139‑1 (MQ=255)                                                                                                        gTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTaaa                                                    >  2:207387/1‑139 (MQ=255)                                                                                                           cgcgGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTTGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTgg                                       <  5:398636/149‑1 (MQ=255)                                                                                                           cgcgGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTTGACGATACGGTTGGCAATGCGACGAGGCGCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTgg                                       <  5:404717/149‑1 (MQ=255)                                                                                                              gggCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTa                                    <  6:251789/149‑1 (MQ=255)                                                                                                                 ccGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTTGTTATGa                                 <  8:399140/149‑1 (MQ=255)                                                                                                                   gCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTaaa                                                    <  1:334414/128‑1 (MQ=255)                                                                                                                   gCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTaaa                                                    >  2:334414/1‑128 (MQ=255)                                                                                                                    cacaAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGc                              <  7:190779/149‑1 (MQ=255)                                                                                                                     acaAATGTTAGCGGCAATCGACGCTTCTTCGACGAAACGCTAGGCAATGCGACGAGGCACGGCGCCGATATCCCGCCCTGCGCCTTTTTCGCCGCGCGTAGAGCGGGGATCCGGGCGCGCTGTAAACAGCAGTGCGTGCGTATGACGcc                             >  8:171625/1‑149 (MQ=255)                                                                                                                     acaAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCga                                                                                   <  1:90067/95‑1 (MQ=255)                                                                                                                     acaAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCga                                                                                   >  2:90067/1‑95 (MQ=255)                                                                                                                      caAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCa                            <  6:343808/149‑1 (MQ=255)                                                                                                                       aaaTGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCAc                           >  6:359549/1‑149 (MQ=255)                                                                                                                                cGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCgcgg                  >  2:156234/1‑149 (MQ=255)                                                                                                                                     atcatcGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTg             <  5:154164/149‑1 (MQ=255)                                                                                                                                      tcatcGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTa                                    >  7:375786/1‑125 (MQ=255)                                                                                                                                      tcatcGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTa                                    <  8:375786/125‑1 (MQ=255)                                                                                                                                           gcgtcttcGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCaaa        <  5:116306/145‑1 (MQ=255)                                                                                                                                           gCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCaaa        <  3:103686/148‑1 (MQ=255)                                                                                                                                             ttcttcGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAAttt     >  2:418251/1‑149 (MQ=255)                                                                                                                                             ttcttcGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAAttt     >  2:71431/1‑149 (MQ=255)                                                                                                                                              tcttcGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAATTTg    >  3:160676/1‑149 (MQ=255)                                                                                                                                                ttcGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAATTTgcg  >  1:140554/1‑149 (MQ=255)                                                                                                                                                  cgacgaTACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAATTTgc   <  1:418251/146‑1 (MQ=255)                                                                                                                                                  cgacgaTACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAATTTgc   <  1:71431/146‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   ACATGCAGACCGTGCGTTTTCAGCAACGCTGCCAGCGCGTCGGCATTCGCCGGATCAAAGCCCATATGCACGTTATAGATGCCAAAACCGAGCTTATCGCGCAGTACGCGGGCCGCACAAATGTTAGCGGCAATCATCGCTTCTTCGACGATACGGTTGGCAATGCGACGAGGCTCGGCGACGATATCCAGCACTTCACCTTTTTCACCGAGAATAAAGCGGTAATCCGGGCGATCTTTAAACACCAGTGCGTGGTTATGACGCCACTCGCCGCGGCGTTGGCAAATTTGCG  >  NC_000913/1347600‑1347891

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