Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,636,979	A→T	D144E (GAT→GAA)	yhiN ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,636,979	0	A	T	99.0%	333.6 / ‑4.7	95	D144E (GAT→GAA)	yhiN	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (42/52);  total (42/53)
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.68e-01

AATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCAT  >  NC_000913/3636883‑3637074                                                                                                 |                                                                                                aaTCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCgttcagttca                                                                                              <  1:798393/100‑1 (MQ=255)  aTCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCaa                                                                                             <  1:1169534/100‑1 (MQ=255)  aTCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCaa                                                                                             >  2:1457557/1‑100 (MQ=255)   tCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAg                                                                                            <  1:1359503/100‑1 (MQ=255)   tCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAg                                                                                            <  2:892062/100‑1 (MQ=255)   tCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAg                                                                                            >  1:929/1‑100 (MQ=255)       ataACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTg                                                                                        <  2:1046181/100‑1 (MQ=255)       ataACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTg                                                                                        >  1:1266942/1‑100 (MQ=255)          aCCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAg                                                                                     >  1:1927389/1‑100 (MQ=255)            cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGcc                                                                                   >  2:1914794/1‑100 (MQ=255)            cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGcc                                                                                   >  1:1829850/1‑100 (MQ=255)            cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGcc                                                                                   <  2:285600/100‑1 (MQ=255)            cAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGcc                                                                                   >  1:1765479/1‑100 (MQ=255)              aaCGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTg                                                                                 <  2:1141021/100‑1 (MQ=255)              aaCGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTg                                                                                 >  1:2171806/1‑100 (MQ=255)               aCGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGt                                                                                >  1:1286329/1‑100 (MQ=255)                cGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGtt                                                                               >  2:2270708/1‑100 (MQ=255)                 ggCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGttt                                                                              <  1:1062588/100‑1 (MQ=255)                  gCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTc                                                                             >  2:241/1‑100 (MQ=255)                  gCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTc                                                                             <  2:1879840/100‑1 (MQ=255)                    gACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCAt                                                                           <  2:2322283/100‑1 (MQ=255)                    gACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCAt                                                                           >  2:110495/1‑100 (MQ=255)                     aCGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATc                                                                          <  1:1914801/100‑1 (MQ=255)                      cgcgCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATcc                                                                         >  1:75881/1‑100 (MQ=255)                      cgcgCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATcc                                                                         <  1:938282/100‑1 (MQ=255)                      cgcgCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATcc                                                                         >  1:973707/1‑100 (MQ=255)                        cgcCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCtt                                                                       >  1:1526477/1‑100 (MQ=255)                            cAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGcc                                                                   <  2:780248/100‑1 (MQ=255)                              gCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCcac                                                                 <  2:1390509/100‑1 (MQ=255)                                cccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCacac                                                               >  1:34746/1‑100 (MQ=255)                                cccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCacac                                                               <  1:399194/100‑1 (MQ=255)                                cccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCacac                                                               <  2:629955/100‑1 (MQ=255)                                 ccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACt                                                              >  1:605551/1‑100 (MQ=255)                                 ccGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACt                                                              <  2:1852392/100‑1 (MQ=255)                                   ggCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCa                                                            <  2:1239030/100‑1 (MQ=255)                                     cATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGc                                                          <  2:674792/100‑1 (MQ=255)                                      aTTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCa                                                         >  2:290430/1‑100 (MQ=255)                                       ttGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCAc                                                        >  1:2062667/1‑100 (MQ=255)                                        tGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACt                                                       >  2:2050063/1‑100 (MQ=255)                                        tGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACt                                                       >  2:2250190/1‑100 (MQ=255)                                         gACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACtt                                                      >  1:194450/1‑100 (MQ=255)                                          aCAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTc                                                     <  1:1753239/100‑1 (MQ=255)                                            aGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCTACACTCAGCACTTCGc                                                   <  1:853103/100‑1 (MQ=255)                                             gCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCt                                                  >  2:1875159/1‑100 (MQ=255)                                             gCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCt                                                  <  1:772645/100‑1 (MQ=255)                                             gCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCt                                                  <  1:314516/100‑1 (MQ=255)                                              cccACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTa                                                 >  1:730243/1‑100 (MQ=255)                                              cccACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTa                                                 >  1:1246826/1‑100 (MQ=255)                                                caccaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACg                                               >  2:2152153/1‑100 (MQ=255)                                                  ccaGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCa                                             <  1:457290/100‑1 (MQ=255)                                                    aGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAAt                                           >  1:278618/1‑100 (MQ=255)                                                     gTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATc                                          <  2:1927382/100‑1 (MQ=255)                                                      tCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCt                                         <  1:2152160/100‑1 (MQ=255)                                                      tCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCt                                         <  2:1725929/100‑1 (MQ=255)                                                      tCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCt                                         <  2:973705/100‑1 (MQ=255)                                                       cgcgATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTg                                        <  1:1777444/100‑1 (MQ=255)                                                        gcgATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGa                                       >  2:1030046/1‑100 (MQ=255)                                                        gcgATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGa                                       <  1:241/100‑1 (MQ=255)                                                         cgATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGaa                                      >  1:1441480/1‑100 (MQ=255)                                                         cgATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGaa                                      <  1:290432/100‑1 (MQ=255)                                                           aTGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAgg                                    >  2:1873737/1‑100 (MQ=255)                                                             gACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTc                                  <  2:2187562/100‑1 (MQ=255)                                                              aCCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTca                                 <  2:1246823/100‑1 (MQ=255)                                                              aCCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTca                                 <  2:730242/100‑1 (MQ=255)                                                               ccAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTcac                                >  2:1478865/1‑100 (MQ=255)                                                               ccAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTcac                                >  1:673396/1‑100 (MQ=255)                                                               ccAGCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTcac                                >  1:1773695/1‑100 (MQ=255)                                                                  gCTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACAtt                             >  1:647387/1‑100 (MQ=255)                                                                   cTTTTCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTg                            >  2:2221833/1‑100 (MQ=255)                                                                      ttCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGccc                         >  2:1604289/1‑100 (MQ=255)                                                                      ttCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGccc                         <  1:2190887/100‑1 (MQ=255)                                                                       tCGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCt                        >  2:941856/1‑100 (MQ=255)                                                                        cGCAACCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCtt                       >  2:285824/1‑100 (MQ=255)                                                                            aCCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTtctca                   <  2:1125202/100‑2 (MQ=255)                                                                            aCCGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCg                   <  1:1478869/100‑1 (MQ=255)                                                                             ccGACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGc                  <  1:110497/100‑1 (MQ=255)                                                                               gACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCAc                <  1:1901055/100‑1 (MQ=255)                                                                               gACAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCAc                <  2:658372/100‑1 (MQ=255)                                                                                aCAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACt               <  1:66362/100‑1 (MQ=255)                                                                                aCAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCGTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACt               <  1:1875165/100‑1 (MQ=255)                                                                                 cAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTc              <  1:640391/100‑1 (MQ=255)                                                                                 cAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTc              <  2:1240948/100‑1 (MQ=255)                                                                                 cAGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTc              >  2:705479/1‑100 (MQ=255)                                                                                  aGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCGCTTCTCGCACTCa             <  2:1441476/100‑1 (MQ=255)                                                                                  aGTCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCa             <  2:2278116/100‑1 (MQ=255)                                                                                    tCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCa             <  2:2062660/98‑1 (MQ=255)                                                                                    tCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATc           <  1:1457561/100‑1 (MQ=255)                                                                                    tCATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATc           <  2:34745/100‑1 (MQ=255)                                                                                     cATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATcc          <  2:1266939/100‑1 (MQ=255)                                                                                     cATGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATcc          >  2:1616131/1‑100 (MQ=255)                                                                                      aTGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATcca         >  2:1668391/1‑100 (MQ=255)                                                                                      aTGCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATcca         <  2:2141351/100‑1 (MQ=255)                                                                                        gCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccacc       <  2:2289526/100‑1 (MQ=255)                                                                                        gCCGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccacc       >  2:1258500/1‑100 (MQ=255)                                                                                         ccGTTCAGTTCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATccacca      <  1:1604294/100‑1 (MQ=255)                                                                                           gttcagttcaAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGc    >  1:1523991/1‑100 (MQ=255)                                                                                             tCAGTTCAAGCGTGCAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCAt  <  1:2221841/100‑1 (MQ=255)                                                                                                 |                                                                                                AATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCGATGACCAGCTTTTCGCAACCGACAGTCATGCCGTTCAGATCAAGCGTGAAGCCTGTTTCATCCTTCGCCACACTCAGCACTTCGCTACGCAATCTGAAGGTCACATTGCCCTTCTCGCACTCATCCACCAGCAT  >  NC_000913/3636883‑3637074

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