Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,419	T→A	5.7%	intergenic (‑91/+119)	tyrV ← / ← tyrT	tRNA‑Tyr/tRNA‑Tyr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,419	0	T	A	5.7%	245.4 / 2.9	70	intergenic (‑91/+119)	tyrV/tyrT	tRNA‑Tyr/tRNA‑Tyr
Reads supporting (aligned to +/- strand):  ref base T (35/31);  new base A (2/2);  total (37/33)
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.85e-01

GATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCG  >  NC_000913/1287278‑1287566                                                                                                                                              |                                                                                                                                                     gATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCAGTTGCCTGTGGCaa                                                                                                                                               >  2:665257/1‑149 (MQ=255)  aTGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAAc                                                                                                                                              <  4:512105/149‑1 (MQ=255)  aTGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTCGAGCTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAAc                                                                                                                                              >  2:472881/1‑149 (MQ=255)        gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctc                                                                                                                                        >  4:160741/1‑149 (MQ=255)        gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACgg                                                                                                                                            <  2:314589/145‑2 (MQ=255)        gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACgg                                                                                                                                            >  1:314589/1‑144 (MQ=255)            aTTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGc                                                                                                                                    <  2:451073/149‑1 (MQ=255)             tttACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCt                                                                                                                                   >  4:315864/1‑149 (MQ=255)                    tCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtc                                                                                                                            >  1:339569/1‑149 (MQ=255)                    tCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtc                                                                                                                            >  2:416491/1‑149 (MQ=255)                           ccTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcga                                                                                                                     >  4:361506/1‑149 (MQ=255)                            cTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcgaa                                                                                                                    <  1:537289/149‑1 (MQ=255)                             tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaa                                                                                                                         >  3:342576/1‑143 (MQ=255)                             tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAc                                                                                                                   <  3:315864/149‑1 (MQ=255)                             tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAc                                                                                                                   >  3:680828/1‑149 (MQ=255)                             tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGCTCTCTCGCTTTCGCtcgaa                                                                                                                         <  4:342576/143‑1 (MQ=255)                                   cGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGt                                                                                                             <  3:361506/149‑1 (MQ=255)                                        gggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAg                                                                                                        <  3:160741/149‑1 (MQ=255)                                        gggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAg                                                                                                        <  4:66979/149‑1 (MQ=255)                                         ggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAgg                                                                                                       >  2:439298/1‑149 (MQ=255)                                              cccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTtctc                                                                                                  >  1:14996/1‑149 (MQ=255)                                                           aTTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGAt                                                                                     <  2:14996/149‑1 (MQ=255)                                                            ttCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATg                                                                                    >  3:681058/1‑149 (MQ=255)                                                            ttCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGAGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATg                                                                                    >  3:571094/1‑149 (MQ=255)                                                             tCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGa                                                                                   <  1:439298/149‑1 (MQ=255)                                                                  ttttGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCa                                                                              >  2:84143/1‑149 (MQ=255)                                                                         gTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTc                                                                       >  1:463767/1‑149 (MQ=255)                                                                         gTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTc                                                                       <  4:680828/149‑1 (MQ=255)                                                                           aaTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcac                                                                     >  1:481073/1‑149 (MQ=255)                                                                           aaTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcac                                                                     >  2:466239/1‑149 (MQ=255)                                                                                ttGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAAtct                                                                >  1:450263/1‑149 (MQ=255)                                                                                     aTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACttt                                                                        >  3:633539/1‑136 (MQ=255)                                                                                     aTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACttt                                                                        <  4:633539/136‑1 (MQ=255)                                                                                       ggtggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑aa                                                        >  4:101755/1‑149 (MQ=255)                                                                                          ggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGtt                                                     <  4:681058/149‑1 (MQ=255)                                                                                          ggtggGAGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGtt                                                     <  4:571094/149‑1 (MQ=255)                                                                                              ggggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACca                                                 <  1:665257/149‑1 (MQ=255)                                                                                                ggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACcaca                                               >  4:503775/1‑149 (MQ=255)                                                                                                     gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                          >  3:232422/1‑149 (MQ=255)                                                                                                     gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                          <  3:503775/149‑1 (MQ=255)                                                                                                        tattCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCtgtg                                       <  1:416491/149‑1 (MQ=255)                                                                                                          ttCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTgtggt                                     >  1:291376/1‑149 (MQ=255)                                                                                                             gtcgCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGaa                                  >  2:629146/1‑149 (MQ=255)                                                                                                             gtcgCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGaa                                  <  2:339569/149‑1 (MQ=255)                                                                                                                gcttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTca                                                                      <  1:131359/111‑1 (MQ=255)                                                                                                                gcttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTca                                                                      >  2:131359/1‑111 (MQ=255)                                                                                                                 cttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAt                              <  2:291376/149‑1 (MQ=255)                                                                                                                  ttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtg                             <  2:463767/149‑1 (MQ=255)                                                                                                                    cgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                           <  1:84143/149‑1 (MQ=255)                                                                                                                        cctcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggtggtg                       <  2:481073/149‑1 (MQ=255)                                                                                                                         ctcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggtggtgg                      <  2:450263/149‑1 (MQ=255)                                                                                                                          tcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGATCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTggg                     >  1:181854/1‑149 (MQ=255)                                                                                                                           cACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTgggg                    >  1:429576/1‑149 (MQ=255)                                                                                                                            aCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                          <  2:289795/126‑1 (MQ=255)                                                                                                                            aCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                          >  1:289795/1‑126 (MQ=255)                                                                                                                              ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                           >  2:305009/1‑139 (MQ=255)                                                                                                                              ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                           <  1:305009/139‑1 (MQ=255)                                                                                                                              ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaa                 >  3:326722/1‑149 (MQ=255)                                                                                                                              ccTTCGGGACGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaa                 >  3:636653/1‑149 (MQ=255)                                                                                                                               cTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAg                >  3:671747/1‑149 (MQ=255)                                                                                                                                  gagaacGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaa                 <  4:600505/140‑1 (MQ=255)                                                                                                                                  gagaacGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGCGAATTATGGTGGTGGGGGaa                 >  3:600505/6‑145 (MQ=255)                                                                                                                                  cGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGAt             <  2:429576/149‑1 (MQ=255)                                                                                                                                    ggCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCATAATTCACCACAGCGATgtggt                                     <  1:312655/124‑1 (MQ=255)                                                                                                                                      acGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGa         >  4:242893/2‑149 (MQ=255)                                                                                                                                       cGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGaa        <  4:671747/149‑1 (MQ=255)                                                                                                                                         ttGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAAcc      <  3:101755/149‑1 (MQ=255)                                                                                                                                            ccTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTc   <  1:466239/149‑1 (MQ=255)                                                                                                                                             cTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCg  >  2:144272/1‑149 (MQ=255)                                                                                                                                             cTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGATCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAGGGTGGTGGGGGAAGGATTCGAACCTTCg  >  2:642880/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     GATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCG  >  NC_000913/1287278‑1287566

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