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_0009133,939,2190GA100.0% 314.4 / NA 96*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 (44/52);  total (44/52)

ACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTT  >  NC_000913/3939088‑3939351
                                                                                                                                   |                                                                                                                                    
aCTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  2:1005681/139‑1 (MQ=255)
 cTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACCACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              <  2:225857/139‑1 (MQ=255)
 cTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              >  1:275629/1‑139 (MQ=255)
 cTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              >  1:1581397/1‑139 (MQ=255)
 cTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              <  2:812579/139‑1 (MQ=255)
  tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                             <  2:1264815/139‑1 (MQ=255)
  tATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                             <  1:304782/139‑1 (MQ=255)
       accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCtt                                                                                                                        >  1:280382/1‑139 (MQ=255)
         caACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCtttt                                                                                                                      >  2:1583261/1‑139 (MQ=255)
         caACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACA‑‑‑‑ACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATc                                                                                                                  <  1:950143/139‑1 (MQ=255)
            ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                                   >  2:876581/1‑139 (MQ=255)
                   aTGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAatta                                                                                                            <  2:275629/139‑1 (MQ=255)
                   aTGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAatta                                                                                                            >  1:653541/1‑139 (MQ=255)
                      aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCa                                                                                                                 >  1:765741/1‑131 (MQ=255)
                      aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCa                                                                                                                 <  2:765741/131‑1 (MQ=255)
                      aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAAttatta                                                                                                         >  1:1280362/1‑139 (MQ=255)
                      aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAAttatta                                                                                                         >  1:905514/1‑139 (MQ=255)
                      aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAAttatta                                                                                                         >  1:1578459/1‑139 (MQ=255)
                            gggAGCTGGCGATTGATGTACTCATCCATCGGATAACCCACCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCg                                                                                                   >  1:1307930/1‑139 (MQ=255)
                             ggAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTTCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                                  <  1:352893/139‑1 (MQ=255)
                             ggAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                                  <  2:1307930/139‑1 (MQ=255)
                             ggAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACACTTCAGCAACAACGATTTCAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCAACgg                                                                                                  <  1:1036370/139‑1 (MQ=255)
                                gCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                                  <  2:463267/136‑1 (MQ=255)
                                gCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                                  >  1:463267/1‑136 (MQ=255)
                                      gATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                         <  1:183128/139‑1 (MQ=255)
                                      gATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                         <  1:1332841/139‑1 (MQ=255)
                                        ttGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                       <  1:1206414/139‑1 (MQ=255)
                                           aTGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                                    <  1:876581/139‑1 (MQ=255)
                                           aTGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                                    <  2:1578459/139‑1 (MQ=255)
                                                   atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                                  <  1:529221/117‑1 (MQ=255)
                                                   atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                                  >  2:529221/1‑117 (MQ=255)
                                                   atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACAt                                                                            >  2:667228/1‑139 (MQ=255)
                                                   atccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACAt                                                                            >  2:1343665/1‑139 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  2:73279/86‑1 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               >  1:73279/1‑86 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               >  2:1090944/1‑86 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  1:1090944/86‑1 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                          <  1:763049/139‑1 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                          >  2:938712/1‑139 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                          <  1:1583261/139‑1 (MQ=255)
                                                      catcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                         <  2:367838/123‑1 (MQ=255)
                                                      catcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                         >  1:367838/1‑123 (MQ=255)
                                                        tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              >  2:204974/1‑84 (MQ=255)
                                                        tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              <  1:204974/84‑1 (MQ=255)
                                                        tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctg                                                                       >  1:1333751/1‑139 (MQ=255)
                                                        tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctg                                                                       >  2:826036/1‑139 (MQ=255)
                                                         cGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                                      >  1:526259/1‑139 (MQ=255)
                                                              aaCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAg                                                                 <  2:320890/139‑1 (MQ=255)
                                                               aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                                <  1:1248743/139‑1 (MQ=255)
                                                               aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                                <  2:905514/139‑1 (MQ=255)
                                                                cccAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCa                                                               >  2:653256/1‑139 (MQ=255)
                                                                 ccAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCaa                                                              <  1:1343665/139‑1 (MQ=255)
                                                                  cAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                             <  2:1372242/75‑1 (MQ=255)
                                                                  cAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                             >  1:1372242/1‑75 (MQ=255)
                                                                   aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                          <  1:12476/125‑1 (MQ=255)
                                                                   aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                          >  2:12476/1‑125 (MQ=255)
                                                                   aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                             <  2:889389/138‑1 (MQ=255)
                                                                   aGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                             >  1:889389/1‑138 (MQ=255)
                                                                    gCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                           <  2:1280362/139‑1 (MQ=255)
                                                                    gCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                           <  2:526259/139‑1 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGaa                                                                               >  1:185559/1‑118 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGaa                                                                               <  2:185559/118‑1 (MQ=255)
                                                                       gACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaa                                                        >  2:44589/1‑139 (MQ=255)
                                                                        aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                       <  2:653541/139‑1 (MQ=255)
                                                                           cTTCAGCAACAACGATTACAAGTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCGTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTaataat                                                    <  1:938712/139‑1 (MQ=255)
                                                                            ttCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATg                                                   <  1:1513952/139‑1 (MQ=255)
                                                                               aGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTatg                                                                                                                               >  1:382118/1‑58 (MQ=255)
                                                                               aGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               >  1:139354/1‑60 (MQ=255)
                                                                               aGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  2:139354/60‑1 (MQ=255)
                                                                               aGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  2:382118/60‑1 (MQ=255)
                                                                                gCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCa                                               <  1:667228/139‑1 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                                    <  1:1568320/98‑1 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGc                                                                                    >  2:1568320/1‑98 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATtgt                                           <  1:653256/139‑1 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCGGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGTTGCTGAAGCAACAGTAATAATGCCCATtgt                                           >  1:627486/1‑139 (MQ=255)
                                                                                      acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  2:152849/53‑1 (MQ=255)
                                                                                      acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               >  1:152849/1‑53 (MQ=255)
                                                                                       cGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGAt                                        <  2:543799/139‑1 (MQ=255)
                                                                                           tACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAAc                                    >  2:523206/1‑139 (MQ=255)
                                                                                               aCTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               >  2:1208883/1‑44 (MQ=255)
                                                                                               aCTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                               <  1:1208883/44‑1 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                                   <  2:1473480/52‑1 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                                   >  1:1473480/1‑52 (MQ=255)
                                                                                                    cTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              <  2:466249/40‑1 (MQ=255)
                                                                                                    cTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                              >  1:466249/1‑40 (MQ=255)
                                                                                                               aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTgctgct                                                                     <  2:1065164/86‑1 (MQ=255)
                                                                                                               aTGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTgctgct                                                                     >  1:1065164/1‑86 (MQ=255)
                                                                                                                   aaCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                                >  1:1310370/1‑87 (MQ=255)
                                                                                                                   aaCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                                <  2:1310370/87‑1 (MQ=255)
                                                                                                                     cgcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAattatt                                                                                                          >  1:21166/1‑43 (MQ=255)
                                                                                                                     cgcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAattatt                                                                                                          <  2:21166/43‑1 (MQ=255)
                                                                                                                         gTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGa                                         <  2:1515018/104‑1 (MQ=255)
                                                                                                                         gTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGa                                         >  1:1515018/1‑104 (MQ=255)
                                                                                                                         gTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTtgtg      >  1:717120/1‑139 (MQ=255)
                                                                                                                         gTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTtgtg      <  1:826036/139‑1 (MQ=255)
                                                                                                                             ggCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCtt  >  1:1194880/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                    
ACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTTGTGCCTT  >  NC_000913/3939088‑3939351

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