Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,304,838	A→G	intergenic (+445/+270)	eco → / ← mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,304,838	0	A	G	91.5%	227.8 / ‑0.3	95	intergenic (+445/+270)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (6/2);  new base G (41/45);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAT  >  NC_000913/2304741‑2304897                                                                                                  |                                                            gAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCggg                                                           >  2:1602195/1‑100 (MQ=255) gAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCggg                                                           <  2:662707/100‑1 (MQ=255)  aaGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGc                                                          <  1:2087715/100‑1 (MQ=255)  aaGTGCTCCACCCCGTAGGTTGGACAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGc                                                          <  1:1992249/100‑1 (MQ=255)  aaGTGCTCCACCCCGTAGCTTGGATAAGGCGGTGACGCCGCCTCCGACCTCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGc                                                          <  1:1992610/100‑1 (MQ=11)   aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGcc                                                         >  2:502668/1‑100 (MQ=255)   aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGcc                                                         >  2:138148/1‑100 (MQ=255)   aGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         <  1:1955426/100‑1 (MQ=255)   aGTCCTCCACCCCGTACGTCGGATAAGGCGTTTACGCTGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGcc                                                         <  2:2359869/100‑1 (MQ=14)    gTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCt                                                        <  1:220112/100‑1 (MQ=255)    gTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCt                                                        >  1:1344031/1‑100 (MQ=255)    gTGCTCCACCCCCCAGGTTGGATACGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCt                                                        <  2:2000567/100‑1 (MQ=14)     tGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTa                                                       >  1:2571850/1‑100 (MQ=255)     ggCTCCACCTCGTGGGTGGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTa                                                       <  1:1934795/99‑1 (MQ=11)      gCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTac                                                      >  1:1894595/1‑100 (MQ=255)      gCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTac                                                      >  1:1894604/1‑100 (MQ=255)       cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTaca                                                     <  2:781927/100‑1 (MQ=255)       cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCTaca                                                     >  2:838737/1‑100 (MQ=255)        tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTacac                                                    >  1:2377397/1‑100 (MQ=255)        tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                    >  2:761590/1‑100 (MQ=255)         ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAGCGGGCCTACGcc                                                   >  2:1695497/1‑100 (MQ=255)         ccACCCCGTAGGTCGGACAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACAcc                                                   <  2:1487417/100‑1 (MQ=11)           aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGc                                                 >  1:41050/1‑100 (MQ=255)           aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGc                                                 >  2:867702/1‑100 (MQ=255)           aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGc                                                 <  1:1132297/100‑1 (MQ=255)            ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCt                                                <  2:2018090/100‑1 (MQ=255)            ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCt                                                <  2:2017904/100‑1 (MQ=255)            ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCt                                                >  1:2531396/1‑100 (MQ=255)             cccGTCGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGATCCGGTTGTCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtg                                               <  1:2507106/100‑1 (MQ=17)             cccGTAGGTTAGATAAGGCGGCGACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtg                                               <  2:1445787/100‑1 (MQ=11)             cacGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtg                                               <  1:1853588/98‑1 (MQ=255)              ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtgt                                              >  2:1591507/1‑100 (MQ=255)              ccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtgt                                              >  2:356978/1‑100 (MQ=255)               cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCtgtg                                             >  2:2031464/1‑100 (MQ=255)                gTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGa                                            <  2:1955894/100‑1 (MQ=255)                 tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGACCACACCGCTGTGaa                                           >  2:772385/1‑100 (MQ=255)                 tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCcgtggc                                           >  2:772027/1‑98 (MQ=255)                  aGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAg                                          <  1:1717579/100‑1 (MQ=255)                  aGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAg                                          <  2:43812/100‑1 (MQ=255)                   ggTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTCCACCGCTGTGACGt                                         >  2:466767/1‑100 (MQ=255)                   ggTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGt                                         <  2:1574746/100‑1 (MQ=255)                   ggTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGt                                         >  1:917257/1‑100 (MQ=255)                    gTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTg                                        >  1:1907461/1‑100 (MQ=255)                     ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGc                                       <  1:1888001/100‑1 (MQ=255)                      tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCt                                      <  2:1511113/100‑1 (MQ=255)                      tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCt                                      >  2:179817/1‑100 (MQ=255)                      tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCt                                      <  1:754653/100‑1 (MQ=255)                      tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTCCACCGCTGTGCAGTggt                                      >  2:2183086/1‑98 (MQ=11)                       ggATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTc                                     >  2:1225272/1‑100 (MQ=255)                         aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCa                                   >  2:2098283/1‑100 (MQ=255)                          tAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAc                                  >  2:388346/1‑100 (MQ=255)                          tAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAc                                  <  1:518456/100‑1 (MQ=255)                          tAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAc                                  <  2:706490/100‑1 (MQ=255)                          tAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTCGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAc                                  >  2:293739/1‑100 (MQ=255)                            aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAccc                                <  1:2232577/100‑1 (MQ=255)                             ggCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCAcccc                               >  2:1831611/1‑100 (MQ=255)                                ggTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                            <  1:934761/100‑1 (MQ=255)                                ggTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTa                            <  1:2498895/100‑1 (MQ=255)                                 gTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAg                           <  1:1516036/100‑1 (MQ=255)                                 gTTACGCCGCATCCGACAGCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAg                           <  1:2474022/100‑1 (MQ=37)                                  ttACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                          <  1:1570526/100‑1 (MQ=255)                                   tACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTCCACCGCTGTGAAGTGCTCCACCCCGTAGGt                         >  2:479948/1‑100 (MQ=255)                                   tACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                         >  1:524884/1‑100 (MQ=25)                                     cgccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                       >  2:1610917/1‑100 (MQ=21)                                     cgccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                       >  1:744626/1‑100 (MQ=18)                                     cgccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCCGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                       >  2:735937/1‑100 (MQ=18)                                      gccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                      >  2:785521/1‑100 (MQ=18)                                       ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     >  1:2032850/1‑100 (MQ=21)                                       ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     >  1:354675/1‑100 (MQ=21)                                       ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCTCCACCCCGTAGGTCGGa                     >  2:385365/1‑100 (MQ=21)                                         gcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                   <  1:1073546/100‑1 (MQ=18)                                         gcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATa                   >  1:1720526/1‑100 (MQ=21)                                          cATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                  >  1:1115983/1‑100 (MQ=21)                                            tCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                <  1:364171/100‑1 (MQ=21)                                            tCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                <  1:2375462/100‑1 (MQ=21)                                            tCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                <  1:1354376/100‑1 (MQ=18)                                             ncGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTTCACCCCGTAGGTCGGATAAGGc               >  2:2284645/2‑100 (MQ=255)                                             ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               <  1:714641/100‑1 (MQ=21)                                             ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               <  2:1326044/100‑1 (MQ=18)                                             ccGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc               >  2:1652434/1‑100 (MQ=21)                                               gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             >  1:1553787/1‑100 (MQ=18)                                               gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCgg             <  1:87952/100‑1 (MQ=18)                                                aCATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt            <  2:1271587/100‑1 (MQ=18)                                                aCATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGt            >  1:1870632/1‑100 (MQ=21)                                                 cATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt           <  2:2340706/100‑1 (MQ=18)                                                 cATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGtt           <  1:775759/100‑1 (MQ=18)                                                  aTCCAACGTCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa          <  2:708027/100‑1 (MQ=16)                                                   tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         <  1:1256542/100‑1 (MQ=18)                                                   tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         <  1:307639/100‑1 (MQ=18)                                                   tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAc         >  1:1022209/1‑100 (MQ=21)                                                    ccAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACCCCGCTGTGAAGTGCTCTACCCCGTAGGTCGGATAAGGCGGTTAcg        >  2:1802012/1‑100 (MQ=255)                                                    ccAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcg        <  2:152045/100‑1 (MQ=18)                                                    caaaCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcg        <  1:716030/98‑1 (MQ=16)                                                     naaCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGCCCTACACCGCTGTGAAGTGCCCCGCCCCGTAGGTCGGATAAGGCGGTTgcgc       >  2:2224684/2‑100 (MQ=11)                                                        cGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCGCACCCCGCTCTCAGGTGCTACACCCGATGCGTCGGATCAGTCGGTTAcgcggc    >  1:203439/1‑100 (MQ=2)                                                         gTCCGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCa   <  2:1657144/100‑1 (MQ=12)                                                          tCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAt  >  1:1958653/1‑100 (MQ=14)                                                                                                  |                                                            GAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAT  >  NC_000913/2304741‑2304897

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