Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,552,185	G→C	6.7%	A184P (GCG→CCG)	uxuA →	mannonate hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,552,185	0	G	C	6.7%	264.9 / 8.8	76	A184P (GCG→CCG)	uxuA	mannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base G (39/31);  new base C (3/2);  total (43/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.11e-01

TTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGT  >  NC_000913/4552041‑4552330                                                                                                                                                 |                                                                                                                                                  ttCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGga                                                                                                                                               >  4:286350/1‑149 (MQ=255) ttCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCATAACATCATTGCTGGTCTTCCGGGCGCGga                                                                                                                                               >  2:372226/1‑149 (MQ=255)     aaaTGCATATCCTGAAACGCCCAGGCGCGGGAGCGGATTACACCGAAGAACAAATTGCTCAGGCAGCTGAACGCTTCGCCACCATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaagaa                                                                                                                                           <  3:384888/149‑1 (MQ=255)     aaaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGgaagaa                                                                                                                                           <  3:343073/149‑1 (MQ=255)     aaaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagac                                                                                                                                           <  4:141310/149‑2 (MQ=255)      aaTGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGCGCCCGgaagacc                                                                                                                                          <  3:61816/149‑3 (MQ=255)        tGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCTATGAAGATATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAATACAAAGCGCGTCTGACCCGTAGCATCATTGCTGGTCTTCCGGGCTCGGAAGAAGGc                                                                                                                                        >  4:382822/1‑149 (MQ=255)           atatCCTGAAACGCCCAGGCGCGGAAGCGGGTTGCACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGACTTCCGGGCGCGGAAGAAGGCTac                                                                                                                                     >  2:543645/1‑149 (MQ=255)              tCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACAcc                                                                                                                                  >  3:29409/1‑149 (MQ=255)                cTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCt                                                                                                                                >  3:1802/1‑149 (MQ=255)                 tGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTc                                                                                                                               >  3:563628/1‑149 (MQ=255)                           aGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTcc                                                                                                                     <  3:13431/149‑1 (MQ=255)                                     gCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCt                                                                                                           >  4:366707/1‑149 (MQ=255)                                     gCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCt                                                                                                           >  2:24741/1‑149 (MQ=255)                                     gCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCt                                                                                                           <  4:29409/149‑1 (MQ=255)                                      cGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTg                                                                                                          <  4:347915/149‑1 (MQ=255)                                        gATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGa                                                                                                        <  4:563628/149‑1 (MQ=255)                                          ttACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGc                                                                                                      >  4:433900/1‑149 (MQ=255)                                          ttACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGc                                                                                                      >  2:348961/1‑149 (MQ=255)                                           tACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCt                                                                                                     <  1:139096/149‑1 (MQ=255)                                             cacCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGt                                                                                                   <  1:619095/149‑1 (MQ=255)                                              acCGAAGAAGAAATTGATCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAACGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGTAGCTGTa                                                                                                  >  2:146390/1‑149 (MQ=255)                                               ccGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTAc                                                                                                 <  1:306743/149‑1 (MQ=255)                                                cGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACa                                                                                                >  4:101568/1‑149 (MQ=255)                                                cGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACa                                                                                                >  2:501285/1‑149 (MQ=255)                                                 gaagaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACaa                                                                                               >  4:301512/1‑149 (MQ=255)                                                 gaagaagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAAACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACaa                                                                                               >  1:190579/1‑149 (MQ=255)                                                     aagaaATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGat                                                                                           <  1:415379/149‑1 (MQ=255)                                                        aaaTTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATc                                                                                        <  3:58533/149‑1 (MQ=255)                                                                   gCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGaa                                                                             >  3:259130/1‑149 (MQ=255)                                                                        tgtaCGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTgc                                                                        <  3:547550/146‑1 (MQ=255)                                                                        tGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACAccc                                                                                                                                 <  1:95866/92‑1 (MQ=255)                                                                        tGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACAccc                                                                                                                                 >  2:95866/1‑92 (MQ=255)                                                                               ttCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGaaaa                                                                 >  2:596984/1‑149 (MQ=255)                                                                                         tGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTc                                                       >  2:48561/1‑149 (MQ=255)                                                                                          gAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCt                                                                                                     <  1:116050/102‑1 (MQ=255)                                                                                          gAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCt                                                                                                     >  2:116050/1‑102 (MQ=255)                                                                                           aGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCtt                                                     <  3:505489/149‑1 (MQ=255)                                                                                                  aaGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGaaa                                              <  1:372226/149‑1 (MQ=255)                                                                                                     aCAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCg                                           >  1:551498/1‑149 (MQ=255)                                                                                                      cAAAGCGCGTCTGATCCGTAAAATCATTGATTGTCTTCCGGTCGCGGAAGAAGGCTACACCCTCGACCAGATGCGTAAACACCTGGAGCTGTACAAAGATATCGCCAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                          <  3:382822/149‑1 (MQ=255)                                                                                                      cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                          <  1:370410/149‑1 (MQ=255)                                                                                                      cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                          <  4:1802/149‑1 (MQ=255)                                                                                                      cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                          <  3:169748/149‑1 (MQ=255)                                                                                                      cAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGa                                          <  1:467421/149‑1 (MQ=255)                                                                                                       aaaGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGat                                         >  3:258804/1‑149 (MQ=255)                                                                                                              gTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                                  <  3:101568/149‑1 (MQ=255)                                                                                                               tCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAg                                 >  2:126437/1‑149 (MQ=255)                                                                                                                     ccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTg                           <  1:24741/149‑1 (MQ=255)                                                                                                                     ccGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTg                           >  1:484820/1‑149 (MQ=255)                                                                                                                          aCATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTgaagaa                      <  3:286350/149‑1 (MQ=255)                                                                                                                            atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                                  >  3:320099/1‑135 (MQ=255)                                                                                                                            atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCa                                  <  4:320099/135‑1 (MQ=255)                                                                                                                            atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt                    >  3:593865/1‑149 (MQ=255)                                                                                                                            atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt                    >  4:568133/1‑149 (MQ=255)                                                                                                                            atcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGt                    >  3:447801/1‑149 (MQ=255)                                                                                                                             tcatTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTc                   >  3:447089/1‑149 (MQ=255)                                                                                                                              catTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCg                  <  2:551498/149‑1 (MQ=255)                                                                                                                               atTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCgg                 >  2:82637/1‑149 (MQ=255)                                                                                                                                ttGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgc                <  3:301512/149‑1 (MQ=255)                                                                                                                                  gCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt              <  1:282814/149‑1 (MQ=255)                                                                                                                                  gCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgt              >  2:9694/1‑149 (MQ=255)                                                                                                                                    tGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAg                                             <  1:155430/116‑1 (MQ=255)                                                                                                                                    tGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAg                                             <  1:437086/116‑1 (MQ=255)                                                                                                                                    tGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAg                                             >  2:155430/1‑116 (MQ=255)                                                                                                                                    tGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAg                                             >  2:437086/1‑116 (MQ=255)                                                                                                                                    ccttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgc            >  1:345981/4‑149 (MQ=255)                                                                                                                                      ttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcgt          >  2:583521/2‑149 (MQ=255)                                                                                                                                      ttCTTCCGCGCCCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcgt          >  4:211392/2‑149 (MQ=255)                                                                                                                                      gTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                             <  4:497595/98‑1 (MQ=255)                                                                                                                                      gTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACttt                                                             >  3:497595/1‑98 (MQ=255)                                                                                                                                      gTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcgt          >  2:531788/1‑149 (MQ=255)                                                                                                                                      gTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGgcgtgcgt          >  3:301677/1‑149 (MQ=255)                                                                                                                                       tCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGGAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTGGCTGAAGAAGTCGGCGTGCGTa         >  3:379672/1‑149 (MQ=255)                                                                                                                                         ttCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATg       >  3:152440/1‑149 (MQ=255)                                                                                                                                              ggCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGt  >  1:186256/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TTCGAAATGCATATCCTGAAACGCCCAGGCGCGGAAGCGGATTACACCGAAGAAGAAATTGCTCAGGCAGCTGAACGCTTCGCCACTATGAGCGATGAAGACAAAGCGCGTCTGACCCGTAACATCATTGCTGGTCTTCCGGGCGCGGAAGAAGGCTACACCCTCGACCAGTTCCGTAAACACCTGGAGCTGTACAAAGATATCGACAAAGCGAAGCTGCGCGAAAACTTTGCCGTCTTCCTGAAAGCGATTATTCCAGTTGCTGAAGAAGTCGGCGTGCGTATGGCTGT  >  NC_000913/4552041‑4552330

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