Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,760,788	C→G	6.2%	A154A (GCG→GCC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,760,788	0	C	G	6.2%	304.7 / 6.6	91	A154A (GCG→GCC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base C (42/42);  new base G (4/2);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCC  >  NC_000913/3760643‑3760935                                                                                                                                                  |                                                                                                                                                    tggtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAt                                                                                                                                                 >  2:925120/1‑150 (MQ=255)   gtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTg                                                                                                                                               <  1:302381/150‑1 (MQ=255)   gtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTg                                                                                                                                               >  1:936826/1‑150 (MQ=255)   gtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTg                                                                                                                                               >  2:343714/1‑150 (MQ=255)    tCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGt                                                                                                                                              >  2:291178/1‑150 (MQ=255)     cccTACTTCGTGTGGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGtt                                                                                                                                             >  2:1002972/1‑150 (MQ=255)     cccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGtt                                                                                                                                             >  2:900425/1‑150 (MQ=255)      ccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCAGAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTa                                                                                                                                            >  1:1362316/1‑150 (MQ=255)       cTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAt                                                                                                                                           >  1:302351/1‑150 (MQ=255)       cTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAt                                                                                                                                           >  1:575864/1‑150 (MQ=255)         aCTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTg                                                                                                                                         >  1:878994/1‑150 (MQ=255)           ttCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCTTTGTTATTGAc                                                                                                                                       <  2:1322217/150‑1 (MQ=255)            tCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACg                                                                                                                                      >  2:1109371/1‑150 (MQ=255)               tgtAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATa                                                                                                                                   >  2:335165/1‑150 (MQ=255)                 tAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACa                                                                                                                                 >  2:1391044/1‑150 (MQ=255)                   gggTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAgg                                                                                                                               >  2:947578/1‑150 (MQ=255)                           cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTcc                                                                                                                       <  1:480324/150‑1 (MQ=255)                            ctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCg                                                                                                                      <  1:407154/150‑1 (MQ=255)                             tgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGc                                                                                                                     <  2:1106320/150‑1 (MQ=255)                              gcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGcc                                                                                                                    <  1:1159531/150‑1 (MQ=255)                              gcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGcc                                                                                                                    <  2:1068461/150‑1 (MQ=255)                               cctgACGCATAACATCGGGAATACGAAACGCACCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGccc                                                                                                                   <  1:397220/150‑1 (MQ=255)                                 tgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGccccc                                                                                                                 <  2:453176/150‑1 (MQ=255)                                 tgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGccccc                                                                                                                 <  2:690767/150‑1 (MQ=255)                                 tgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGAATTGTTATTGACGATACAGGCATCTTCCGccccc                                                                                                                 >  1:1260328/1‑150 (MQ=255)                                  gACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCg                                                                                                                <  1:1100260/150‑1 (MQ=255)                                  gACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCg                                                                                                                >  1:224858/1‑150 (MQ=255)                                    agCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTa                                                                                                              >  1:134392/2‑150 (MQ=255)                                          aCATCGGGAATACGAAACGCGACGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACgg                                                                                                        <  1:816672/150‑1 (MQ=255)                                           cATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGc                                                                                                       <  2:1153554/150‑1 (MQ=255)                                             tCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGcac                                                                                                     >  1:991435/1‑150 (MQ=255)                                                      cGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTgc                                                                                            >  1:1405083/1‑150 (MQ=255)                                                       gAAACGCGCCGCCTATCTCGACCAGTTCGCCGCGAGATACGACCTCGTCTTTTGCGCTGGCAGTGGCCGCGAAGATCAATAACACCGCCGCGGGTTTGTTATTGACGATACTGGCATGTTCCGCCCCGATAATACGGCACAGCAGCTgcg                                                                                           <  2:93851/150‑1 (MQ=255)                                                       gAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACGTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCATAATACGGCACAGCAGCTgcg                                                                                           <  2:98980/150‑1 (MQ=255)                                                       gAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTgcg                                                                                           >  1:1281008/1‑150 (MQ=255)                                                       gAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTacg                                                                                           <  2:936826/150‑1 (MQ=255)                                                        aaaCGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTgcgc                                                                                          >  1:1300725/1‑150 (MQ=255)                                                         aaCGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGcc                                                                                         <  2:253575/150‑1 (MQ=255)                                                           cgcgCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAg                                                                                       >  1:979330/1‑150 (MQ=255)                                                                cgccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACACCAGCTGCGCCAGCGccc                                                                                  >  2:580402/1‑150 (MQ=255)                                                                 gccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCg                                                                                 <  2:1195541/150‑1 (MQ=255)                                                                       ctcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATcgcg                                                                           <  2:979330/150‑1 (MQ=255)                                                                        tcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGa                                                                          <  1:1097805/150‑1 (MQ=255)                                                                        tcCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGa                                                                          <  1:1391044/150‑1 (MQ=255)                                                                         ccaccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGAt                                                                         >  1:616595/1‑150 (MQ=255)                                                                          caccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTTGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATg                                                                        >  1:412898/1‑150 (MQ=255)                                                                                ttCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCcgcg                                                                  <  1:335165/150‑1 (MQ=255)                                                                                     cgcgAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCgg                                                             <  2:1084698/150‑1 (MQ=255)                                                                                      gcgAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGc                                                            >  1:1176723/1‑150 (MQ=255)                                                                                       cgAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcg                                                           >  2:223369/1‑150 (MQ=255)                                                                                          gaTACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgtcg                                                        <  2:384505/150‑1 (MQ=255)                                                                                            tACCACCACCTCTTTTCCGCTGGCACTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGcgacctc                                                      >  2:275357/1‑145 (MQ=255)                                                                                             accaccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                                     >  1:471126/1‑150 (MQ=255)                                                                                             accaccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTcc                                                     >  2:831070/1‑150 (MQ=255)                                                                                              ccaccaccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCa                                                    >  2:320667/1‑150 (MQ=255)                                                                                                 ccaccTCTTTTCCACTGGCAGTGGCCGCCAACATCAATAACACCGCCGCAGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACTGCACAGCAGCTGCGCCAGCTCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGAt                                                 <  2:1021032/150‑1 (MQ=255)                                                                                                   accTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAtt                                                                                                                                          >  2:1432759/1‑59 (MQ=255)                                                                                                   accTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAtt                                                                                                                                          <  1:1432759/59‑1 (MQ=255)                                                                                                   accTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCa                                               <  1:925120/150‑1 (MQ=255)                                                                                                   accTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCa                                               >  2:720317/1‑150 (MQ=255)                                                                                                    ccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCAt                                              <  1:947578/150‑1 (MQ=255)                                                                                                    ccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCAt                                              <  1:223369/150‑1 (MQ=255)                                                                                                    ccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCAt                                              >  2:1376590/1‑150 (MQ=255)                                                                                                    ccTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCAt                                              <  1:550287/150‑1 (MQ=255)                                                                                                          ttCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGa                                        <  2:134392/150‑1 (MQ=255)                                                                                                          ttCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGa                                        >  1:252896/1‑150 (MQ=255)                                                                                                           tCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAg                                       <  2:878994/150‑1 (MQ=255)                                                                                                             cGCTGGCAGTTGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAggg                                     <  2:412898/150‑1 (MQ=255)                                                                                                             cGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAggg                                     >  2:841322/1‑150 (MQ=255)                                                                                                              gCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGt                                    >  1:215783/1‑150 (MQ=255)                                                                                                               cTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTc                                   <  1:720317/150‑1 (MQ=255)                                                                                                                 ggCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                 >  2:116344/1‑150 (MQ=255)                                                                                                                    aGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCa                                                    >  1:615995/1‑128 (MQ=255)                                                                                                                    aGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCa                                                    <  2:615995/128‑1 (MQ=255)                                                                                                                     gTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                             <  1:831070/150‑1 (MQ=255)                                                                                                                      tGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCg                            <  2:302351/150‑1 (MQ=255)                                                                                                                      tGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCg                            <  2:575864/150‑1 (MQ=255)                                                                                                                        gccgccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGaa                          >  2:1103969/1‑150 (MQ=255)                                                                                                                           gccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGc                       <  1:1109371/150‑1 (MQ=255)                                                                                                                           gccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGc                       <  1:900425/150‑1 (MQ=255)                                                                                                                            ccAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                      <  2:1362316/150‑1 (MQ=255)                                                                                                                              aaCATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCTATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATa                    <  2:1260328/150‑1 (MQ=255)                                                                                                                                cgtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGc                  >  2:857674/3‑150 (MQ=255)                                                                                                                                 gtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATa                    >  1:662567/2‑147 (MQ=255)                                                                                                                                 gtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATa                    <  2:662567/146‑1 (MQ=255)                                                                                                                                 gtCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                 >  1:1170999/2‑150 (MQ=255)                                                                                                                                  tCAATAACAATGCGGCGGCGGTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCt                >  1:591610/1‑150 (MQ=255)                                                                                                                                   cAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCGCTGGCGAACGCATAGCCTg               >  2:312775/1‑150 (MQ=255)                                                                                                                                      tAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc            <  2:1176723/150‑1 (MQ=255)                                                                                                                                        acacCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa          >  2:1414273/1‑150 (MQ=255)                                                                                                                                          acCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACg        <  2:1281008/150‑1 (MQ=255)                                                                                                                                              ccgccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCCGCGTCGTCCACATCATACTCGAGGGACACTGGCGAACGCATAGCCTCCGCAACGGCtt    >  1:143187/1‑150 (MQ=255)                                                                                                                                                gccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTcc  <  1:89347/150‑1 (MQ=255)                                                                                                                                                gccgcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTcc  <  1:1103969/150‑1 (MQ=255)                                                                                                                                                gccgcATTGTTATTGACGATACAGGCATCCTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGGCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTcc  <  1:585523/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCC  >  NC_000913/3760643‑3760935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑