Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,939,219	G→A	*331* (TAG→TAA)  S465F (TCT→TTT)	rbsR → hsrA ←	transcriptional repressor of ribose metabolism putative multidrug or homocysteine efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,939,219	0	G	A	99.0%	310.6 / NA	95	*331* (TAG→TAA)  S465F (TCT→TTT)	rbsR hsrA	transcriptional repressor of ribose metabolism putative multidrug or homocysteine efflux system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (45/49);  minor base C (1/0);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTG  >  NC_000913/3939090‑3939345                                                                                                                                  |                                                                                                                               tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  2:11784/137‑1 (MQ=255) tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  1:11784/1‑137 (MQ=255) tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  1:1221788/1‑139 (MQ=255)     caccaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCt                                                                                                                   >  1:487468/1‑139 (MQ=255)      accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  2:1406556/1‑132 (MQ=255)      accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  1:1406556/132‑1 (MQ=255)        caACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAGCGCGGTTCGGCTTAAATTTACGCTGTCtttt                                                                                                                >  1:1186686/1‑139 (MQ=255)          aCCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGa                                                                                                              >  1:1368794/1‑139 (MQ=255)           ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  2:900296/1‑129 (MQ=255)           ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       <  1:900296/129‑1 (MQ=255)           ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             >  1:651759/1‑139 (MQ=255)           ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             <  2:487468/139‑1 (MQ=255)           ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             <  1:931851/139‑1 (MQ=255)           ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             >  2:931851/1‑139 (MQ=255)            cGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATc                                                                                                            <  1:290947/139‑1 (MQ=255)               aaGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCaaa                                                                                                         >  2:17823/1‑139 (MQ=255)                     aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAAttatta                                                                                                   >  1:1610439/1‑139 (MQ=255)                     aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTT‑CTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTAc                                                                                                  >  2:1389745/1‑139 (MQ=255)                      aCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTAc                                                                                                  >  1:1276524/1‑139 (MQ=255)                        tGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCa                                                                                                >  1:1583176/1‑139 (MQ=255)                           gggAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTT‑CTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                            <  1:1389745/139‑1 (MQ=255)                               gCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgttg                                                                                         >  1:1375889/1‑139 (MQ=255)                                cTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgttgt                                                                                        >  1:1076908/1‑139 (MQ=255)                                     gATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   <  2:651759/139‑1 (MQ=255)                                       ttGATGTACTCATCCATCGGATAACCTAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                 <  2:1076908/139‑1 (MQ=255)                                       ttGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAGCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                 <  2:1186686/139‑1 (MQ=255)                                          aTGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                              <  2:1566496/139‑1 (MQ=255)                                             tACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  1:1663347/94‑1 (MQ=255)                                             tACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        >  2:1663347/1‑94 (MQ=255)                                                 catccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             >  1:574255/1‑101 (MQ=255)                                                 catccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             <  2:574255/101‑1 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTa                                                                                                                           <  1:436940/86‑1 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTa                                                                                                                           >  2:436940/1‑86 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        >  2:171403/1‑89 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  1:171403/89‑1 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       <  1:254409/90‑1 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  2:254409/1‑90 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGCTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   <  1:542703/126‑1 (MQ=255)                                                  atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGCTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   >  2:542703/1‑126 (MQ=255)                                                   tccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  1:530709/1‑89 (MQ=255)                                                   tccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       <  2:530709/89‑1 (MQ=255)                                                    ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       <  2:972127/88‑1 (MQ=255)                                                    ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  1:972127/1‑88 (MQ=255)                                                    ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                    <  1:800958/139‑1 (MQ=255)                                                           aTAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGa                                                             <  1:17823/139‑1 (MQ=255)                                                              aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  1:887019/76‑1 (MQ=255)                                                              aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  2:887019/1‑76 (MQ=255)                                                              aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          >  1:1613969/1‑139 (MQ=255)                                                               cccAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCa                                                         <  2:1375889/139‑1 (MQ=255)                                                                 cAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                       <  2:1613969/139‑1 (MQ=255)                                                                  aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  2:202815/1‑72 (MQ=255)                                                                  aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  1:202815/72‑1 (MQ=255)                                                                  aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACa                                                      <  2:1583176/139‑1 (MQ=255)                                                                   gCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                     <  2:1221788/139‑1 (MQ=255)                                                                    ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  1:327864/70‑1 (MQ=255)                                                                    ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  2:327864/1‑70 (MQ=255)                                                                    ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAt                                                                                               <  1:335754/96‑1 (MQ=255)                                                                    ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAt                                                                                               >  2:335754/1‑96 (MQ=255)                                                                       aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          >  2:1040360/1‑130 (MQ=255)                                                                       aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          <  1:1040360/130‑1 (MQ=255)                                                                       aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                 <  2:1276524/139‑1 (MQ=255)                                                                       aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                 <  1:1666714/139‑1 (MQ=255)                                                                        cccTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   <  1:1217601/104‑1 (MQ=255)                                                                        cccTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   >  2:1217601/1‑104 (MQ=255)                                                                          cTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTaataat                                              <  1:1403718/139‑1 (MQ=255)                                                                            tCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGc                                            >  2:971564/1‑139 (MQ=255)                                                                             cAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   >  1:350175/1‑99 (MQ=255)                                                                             cAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   <  2:350175/99‑1 (MQ=255)                                                                                 aacaacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAtt                                       <  2:1368794/139‑1 (MQ=255)                                                                                     acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  1:1070245/53‑1 (MQ=255)                                                                                     acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  2:1070245/1‑53 (MQ=255)                                                                                      cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                              <  1:7146/95‑1 (MQ=255)                                                                                      cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                              >  2:7146/1‑95 (MQ=255)                                                                                       gATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATa                                 <  2:301108/139‑1 (MQ=255)                                                                                                   cTCCGATTCTGATGGAACGCGGTTCGGCTTACATTTACGCTGTCttt                                                                                                                 >  2:1459739/1‑47 (MQ=255)                                                                                                   cTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCttt                                                                                                                 <  1:1459739/47‑1 (MQ=255)                                                                                                    tCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAt                                                                                               >  2:1105353/1‑64 (MQ=255)                                                                                                    tCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAt                                                                                               <  1:1105353/64‑1 (MQ=255)                                                                                                         ttgtGATGGAATGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                                >  2:980047/4‑90 (MQ=255)                                                                                                         ttCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                                <  1:980047/90‑1 (MQ=255)                                                                                                              aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                                <  1:960141/85‑1 (MQ=255)                                                                                                              aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                                >  2:960141/1‑85 (MQ=255)                                                                                                                  aaCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          >  2:709572/1‑87 (MQ=255)                                                                                                                  aaCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          <  1:709572/87‑1 (MQ=255)                                                                                                                   aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                              <  1:913615/66‑1 (MQ=255)                                                                                                                   aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                              >  2:913615/1‑66 (MQ=255)                                                                                                                   aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                     >  2:36654/1‑91 (MQ=255)                                                                                                                   aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                     <  1:36654/91‑1 (MQ=255)                                                                                                                      cgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAt                                                                                               <  2:987914/46‑1 (MQ=255)                                                                                                                      cgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAt                                                                                               >  1:987914/1‑46 (MQ=255)                                                                                                                      cgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          <  2:712846/83‑1 (MQ=255)                                                                                                                      cgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          >  1:712846/1‑83 (MQ=255)                                                                                                                      cgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTtg  >  2:1308236/1‑139 (MQ=255)                                                                                                                           cGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATg                                             <  2:570861/91‑1 (MQ=255)                                                                                                                           cGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATg                                             >  1:570861/1‑91 (MQ=255)                                                                                                                                  |                                                                                                                               TATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTG  >  NC_000913/3939090‑3939345

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