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	100.0%	273.5 / NA	84	*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);  new base A (41/43);  total (41/43)

CTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTG  >  NC_000913/3939089‑3939345                                                                                                                                    |                                                                                                                               cTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  1:45000/139‑1 (MQ=255) cTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  2:216741/139‑1 (MQ=255)  tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  2:842590/1‑139 (MQ=255)  caTCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  2:821190/2‑139 (MQ=255)    tCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  1:474561/135‑1 (MQ=255)    tCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  2:474561/1‑135 (MQ=255)       accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCtt                                                                                                                  >  2:385195/1‑139 (MQ=255)       accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCGTAAATTTACGCTGTCtt                                                                                                                  <  1:821190/139‑1 (MQ=255)            ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             <  1:569742/139‑1 (MQ=255)                            gggAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCg                                                                                             >  1:437807/1‑139 (MQ=255)                             ggAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                            <  2:1240384/139‑1 (MQ=255)                               aGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgtt                                                                                          <  1:385195/139‑1 (MQ=255)                                gCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCttt                                                                                                                 >  2:366961/1‑115 (MQ=255)                                gCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCttt                                                                                                                 <  1:366961/115‑1 (MQ=255)                                gCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             >  1:1222930/1‑119 (MQ=255)                                gCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             <  2:1222930/119‑1 (MQ=255)                                gCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgttg                                                                                         >  2:345070/1‑139 (MQ=255)                                   ggCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGttt                                                                                      >  1:1138579/1‑139 (MQ=255)                                      gATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   <  1:705369/139‑1 (MQ=255)                                        tggatgagtacatcaatcgGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAatta                                                                                                      >  1:780509/11‑119 (MQ=255)                                        tggatgagtacatcaatcgGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAatta                                                                                                      <  2:780509/109‑1 (MQ=255)                                                cTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGaa                                                                         >  1:261325/1‑139 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         <  2:1073462/88‑1 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                         >  1:1073462/1‑88 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       <  1:586457/90‑1 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  2:586457/1‑90 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             <  2:109462/100‑1 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                             >  1:109462/1‑100 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACAt                                                                      >  2:1282074/1‑139 (MQ=255)                                                   atccatc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACAt                                                                      >  2:1019830/1‑139 (MQ=255)                                                      catc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  2:985212/86‑1 (MQ=255)                                                      catc‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        >  1:985212/1‑86 (MQ=255)                                                         c‑GGATAACCTAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                                <  1:842590/139‑1 (MQ=255)                                                           ggATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTgctgct                                                               >  1:405806/1‑139 (MQ=255)                                                              tAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  1:958119/79‑1 (MQ=255)                                                              tAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        >  2:958119/1‑79 (MQ=255)                                                                aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          <  1:1019830/139‑1 (MQ=255)                                                                aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          >  1:958704/1‑139 (MQ=255)                                                                    aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACa                                                      >  1:574872/1‑139 (MQ=255)                                                                      ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          <  1:1078463/133‑1 (MQ=255)                                                                      ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                          >  2:1078463/1‑133 (MQ=255)                                                                         aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                 <  1:1037024/139‑1 (MQ=255)                                                                         aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                 >  2:1198380/1‑139 (MQ=255)                                                                          cccTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   >  2:464773/1‑104 (MQ=255)                                                                          cccTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                   <  1:464773/104‑1 (MQ=255)                                                                                aGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGccc                                          <  2:574872/139‑1 (MQ=255)                                                                                aGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGccc                                          >  1:280759/1‑139 (MQ=255)                                                                                 gCAATAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAATCAACAGTAATAATGCCCa                                         <  1:865696/139‑1 (MQ=255)                                                                                 gCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCa                                         <  2:261325/139‑1 (MQ=255)                                                                                 gCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCa                                         <  1:1282074/139‑1 (MQ=255)                                                                                  caacaaCGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAt                                        <  2:1138579/139‑1 (MQ=255)                                                                                  caacaaCGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAt                                        <  2:461261/139‑1 (MQ=255)                                                                                     caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATtgt                                     >  1:408362/1‑139 (MQ=255)                                                                                       acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        >  2:694005/1‑54 (MQ=255)                                                                                       acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                        <  1:694005/54‑1 (MQ=255)                                                                                       acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGa                                   <  2:405806/139‑1 (MQ=255)                                                                                        cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       <  1:731558/54‑1 (MQ=255)                                                                                        cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                       >  2:731558/1‑54 (MQ=255)                                                                                        cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCa                                         <  2:1262865/132‑1 (MQ=255)                                                                                        cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCa                                         >  1:1262865/1‑132 (MQ=255)                                                                                          aTTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATaa                                <  1:1063910/139‑1 (MQ=255)                                                                                          aTTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATaa                                >  2:1063910/1‑139 (MQ=255)                                                                                                 cTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAt                                        >  2:26237/1‑124 (MQ=255)                                                                                                 cTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAt                                        <  1:26237/124‑1 (MQ=255)                                                                                                  ttACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAg                        <  2:437807/139‑1 (MQ=255)                                                                                                      tCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCAGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGa                    >  2:894612/1‑139 (MQ=255)                                                                                                           ttCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAt                                        <  1:1144304/114‑1 (MQ=255)                                                                                                           ttCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAt                                        >  2:1144304/1‑114 (MQ=255)                                                                                                           ttCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGt               >  2:946238/1‑139 (MQ=255)                                                                                                                aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGaa                                                            >  2:1138770/1‑89 (MQ=255)                                                                                                                aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGaa                                                            <  1:1138770/89‑1 (MQ=255)                                                                                                                aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGAc          <  1:946238/139‑1 (MQ=255)                                                                                                                aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGAc          <  1:114169/139‑1 (MQ=255)                                                                                                                    aaCGCGGTTCGGCTTAAACTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGAAAGTc      >  1:865642/1‑139 (MQ=255)                                                                                                                     aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATa                                 <  2:474130/111‑1 (MQ=255)                                                                                                                     aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATa                                 >  1:474130/1‑111 (MQ=255)                                                                                                                      cgcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgtt                                                                                          <  1:192705/53‑1 (MQ=255)                                                                                                                      cgcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgtt                                                                                          >  2:192705/1‑53 (MQ=255)                                                                                                                      cgcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTct                                                                                   >  2:569825/1‑59 (MQ=255)                                                                                                                      cgcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTct                                                                                   <  1:569825/60‑2 (MQ=255)                                                                                                                       gcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAtt                                       <  1:1091374/103‑1 (MQ=255)                                                                                                                       gcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCAtt                                       >  2:1091374/1‑103 (MQ=255)                                                                                                                       gcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGtt   <  1:1119070/139‑1 (MQ=255)                                                                                                                        cgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTtg  >  2:748303/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                               CTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATC‑GGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTG  >  NC_000913/3939089‑3939345

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