Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,155,440	A→G	6.4%	Q106R (CAG→CGG)	tmk →	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,155,440	0	A	G	6.4%	247.1 / 5.3	78	Q106R (CAG→CGG)	tmk	thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base A (42/31);  new base G (2/3);  total (44/34)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

TCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCG  >  NC_000913/1155293‑1155586                                                                                                                                                    |                                                                                                                                                   tCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCAGCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCa                                                                                                                                                    >  3:732454/1‑148 (MQ=255)  cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAgg                                                                                                                                                  >  2:284432/1‑149 (MQ=255)  cAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATATTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAg                                                                                                                                                   <  4:96483/148‑1 (MQ=255)    aaTCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGc                                                                                                                                                 >  4:221318/1‑148 (MQ=255)     aTCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGt                                                                                                                                               >  2:633825/1‑149 (MQ=255)         gTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATc                                                                                                                                            <  3:361038/148‑1 (MQ=255)           aGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAgg                                                                                                                                         <  2:258082/149‑1 (MQ=255)             gCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAggg                                                                                                                                        >  4:7188/1‑148 (MQ=255)                aTGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGggcggc                                                                                                                                    <  1:821698/149‑1 (MQ=255)                aTGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGggcgg                                                                                                                                     <  4:711767/148‑1 (MQ=255)                       cATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTg                                                                                                                             <  1:382299/149‑1 (MQ=255)                        aTTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTgg                                                                                                                            <  1:1069195/149‑1 (MQ=255)                         ttACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTGTCAGGGCGGCGGACGTgg                                                                                                                            >  3:446816/1‑148 (MQ=255)                         ttACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTgg                                                                                                                            >  2:200130/1‑148 (MQ=255)                         ttACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTgg                                                                                                                            <  1:200130/148‑1 (MQ=255)                          tACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGt                                                                                                                           >  4:350104/1‑148 (MQ=255)                            ccGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTAt                                                                                                                         <  4:311962/148‑1 (MQ=255)                              gATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGa                                                                                                                      >  2:980139/1‑149 (MQ=255)                               aTAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGAc                                                                                                                     >  2:372181/1‑149 (MQ=255)                                tAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGAc                                                                                                                     >  3:195493/1‑148 (MQ=255)                                  aaGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCa                                                                                                                   >  4:616904/1‑148 (MQ=255)                                      cGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCGGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACata                                                                                                              <  1:702799/149‑1 (MQ=255)                                       gAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAgg                                                                                                                                                  <  4:349059/112‑1 (MQ=255)                                       gAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAgg                                                                                                                                                  >  3:349059/1‑112 (MQ=255)                                           ttCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTg                                                                                                         <  1:372181/149‑1 (MQ=255)                                            tCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCGGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTg                                                                                                         >  3:281943/1‑148 (MQ=255)                                            tCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTg                                                                                                         <  3:7188/148‑1 (MQ=255)                                               gATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTgg                                                                                                                            <  4:45130/126‑1 (MQ=255)                                               gATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTgg                                                                                                                            >  3:45130/1‑126 (MQ=255)                                               gATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCa                                                                                                      >  4:403271/1‑148 (MQ=255)                                               gATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCa                                                                                                      >  4:553239/1‑148 (MQ=255)                                                aTGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCaa                                                                                                     >  4:1030507/1‑148 (MQ=255)                                                     ttATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTg                                                                                                <  4:517062/148‑1 (MQ=255)                                                     ttATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGc                                                                                               <  1:690414/149‑1 (MQ=255)                                                       aTGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCg                                                                                              <  4:779317/148‑1 (MQ=255)                                                              gcgcgTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTg                                                                                      >  2:799935/1‑149 (MQ=255)                                                                 cgTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGtt                                                                                    >  3:839578/1‑148 (MQ=255)                                                                   ttCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCGGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTtct                                                                                  >  3:477255/1‑148 (MQ=255)                                                                      aaCTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCggg                                                                              >  2:534624/1‑149 (MQ=255)                                                                       aCTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCggg                                                                              >  4:61470/1‑148 (MQ=255)                                                                         tGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGAt                                                                           >  2:358293/1‑149 (MQ=255)                                                                           gTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGAtt                                                                          <  3:506619/148‑1 (MQ=255)                                                                            tAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGAc                                                                                                                                <  2:119819/93‑1 (MQ=255)                                                                            tAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGAc                                                                                                                                >  1:119819/1‑93 (MQ=255)                                                                            tAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGAttt                                                                         >  3:766137/1‑148 (MQ=255)                                                                             aGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGAtttt                                                                        >  3:379120/1‑148 (MQ=255)                                                                              gAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCg                                                                      >  2:721110/1‑149 (MQ=255)                                                                                aaCGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTg                                                                                                         >  4:750209/1‑112 (MQ=255)                                                                                aaCGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTg                                                                                                         <  3:750209/112‑1 (MQ=255)                                                                                     tcatcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGa                                                                >  4:212855/1‑148 (MQ=255)                                                                                       atcaAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACtt                                                             >  1:521564/1‑149 (MQ=255)                                                                                          aaaCCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTaa                                                           <  3:1030507/148‑1 (MQ=255)                                                                                              cAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct                                                       >  4:253114/1‑148 (MQ=255)                                                                                              cAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGct                                                       <  3:567618/148‑1 (MQ=255)                                                                                                      gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCg                                               <  4:538441/148‑1 (MQ=255)                                                                                                      gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGa                                              <  1:284432/149‑1 (MQ=255)                                                                                                      gCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGa                                              <  2:45313/149‑1 (MQ=255)                                                                                                       cTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGa                                              <  3:253114/148‑1 (MQ=255)                                                                                                       cTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGAt                                             <  1:633825/149‑1 (MQ=255)                                                                                                             gCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTAc                                        >  4:673397/1‑148 (MQ=255)                                                                                                             gCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTAc                                        >  4:722656/1‑148 (MQ=255)                                                                                                              cACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTAccc                                      >  2:668459/1‑149 (MQ=255)                                                                                                                    ggTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAg                                >  1:23964/1‑149 (MQ=255)                                                                                                                      tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGt                               <  4:195493/148‑1 (MQ=255)                                                                                                                      tGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGt                               <  3:553239/148‑1 (MQ=255)                                                                                                                        aTTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTg                             >  3:515203/1‑148 (MQ=255)                                                                                                                        aTTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATAGGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTg                             <  3:403271/148‑1 (MQ=255)                                                                                                                         ttGGCGATCGCCACGATCTCTCCACTCGGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACAGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTgg                            <  4:984955/148‑1 (MQ=255)                                                                                                                         ttGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGc                           >  2:1006805/1‑149 (MQ=255)                                                                                                                         ttGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGc                           >  1:961907/1‑149 (MQ=255)                                                                                                                              gATCGCCACGATCTCTCCACTCGGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTaa                       <  4:281943/148‑1 (MQ=255)                                                                                                                               aTCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTaaa                      >  3:181465/1‑148 (MQ=255)                                                                                                                                     aCGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAAcgcgc                >  3:602211/1‑148 (MQ=255)                                                                                                                                     aCGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAAcgcgc                >  3:601159/1‑148 (MQ=255)                                                                                                                                               cACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGcgcgcg     >  1:1077370/1‑149 (MQ=255)                                                                                                                                                 cTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCgcgg    <  4:181465/148‑1 (MQ=255)                                                                                                                                                  tCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgtt   <  2:61042/148‑3 (MQ=255)                                                                                                                                                  tCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgtt   >  1:61042/1‑146 (MQ=255)                                                                                                                                                   cAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCgcggcg  >  4:126007/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCAAATCGGTAGGCGATGAAGTCATTACCGATAAAGCCGAAGTTCTGATGTTTTATGCCGCGCGCGTTCAACTGGTAGAAACGGTCATCAAACCAGCTCTGGCTAACGGCACCTGGGTGATTGGCGATCGCCACGATCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCG  >  NC_000913/1155293‑1155586

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: ‑