Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,992,083	G→C	8.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	8.4%	316.3 / 13.7	94	noncoding (35/76 nt)	glyW	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (40/46);  new base C (3/5);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTT  >  NC_000913/1991935‑1992228                                                                                                                                                     |                                                                                                                                                  aaCCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTgg                                                                                                                                                   >  3:720194/1‑149 (MQ=255)      gACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAgg                                                                                                                                              >  1:622360/1‑149 (MQ=255)        cTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                                                                                                                            <  3:227644/149‑1 (MQ=255)        cTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTc                                                                                                                                            <  4:360201/149‑1 (MQ=255)           gACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                                                                                                                         <  1:754049/149‑3 (MQ=255)           gACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                                                                                                                         <  2:228164/149‑3 (MQ=255)           gACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCggg                                                                                                                                         <  4:200664/149‑3 (MQ=255)              ggATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                                                                                                                      >  4:289038/1‑149 (MQ=255)              ggATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctc                                                                                                                                      >  4:117309/1‑149 (MQ=255)               gATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAATTCGCGACCCCGACCTTGGCAAGGTCGTGatct                                                                                                                                     <  4:485203/149‑1 (MQ=255)               gATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctct                                                                                                                                     <  3:695852/149‑1 (MQ=255)               gATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctct                                                                                                                                     <  1:587057/149‑1 (MQ=255)                aTTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCNACCCCGACCTTGGCNAGNNNGTGCTCTa                                                                                                                                    <  2:241442/149‑1 (MQ=255)                aTTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTa                                                                                                                                    <  1:418950/149‑1 (MQ=255)                 tttGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCgggg                                                                                                                                        <  3:26601/144‑4 (MQ=255)                 tttGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGCCAAGGTCgggg                                                                                                                                        >  4:26601/1‑141 (MQ=255)                    gCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCaa                                                                                                                                >  4:65154/1‑149 (MQ=255)                        tCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            <  4:210491/149‑1 (MQ=255)                        tCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGa                                                                                                                            <  1:504054/149‑1 (MQ=255)                          cGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGc                                                                                                                          <  1:669615/149‑1 (MQ=255)                           gCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCt                                                                                                                         >  4:551397/1‑149 (MQ=255)                                aTAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTcc                                                                                                                    <  2:622360/149‑1 (MQ=255)                                   aCCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGc                                                                                                                 <  4:325364/149‑1 (MQ=255)                                     cGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCAt                                                                                                               >  2:718456/1‑149 (MQ=255)                                        tttGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATtca                                                                                                            >  1:206283/1‑149 (MQ=255)                                            ttAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCaa                                                                                                        <  1:537327/149‑1 (MQ=255)                                            ttAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCaa                                                                                                        >  1:472855/1‑149 (MQ=255)                                             tAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAg                                                                                                       >  1:62966/1‑149 (MQ=255)                                              aaCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctct                                                                                                                                     <  1:442618/118‑1 (MQ=255)                                              aaCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGctct                                                                                                                                     >  2:442618/1‑118 (MQ=255)                                              aaCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGc                                                                                                      >  3:90900/1‑149 (MQ=255)                                               aCGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCa                                                                                                     >  1:686912/1‑149 (MQ=255)                                                 gcgcCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAAt                                                                                                   >  4:341234/1‑149 (MQ=255)                                                 gcgcCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATCCCCGCATTCATCAAGCAAt                                                                                                   <  3:445496/149‑1 (MQ=255)                                                 gcgcCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGGCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAAt                                                                                                   <  4:260830/149‑1 (MQ=255)                                                     cAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGt                                                                                               >  3:119779/1‑149 (MQ=255)                                                        aTTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCAATTCCCGCATTCATCAAGCAATCAGTTaa                                                                                            <  1:615185/149‑1 (MQ=255)                                                         ttcttcAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAAt                                                                                           <  2:62966/149‑1 (MQ=255)                                                            tacAGGCGTTTCAGCCAGATATCCGCATGACGCCGATGTCTTTTAAACTGGCACGGGAAACGAGACTCGAACTCTCGACCCCGACCTTTGAAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCAc                                                                                        >  2:85570/3‑149 (MQ=25)                                                                    tttCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCaa                                                                                                        >  4:753207/1‑125 (MQ=255)                                                                    tttCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCaa                                                                                                        <  3:753207/125‑1 (MQ=255)                                                                    tttCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAtttt                                                                                >  1:775421/1‑149 (MQ=255)                                                                     ttCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTtta                                                                               <  2:281256/149‑1 (MQ=255)                                                                     gtttaaaaGACATCGGCGTCAAGCGGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTtta                                                                               >  1:372477/8‑149 (MQ=21)                                                                         gCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTtattat                                                                           >  2:325650/1‑149 (MQ=255)                                                                          ccAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATc                                                                          <  4:414000/149‑1 (MQ=255)                                                                          ccAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATc                                                                          >  3:657758/1‑149 (MQ=255)                                                                          ccAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATc                                                                          <  1:750198/149‑1 (MQ=255)                                                                            aGACATCGGCGTCAAGCGGATGTCTTTTCAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGt                                                                        >  3:651507/1‑149 (MQ=21)                                                                            aGACATCGGCGTCAAGCGGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTaa                                                                                            <  3:81603/129‑1 (MQ=18)                                                                              aCATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCt                                                                      <  4:119779/149‑1 (MQ=255)                                                                               cATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTg                                                                     <  3:117309/149‑1 (MQ=255)                                                                               cATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTg                                                                     <  3:23509/149‑1 (MQ=255)                                                                               cATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTg                                                                     <  3:242147/149‑1 (MQ=255)                                                                               cATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCACGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTg                                                                     <  3:26017/149‑1 (MQ=255)                                                                                  ccGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCa                                                                  >  4:607044/1‑149 (MQ=255)                                                                                   cGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCaa                                                                 <  4:90900/149‑1 (MQ=255)                                                                                        gACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCt                                                                      >  4:733708/1‑139 (MQ=255)                                                                                        gACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCt                                                                      <  3:733708/139‑1 (MQ=255)                                                                                          cGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgc                                                          <  3:391161/149‑1 (MQ=255)                                                                                           gCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgcc                                                         >  2:132070/1‑149 (MQ=255)                                                                                           gCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTgcc                                                         >  2:703346/1‑149 (MQ=255)                                                                                               aTGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCATCTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATAACTTGATTTTATTATCGTCTGGCAATCAGTGccgccg                                                     >  1:759173/1‑149 (MQ=255)                                                                                                    ttttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGa                                                >  1:433957/1‑149 (MQ=255)                                                                                                    ttttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGa                                                >  2:386699/1‑149 (MQ=255)                                                                                                      ttAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATg                                              <  2:317817/149‑1 (MQ=255)                                                                                                          aCTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGtt                                          <  1:148987/149‑1 (MQ=255)                                                                                                          aCTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGtt                                          <  2:177580/149‑1 (MQ=255)                                                                                                            tGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGa                                                                                    <  1:639560/105‑1 (MQ=255)                                                                                                            tGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGa                                                                                    >  2:639560/1‑105 (MQ=255)                                                                                                            tGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCACTACCAACTGAGCTATTCCCGCACTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGc                                        <  2:206283/149‑1 (MQ=255)                                                                                                                gCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTc                                    <  3:289038/149‑1 (MQ=255)                                                                                                                 cGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCt                                   <  3:3914/149‑1 (MQ=255)                                                                                                                  gggAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTa                                  <  2:686912/149‑1 (MQ=255)                                                                                                                      aaCGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTa                              >  4:626393/1‑149 (MQ=255)                                                                                                                         gagaCTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCt                           >  2:706825/1‑149 (MQ=255)                                                                                                                              tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAtt                                                                                  >  3:502978/1‑89 (MQ=255)                                                                                                                              tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGAtt                                                                                  <  4:502978/89‑1 (MQ=255)                                                                                                                              tCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCtagcat                      >  4:460487/1‑144 (MQ=255)                                                                                                                                     cgcgACCCCGACCTTGGCAAGGTCGTGCTCTATCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGACGCGATGAGt               >  4:649771/1‑149 (MQ=255)                                                                                                                                     cgcgACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGt               <  3:391967/149‑1 (MQ=255)                                                                                                                                      gcaaCCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTc              >  1:252468/4‑149 (MQ=255)                                                                                                                                        gACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTc                                    <  3:152119/125‑1 (MQ=255)                                                                                                                                        gACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTc                                    >  4:152119/1‑125 (MQ=255)                                                                                                                                          gcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGAt                   >  3:186216/4‑140 (MQ=255)                                                                                                                                          gcacGACCTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGAt                   <  4:186216/137‑1 (MQ=255)                                                                                                                                           cccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGa         >  2:327963/1‑149 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGat        >  1:712343/1‑149 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGat        >  1:670082/1‑149 (MQ=255)                                                                                                                                            ccGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGat        >  2:559597/1‑149 (MQ=255)                                                                                                                                              gACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGatat      <  3:341550/149‑1 (MQ=255)                                                                                                                                                ccTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAttt    <  4:720194/149‑1 (MQ=255)                                                                                                                                                ccTTGCCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAGGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAttt    >  4:303864/1‑149 (MQ=255)                                                                                                                                                  ttGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAttttt  <  1:88257/149‑1 (MQ=255)                                                                                                                                                  ttGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATAttttt  <  1:706825/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  AACCGGACTAGACGGATTTGCAATCCGCTACATAACCGCTTTGTTAACGCGCCAAATTCTTCAGGCCTTTCAGCCAGACATCCGCTTGACGCCGATGTCTTTTAAACTGGAGCGGGAAACGAGACTCGAACTCGCGACCCCGACCTTGGCAAGGTCGTGCTCTACCAACTGAGCTATTCCCGCATTCATCAAGCAATCAGTTAATCACTTGATTTTATTATCGTCTGGCAATCAGTGCCGCCGTTCGATGCGTTGCATTCTACTTACCTGGCGCGATGAGTCAACGATATTTTT  >  NC_000913/1991935‑1992228

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