Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	95.4%	306.8 / 8.9	87	coding (10/1353 nt)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (47/36);  total (49/38)
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.96e-01

GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281                                                                                                     |                                                                                                 gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                   >  2:775621/1‑100 (MQ=255) gCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAt                                                                                                   <  1:2524277/100‑1 (MQ=255)  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                 >  1:522219/1‑100 (MQ=255)  ccTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑c                                                                                                 >  1:102600/1‑100 (MQ=255)   cTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc                                                                                                >  1:2717901/1‑100 (MQ=255)    tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑c                                                                                              <  2:543687/100‑1 (MQ=255)      gCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccc                                                                                            >  1:2368649/1‑100 (MQ=255)       cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑cccc                                                                                           >  2:1067856/1‑100 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccc                                                                                          <  2:1214692/100‑2 (MQ=255)        ttCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccc                                                                                          <  2:468492/100‑2 (MQ=255)          cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccctg                                                                                        <  1:775626/100‑4 (MQ=255)          cGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑cc‑ccccctg                                                                                        >  1:2615153/1‑97 (MQ=255)             cTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTc                                                                                     >  1:1030559/1‑100 (MQ=255)                 aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGa                                                                                 <  2:1561824/100‑1 (MQ=255)                     tCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAAt                                                                             <  2:2474712/100‑1 (MQ=255)                      cACTTCAAATCGCACCAAATTATTTCCACAGAGTCTTCGTATCCTTCTACCACTGAAAATGAAAAGAATATGTAGAAAT‑CC‑CCCCATGTTCTGGAAAAt                                                                             <  2:1376803/99‑1 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGAc                                                                          >  2:1563644/1‑100 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGAc                                                                          <  1:1546735/100‑1 (MQ=255)                        cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCCGGAAAATGAc                                                                          <  1:2397704/100‑1 (MQ=255)                         ttCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCC‑TGTTCTGGAAAATGACTa                                                                        >  2:1687786/1‑100 (MQ=255)                          tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACt                                                                         >  2:516429/1‑99 (MQ=255)                          tCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGGTGAAAT‑CCGCCCC‑TGTTCTGGAAAATGACTa                                                                        >  2:1859151/1‑100 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAg                                                                       <  2:1694825/100‑1 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAg                                                                       >  1:2181305/1‑100 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAg                                                                       >  2:571933/1‑100 (MQ=255)                           cacaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCATGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTa                                                                        >  2:1697992/1‑99 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGcc                                                                     >  1:2431519/1‑100 (MQ=255)                             caTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGcc                                                                     <  1:1687222/100‑1 (MQ=255)                                 gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGtt                                                                 >  1:1385062/1‑100 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:914722/1‑100 (MQ=255)                                     cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGt                                                             >  2:1350981/1‑100 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  2:176123/100‑1 (MQ=255)                                      aaaTTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTg                                                            <  1:2098511/100‑1 (MQ=255)                                         ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  2:2542882/100‑1 (MQ=255)                                         ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAg                                                         <  1:83841/100‑1 (MQ=255)                                           cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTa                                                       >  2:1951991/1‑100 (MQ=255)                                              cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                    <  2:568900/100‑1 (MQ=255)                                               ccACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   <  2:1551323/100‑1 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                   >  1:1538294/1‑98 (MQ=255)                                                 acaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                 >  2:2287222/1‑100 (MQ=255)                                                  caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                >  2:1571115/1‑100 (MQ=255)                                                  caTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgt                                                >  2:2057596/1‑100 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              >  1:997471/1‑100 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              >  1:664842/1‑100 (MQ=255)                                                    tAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttg                                              <  1:2088379/100‑1 (MQ=255)                                                     aGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtgttgt                                             <  1:785545/100‑1 (MQ=255)                                                        cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTgttgtt                                            <  2:1910262/98‑1 (MQ=255)                                                        cTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGa                                          >  1:1058647/1‑100 (MQ=255)                                                         ttCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                         <  2:2724342/100‑1 (MQ=255)                                                           cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  1:2005532/100‑1 (MQ=255)                                                           cGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTc                                       <  2:2717882/100‑1 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      <  2:1934399/100‑1 (MQ=255)                                                            gTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTcc                                      >  1:2109237/1‑100 (MQ=255)                                                             tATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCg                                     >  1:943871/1‑100 (MQ=255)                                                              aTCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGa                                    >  2:174576/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:2743134/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  1:1876606/1‑100 (MQ=255)                                                               tCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGAt                                   >  2:512167/1‑100 (MQ=255)                                                                ccTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATa                                  >  2:570644/1‑100 (MQ=255)                                                                  tgctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAAt                                >  1:780709/1‑100 (MQ=255)                                                                    ctgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                              >  2:76846/1‑100 (MQ=255)                                                                     tgcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCt                             <  2:2208440/100‑1 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            <  2:1182276/100‑1 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:299844/1‑100 (MQ=255)                                                                      gcCATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTg                            >  2:299637/1‑100 (MQ=255)                                                                       ccATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGa                           <  2:2615134/100‑1 (MQ=255)                                                                        cATTGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAt                          >  2:658680/1‑100 (MQ=255)                                                                          ttGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTc                        <  2:1043713/100‑1 (MQ=255)                                                                           tGCAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                       <  2:522215/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGGTTCGGc                     <  2:2245536/100‑1 (MQ=255)                                                                             cAAAGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGc                     <  1:329571/100‑1 (MQ=255)                                                                              aaaGGAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCa                    <  1:659227/100‑1 (MQ=255)                                                                                 ggAGAAGACTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGt                 <  1:2672706/100‑1 (MQ=255)                                                                                   agaagaCTATGTTGAAAT‑CC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  1:1687799/100‑1 (MQ=255)                                                                                    gaagaCTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt              <  1:2098733/100‑1 (MQ=255)                                                                                      ccgatctGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  2:828580/94‑1 (MQ=255)                                                                                      agaCTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            <  2:1385053/100‑1 (MQ=255)                                                                                      agaCTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTg            >  1:2099297/1‑100 (MQ=255)                                                                                       gaCTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt           >  2:985420/1‑100 (MQ=255)                                                                                        aCTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc          <  1:2075939/100‑1 (MQ=255)                                                                                         cTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         >  1:451424/1‑100 (MQ=255)                                                                                         cTATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCg         <  1:1743569/100‑1 (MQ=255)                                                                                          tATGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGa        <  1:1067863/100‑1 (MQ=255)                                                                                            tGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  2:956629/1‑100 (MQ=255)                                                                                            tGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:928294/1‑100 (MQ=255)                                                                                            tGTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc      >  1:111766/1‑100 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTgagatcgg     >  1:828586/1‑93 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  1:810156/1‑100 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2151431/1‑100 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:2038534/1‑100 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     <  1:2743153/100‑1 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1825898/1‑100 (MQ=255)                                                                                             gTTGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACg     >  2:1366369/1‑100 (MQ=255)                                                                                              ttGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  2:592636/1‑100 (MQ=255)                                                                                              ttGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2572712/1‑100 (MQ=255)                                                                                              ttGAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt    >  1:2572540/1‑100 (MQ=255)                                                                                                gAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:2690438/1‑100 (MQ=255)                                                                                                gAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  >  2:1987396/1‑100 (MQ=255)                                                                                                gAAATCCC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGc  <  2:102598/100‑1 (MQ=255)                                                                                                     |                                                                                                 GCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAAGACTATGTTGAAAT‑CC‑CCCC‑TGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGC  >  NC_000913/4638088‑4638281

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