Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,420	G→C	5.7%	intergenic (‑92/+118)	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,420	0	G	C	5.7%	248.5 / 3.3	70	intergenic (‑92/+118)	tyrV/tyrT	tRNA‑Tyr/tRNA‑Tyr
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (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 = 1.00e+00

GATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAA  >  NC_000913/1287278‑1287568                                                                                                                                               |                                                                                                                                                      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                             <  1:305009/139‑1 (MQ=255)                                                                                                                              ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                             >  2:305009/1‑139 (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)                                                                                                                                               gtgGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGaa  <  4:636653/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                      GATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAA  >  NC_000913/1287278‑1287568

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