Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,658,592	C→T	100%	A175V (GCC→GTC)	ynfE →	putative selenate reductase, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,658,592	0	C	T	96.9%	292.2 / ‑0.9	98	A175V (GCC→GTC)	ynfE	putative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (48/46);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACA  >  NC_000913/1658495‑1658687                                                                                                  |                                                                                                cTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCt                                                                                               >  1:2879463/1‑100 (MQ=255)  tGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGtctc                                                                                              >  2:1627571/1‑100 (MQ=255)     agaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggc                                                                                           <  2:1304498/100‑1 (MQ=255)      gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACTCGCTCTTCGCCATCAGTCTcggcg                                                                                          >  2:185421/1‑100 (MQ=255)      gaaAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTcggcg                                                                                          >  1:102392/1‑100 (MQ=255)       aaaaCCGTCGAACAATATGACAATGAGGCTGTATATAATCAGTACGCTTCGGGGATCGTTGGTGGCACTATGACCCGCTCTTCGCCATCAGCCTc                                                                                              >  2:1548438/1‑95 (MQ=255)        aaaCCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGTGGATCGTTAGCGGCAATATGACCCTCTCTACACCATAAGTCTcgacgg                                                                                         >  1:1316498/1‑99 (MQ=255)            cGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCaaa                                                                                    >  2:2799578/1‑100 (MQ=255)             gTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGTCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGACGGTCAAAc                                                                                   >  1:538060/1‑100 (MQ=255)             gTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAAc                                                                                   <  1:378181/100‑1 (MQ=255)              tCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTTTTCGCCATCAGTCTCGGCGGTCAAACg                                                                                  <  2:2866997/100‑1 (MQ=255)                gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGcc                                                                                >  1:636555/1‑100 (MQ=255)                gAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGcc                                                                                >  2:1954670/1‑100 (MQ=255)                 aaCAGTATGGCTATGATACTGTATCTATTCAGTTCTCTTCGGGGATCGTTAGAGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGcc                                                                                <  2:162672/99‑1 (MQ=255)                  aCAATATGGCAATGAGGCTGTATATATTCAGTACTCTCCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtg                                                                              >  1:819230/1‑100 (MQ=255)                    aaTATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgat                                                                            >  1:2626014/1‑100 (MQ=255)                     atatGGCAATGAGGCTGTATATATTCTGTACTCTTCGTGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCtgatg                                                                           <  1:601329/100‑1 (MQ=255)                       atGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGaa                                                                         >  2:502336/1‑100 (MQ=255)                        tGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAAc                                                                        >  1:2056736/1‑100 (MQ=255)                        tGGCAATCAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTCGCGGCACTATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATCAAc                                                                        >  2:1226235/1‑100 (MQ=255)                           cAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCAGCTCTTCGCCATCAGTCTCAGCGGTCAAACGCATGATGAACTGc                                                                     >  1:247836/1‑100 (MQ=255)                           cAATGAGGCTGTATATAATCAGTACTCTTCGGAGATCGATGGCGGCAATATGACCCGATCTTACCCAACAGTCTCGGCGGTCACACGAATCATGAAATGc                                                                     >  2:238059/1‑100 (MQ=255)                             aTGAGGCTGTATATATTCAGTAGTCTTCGGGGATCGTGGGGGGCAATATAACCCGCTCTTCGCCGTCAGTCTCGGCGGTCAAACGCCTGATGAACt                                                                       <  2:1261238/96‑1 (MQ=255)                              tGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAc                                                                  <  1:819425/100‑1 (MQ=255)                              tGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAc                                                                  >  1:326367/1‑100 (MQ=255)                                 ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                               >  1:747528/1‑100 (MQ=255)                                 ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggc                                                               <  1:2119224/100‑1 (MQ=255)                                 ggCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGGCAAACGCCTGATGAACTGCTAcggc                                                               <  1:185423/100‑1 (MQ=255)                                  gCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcg                                                              >  1:1029656/1‑100 (MQ=255)                                   cTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTAcggcgg                                                             <  1:981095/100‑1 (MQ=255)                                     gTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGtt                                                           >  1:2072507/1‑100 (MQ=255)                                      tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGtt                                                           <  1:2767342/99‑1 (MQ=255)                                      tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  1:1696931/100‑1 (MQ=255)                                      tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  1:2711836/100‑1 (MQ=255)                                      tatataTTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTc                                                          <  1:2560758/100‑1 (MQ=255)                                        tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                        <  1:177263/100‑1 (MQ=255)                                        tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                        <  1:2090392/100‑1 (MQ=255)                                        tatatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCAc                                                        <  2:326364/100‑1 (MQ=255)                                          tatTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTc                                                      <  2:1024190/100‑1 (MQ=255)                                           atTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCa                                                     >  1:2312593/1‑100 (MQ=255)                                            ttCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCaa                                                    >  2:2973929/1‑100 (MQ=255)                                             tCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGTTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAAc                                                   <  1:1448135/100‑1 (MQ=255)                                               aGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCTATCAGTCTTGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCa                                                 >  1:2529870/1‑100 (MQ=255)                                                gTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAg                                                <  1:1653241/100‑1 (MQ=255)                                                 tACTCTCCGGCCATCGTTGTCGGCAATATGACCCGCTCTTCGCCATCAGTCTCAGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGt                                               <  2:2056713/100‑1 (MQ=255)                                                  aCTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTa                                              <  1:1068161/100‑1 (MQ=255)                                                    tctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                            <  1:1077174/100‑1 (MQ=255)                                                    tctTCGGGGATCGTTGGAGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATg                                            <  2:988106/100‑1 (MQ=255)                                                     ctTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATgg                                           <  1:1667397/100‑1 (MQ=255)                                                      ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          <  1:95292/100‑1 (MQ=255)                                                      ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          >  2:1709088/1‑100 (MQ=255)                                                      ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          <  1:1473944/100‑1 (MQ=255)                                                      ttCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGc                                          <  1:2315888/100‑1 (MQ=255)                                                         ggggATCGTTGGCGGCAATATGACCCGCTCTTCGTCATCAGTCTCGGAGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTc                                        <  1:1155980/99‑1 (MQ=255)                                                         ggggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTcc                                       >  1:905787/1‑100 (MQ=255)                                                          gggATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCt                                      <  1:646892/100‑1 (MQ=255)                                                             aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                   <  1:2414606/100‑1 (MQ=255)                                                             aTCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGATCACTCAACCAGTATGGCTCCTACa                                   >  1:929780/1‑100 (MQ=255)                                                             aTCGTGGGCGGCAACATGACCCGCTCTTCGCCATCAGTCTCGGCGTTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACa                                   <  2:2072484/100‑1 (MQ=255)                                                               cGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGc                                 <  2:431655/100‑1 (MQ=255)                                                                gTTGGCGGCAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCa                                >  1:2002012/1‑100 (MQ=255)                                                                   ggcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACt                              >  2:8727/1‑99 (MQ=255)                                                                    gcggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgc                            >  1:1835594/1‑100 (MQ=255)                                                                     cggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcg                           >  1:1453340/1‑100 (MQ=255)                                                                     cggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGACGGATCACTCAACCAGTATGGCTCCTACAGCACTg                             >  1:1099491/1‑98 (MQ=255)                                                                      ggcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGTCGGTTCACTCAACCAGTATGGCTCCTACAGCAATgcg                           >  1:2591043/1‑99 (MQ=255)                                                                       gcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTgcgc                          >  2:1325721/1‑99 (MQ=255)                                                                       gcAATATGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCa                         <  2:2312567/100‑1 (MQ=255)                                                                          ctatTACCCGCTCTTCGCCATTAATCTCGGCGATAAAACGCCTAATGAACTGCTACCGCATTACACTTACCCATTATAGCTCTTACAGAACTATGCaact                      >  2:703101/2‑98 (MQ=255)                                                                          atatGACCCGCTCTTCGCCATCAGACTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCCCCTACAGCACTGCGCAAAt                      >  2:522767/1‑100 (MQ=255)                                                                           tatGACCCGCTCTTCGCCATCAGTCTCGGCGGTTAATCGCCTGATTAATTGCTACGGCTTTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                     <  1:2539547/100‑1 (MQ=255)                                                                           tatGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAAtt                     >  1:589315/1‑100 (MQ=255)                                                                             tGACCCGCTCTTCGCCATCAGTCTCGGCGGTCGAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                   >  1:1265333/1‑100 (MQ=255)                                                                             tGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCCAATTTc                   >  2:1376713/1‑100 (MQ=255)                                                                             tGACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTc                   <  1:45210/100‑1 (MQ=255)                                                                              gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                  >  1:1683784/1‑100 (MQ=255)                                                                              gACCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTcc                  >  2:1402209/1‑100 (MQ=255)                                                                               aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                 >  2:2457092/1‑100 (MQ=255)                                                                               aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                 >  2:2064354/1‑100 (MQ=255)                                                                               aCCCGCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                 >  1:1424182/1‑100 (MQ=255)                                                                               aCCCGCTCTTCGCCATCAGTCACGGCGCTCAAACGCCTGATGAACTACTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTAcc                  >  2:1167333/1‑99 (MQ=255)                                                                                cccGCTCTTCGCCATCAGTCTCGGCGGTCGAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtg                <  2:1265321/100‑1 (MQ=255)                                                                                 ccGCTCTTCGCCATCAGTCTCGGCGGTCAAACGTCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCt                 >  2:1817550/1‑98 (MQ=255)                                                                                  cGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCtgtg              <  2:2459137/100‑1 (MQ=255)                                                                                   gCTCTTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGc             >  2:950677/1‑100 (MQ=255)                                                                                     tctTCGCCGTCAGTCGCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCTTACAGCACTGCGCAAATTTCCTGTGcc            <  2:332510/99‑1 (MQ=255)                                                                                      ctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt          >  2:2386000/1‑100 (MQ=255)                                                                                      ctTCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCAt          >  1:895381/1‑100 (MQ=255)                                                                                      ccactcctccAGTCTCGGCCGTCAAACGCTTGATGAACTGCTACGGCTGTTCACTCAACCATTATGGCTCCTATAGCACTGCGCAAATTTCCTGTGc             <  2:2997698/88‑1 (MQ=255)                                                                                        tCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGc        <  1:421768/100‑1 (MQ=255)                                                                                        tCGCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTAAAGCACTGCGCAAATTTCCTGTGCCATGc        >  2:2428034/1‑100 (MQ=255)                                                                                          gCCATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc      <  1:2511705/100‑1 (MQ=255)                                                                                           ccATCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCt     >  1:1629153/1‑100 (MQ=255)                                                                                           ccATCAGTCTCGGCGGTCAAACGCCTGATAAACTGCCACGGCGGTTCTCTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGccc      <  2:1043270/99‑1 (MQ=255)                                                                                             aTCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTac   <  1:2457121/100‑1 (MQ=255)                                                                                              tCAGTCTCTGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTa    <  1:2683585/98‑1 (MQ=255)                                                                                              tCAGTCTCGGCGGTCAAACGCCTGATGAACTGTAACGTAGGTTGGCTCAACCTGTATGGCGCCTACAGCAGTGCGCTAATTTCCTGTTCCAGGGCCTaca  <  2:1360255/100‑1 (MQ=255)                                                                                              tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCTAATTTCCTGTGCCATGCCCTaca  <  2:628656/100‑1 (MQ=255)                                                                                              tCAGTCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTaca  <  2:1970487/100‑1 (MQ=255)                                                                                                  |                                                                                                CTGAAGAAAACCGTCGAACAATATGGCAATGAGGCTGTATATATTCAGTACTCTTCGGGGATCGTTGGCGGCAATATGACCCGCTCTTCGCCATCAGCCTCGGCGGTCAAACGCCTGATGAACTGCTACGGCGGTTCACTCAACCAGTATGGCTCCTACAGCACTGCGCAAATTTCCTGTGCCATGCCCTACA  >  NC_000913/1658495‑1658687

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