Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	217,587	G→A	6.0%	P397S (CCG→TCG)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	217,587	0	G	A	6.0%	205.4 / 2.1	66	P397S (CCG→TCG)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (30/32);  new base A (2/2);  total (32/34)
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.76e-01

CGTAGCAACCCATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAACGGTACGGTCAAT  >  NC_000913/217443‑217731                                                                                                                                                 |                                                                                                                                                 cGTAGCAACCCATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCtt                                                                                                                                              <  2:209088/149‑1 (MQ=255) cGTAGCAACCCATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCtt                                                                                                                                              <  2:140186/149‑1 (MQ=255)     gCAACCCATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGAt                                                                                                                                          <  6:335846/149‑1 (MQ=255)         cccATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTg                                                                                                                                             >  8:359313/1‑142 (MQ=255)         cccATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTg                                                                                                                                             <  7:359313/142‑1 (MQ=255)            aTCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCgccagc                                                                                                                                   <  2:12557/149‑1 (MQ=255)                  aGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCcaccac                                                                                                                             >  3:269177/1‑149 (MQ=255)                    gATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGt                                                                                                                           >  7:335109/1‑149 (MQ=255)                         ggTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACg                                                                                                                       <  7:63461/148‑1 (MQ=255)                           ttACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGAt                                                                                                                                          >  7:331510/1‑127 (MQ=255)                           ttACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGAt                                                                                                                                          <  8:331510/127‑1 (MQ=255)                           ttACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACCTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGa                                                                                                                     >  8:118982/1‑148 (MQ=255)                              cGGCCATCTTCACCCGGTACGGCGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGatat                                                                                                                  <  2:299088/148‑1 (MQ=255)                               ggCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCt                                                                                                                >  1:213436/1‑149 (MQ=255)                                 ccATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTg                                                                                                               >  2:134761/1‑148 (MQ=255)                                 ccATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTg                                                                                                               >  4:317331/1‑148 (MQ=255)                                    tCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGAGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAAc                                                                                                           >  2:347821/1‑149 (MQ=255)                                          cccTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCgg                                                                                                     <  8:141734/149‑1 (MQ=255)                                           ccTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCCCGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCggg                                                                                                    >  4:137934/1‑149 (MQ=255)                                                  ggAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCAATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGa                                                                                             >  6:59012/1‑149 (MQ=255)                                                   gAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGAAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGa                                                                                             >  1:232870/1‑148 (MQ=255)                                                            aGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGAt                                                                                    <  6:736/148‑1 (MQ=255)                                                            aGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGAt                                                                                    <  8:335109/148‑1 (MQ=255)                                                            aGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATc                                                                                   >  2:285333/1‑149 (MQ=255)                                                               gCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCa                                                                                <  1:317793/149‑1 (MQ=255)                                                                   cGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGt                                                                              <  4:365181/147‑1 (MQ=255)                                                                   cGGAGTACTTGGTACACAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGtt                                                                             >  4:328402/1‑148 (MQ=255)                                                                    ggAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGa                                                                           <  2:213436/149‑1 (MQ=255)                                                                    ggAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGAGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGCTAGCGACATCGCGATCCCAGTTGa                                                                           <  1:347821/149‑1 (MQ=255)                                                                               tACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagtag                                                                <  5:59012/149‑1 (MQ=255)                                                                                 cccAGCTGGAAGTTGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAgtg                                                              <  7:118982/149‑1 (MQ=255)                                                                                  ccAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAgtgt                                                             <  1:134761/149‑1 (MQ=255)                                                                                  ccAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAgtgt                                                             >  1:350926/1‑149 (MQ=255)                                                                                   cAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                                                            <  4:269177/149‑1 (MQ=255)                                                                                   cAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGtt                                                            >  2:17655/1‑149 (MQ=255)                                                                                        ggAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagta                                                                 <  4:211890/139‑1 (MQ=255)                                                                                        ggAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAagta                                                                 >  3:211890/1‑139 (MQ=255)                                                                                           aGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTTGATACCCAGCCACCACGTAACCGATATCTCCACCTTCCGTGGTAGCGACAACGCGAACCAAGTTGATGCCGAAGTAGTCTTTACCATc                                                     >  4:58181/1‑148 (MQ=255)                                                                                             aTGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt                                                  <  3:328402/149‑1 (MQ=255)                                                                                                tgACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATACGGGCTTGGATCGCCAGCCACCACTTGACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCATTTGAGGCCGAAGAAGTGTTTACCATCGATGt                                                >  6:251725/1‑148 (MQ=255)                                                                                                tgACCAACTTCGATACCACGTTTGATCAGAAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGt                                                <  2:232870/148‑1 (MQ=255)                                                                                                      acctCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCAc                                          >  7:131867/4‑148 (MQ=255)                                                                                                      aCTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCAcc                                         <  2:110079/149‑1 (MQ=255)                                                                                                          cGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACcag                                       <  1:17655/147‑1 (MQ=255)                                                                                                            aTACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCagcagc                                   <  3:317331/149‑1 (MQ=255)                                                                                                                cACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGaaa                               >  6:106100/1‑149 (MQ=255)                                                                                                                         tCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt                                                  <  5:296110/121‑1 (MQ=255)                                                                                                                         tCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGAt                                                  >  6:296110/1‑121 (MQ=255)                                                                                                                             cagcCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATcgccg                  >  1:210214/1‑149 (MQ=255)                                                                                                                               gcCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCCGCGAAATCACTCATcgccgc                 >  3:21130/1‑148 (MQ=255)                                                                                                                                ccTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCaa               <  5:106100/149‑1 (MQ=255)                                                                                                                                    ggcTGGCGATCCAAGCCCGGATGGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCa                            <  8:54883/130‑1 (MQ=255)                                                                                                                                    ggcTGGCGATCCAAGCCCGGATGGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCa                            >  7:54883/3‑132 (MQ=255)                                                                                                                                     ccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggta          <  2:210214/149‑1 (MQ=255)                                                                                                                                     ccTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggta          <  2:350926/149‑1 (MQ=255)                                                                                                                                       tGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggtacg        >  7:167863/1‑149 (MQ=255)                                                                                                                                       tGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggtacg        >  4:266957/1‑149 (MQ=255)                                                                                                                                       tGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggtacg        >  8:176944/1‑149 (MQ=255)                                                                                                                                       tGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggtacg        >  3:137514/1‑149 (MQ=255)                                                                                                                                         gCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggtacg        <  3:266957/147‑1 (MQ=255)                                                                                                                                         gCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAacggtacg        <  7:176944/147‑1 (MQ=255)                                                                                                                                           cATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAACGGTACGGTc     >  3:165160/1‑148 (MQ=255)                                                                                                                                             tCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGCCATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAACGGTACGGTCAAt  >  5:254833/1‑149 (MQ=255)                                                                                                                                             tCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATc                                                     >  5:3174/1‑98 (MQ=255)                                                                                                                                             tCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATc                                                     <  6:3174/98‑1 (MQ=255)                                                                                                                                             tCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAACGGTACGGTCAAt  <  8:131867/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CGTAGCAACCCATCGTCAGGATTTGGTTACGGCCATCTTCACCCTGTACGGAGGCTTTCAGTGCTTCGGAGTACTTGGTACCCAGCTGGAAGATGTGACCAACTTCGATACCACGTTTGATCAGCAGCCTACCCTGGCCATCCGGGCTTGGATCGCCAGCCACCACGTTACGGATATCTGCAACTTCCGGGGTAGCGACATCGCGATCCCAGTTGATGCCGAAGTAGTGTTTACCATCGATGTTAGCACCAGCAGCGAAATCACTCATCGCCGCAACGGTACGGTCAAT  >  NC_000913/217443‑217731

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