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% 278.2 / NA 86*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/45);  total (41/45)

ACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTT  >  NC_000913/3939088‑3939344
                                                                                                                                   |                                                                                                                             
aCTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        <  2:822677/139‑1 (MQ=255)
  tATCCACCAACCGAAAGATGAACTGGGGTAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCt                                                                                                                      <  1:918813/139‑1 (MQ=255)
    tCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGt                                                                                                                    >  1:427090/1‑139 (MQ=255)
    tCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGt                                                                                                                    >  1:1471864/1‑139 (MQ=255)
       accaACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCtt                                                                                                                 >  1:398049/1‑139 (MQ=255)
            ccGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                            <  2:427090/139‑1 (MQ=255)
                      aaCTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAAttatta                                                                                                  >  2:1185053/1‑139 (MQ=255)
                        cTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTAcc                                                                                                <  1:1185053/139‑1 (MQ=255)
                             ggAGCTGGCGATTGATGTACTCGTCCATCGGATACCCCAGCCGACCCTTCAGCAACATCGATTACTACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCgg                                                                                           <  2:1415390/139‑1 (MQ=255)
                                  tGGCGATTGATGTACTCATCCATCGGATAACCCAGTCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGgttgtt                                                                                      >  1:632823/1‑139 (MQ=255)
                                        ttGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                <  2:398049/139‑1 (MQ=255)
                                            tGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCa                                                                            >  1:782238/1‑139 (MQ=255)
                                              tACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                       >  1:860140/1‑94 (MQ=255)
                                              tACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                       <  2:860140/94‑1 (MQ=255)
                                              tACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATg                                                                          >  1:568636/1‑139 (MQ=255)
                                              tACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATg                                                                          >  2:1042406/1‑139 (MQ=255)
                                                 tCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAAc                                                                       <  2:1471864/139‑1 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCttt                                                                                                                >  2:771016/1‑94 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCttt                                                                                                                <  1:771016/94‑1 (MQ=255)
                                                     ccatcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                   <  2:245860/139‑1 (MQ=255)
                                                        tcGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctg                                                                >  2:1322238/1‑139 (MQ=255)
                                                          ggATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        <  2:571601/81‑1 (MQ=255)
                                                          ggATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        >  1:571601/1‑81 (MQ=255)
                                                           gATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        <  2:828503/80‑1 (MQ=255)
                                                           gATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        >  1:828503/1‑80 (MQ=255)
                                                               aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                <  2:835233/116‑1 (MQ=255)
                                                               aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                >  1:835233/1‑116 (MQ=255)
                                                               aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                         >  1:1464259/1‑139 (MQ=255)
                                                               aCCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                         <  1:1348697/139‑1 (MQ=255)
                                                                  cAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                      <  2:1464259/139‑1 (MQ=255)
                                                                    gCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    <  2:568636/139‑1 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                            <  1:386773/82‑1 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                            >  2:386773/1‑82 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                               >  2:592912/1‑127 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtgctgc                                                               <  1:592912/127‑1 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                         >  2:207327/1‑133 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                         <  1:207327/133‑1 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGt                                                   >  2:348048/1‑139 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGt                                                   >  1:1387916/1‑139 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGt                                                   >  2:754523/1‑139 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACAATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                       <  2:1482168/71‑1 (MQ=255)
                                                                     ccGACCCTTCAGCAACAACAATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                       >  1:1482168/1‑71 (MQ=255)
                                                                        aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                <  2:632823/139‑1 (MQ=255)
                                                                        aCCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                <  1:452316/139‑1 (MQ=255)
                                                                           cTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAg                                                                                 >  2:1263541/1‑103 (MQ=255)
                                                                           cTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAg                                                                                 <  1:1263541/103‑1 (MQ=255)
                                                                              cAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        >  2:534691/1‑61 (MQ=255)
                                                                              cAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        <  1:534691/61‑1 (MQ=255)
                                                                              cAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGcc                                          <  2:782238/139‑1 (MQ=255)
                                                                                gCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCa                                        <  1:1436146/139‑1 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        >  2:1179534/1‑55 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACg                                                                                                                        <  1:1179534/55‑1 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                  >  1:960773/1‑93 (MQ=255)
                                                                                    caacGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCa                                                                                  <  2:960773/93‑1 (MQ=255)
                                                                                    caacGATTACAACTCACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCtt                                                                                                                 >  2:780922/1‑62 (MQ=255)
                                                                                    caacGATTACAACTCACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCtt                                                                                                                 <  1:780922/62‑1 (MQ=255)
                                                                                      acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGa                                  >  1:511158/1‑139 (MQ=255)
                                                                                      acGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGa                                  <  1:1322238/139‑1 (MQ=255)
                                                                                         aTTACAACATACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                       <  2:997868/51‑1 (MQ=255)
                                                                                         aTTACAACATACTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGc                                                                                                                       >  1:997868/1‑51 (MQ=255)
                                                                                            aCAACTTACTCCGATTGTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAAGAATGCCCATTGTGATAGACg                            <  1:1042406/139‑1 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                            <  1:667270/52‑1 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGAt                                                                                                            >  2:667270/1‑52 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                      <  1:1048501/106‑1 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                      >  2:1048501/1‑106 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    >  2:137835/1‑108 (MQ=255)
                                                                                                   aCTCCGATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    <  1:137835/108‑1 (MQ=255)
                                                                                                        gATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                <  2:309558/107‑1 (MQ=255)
                                                                                                        gATTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGtaat                                                >  1:309558/1‑107 (MQ=255)
                                                                                                         aTTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATc                                                                                                           <  2:1088580/47‑1 (MQ=255)
                                                                                                         aTTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATc                                                                                                           >  1:1088580/1‑47 (MQ=255)
                                                                                                         aTTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                <  2:428435/74‑1 (MQ=255)
                                                                                                         aTTCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGa                                                                                >  1:428435/1‑74 (MQ=255)
                                                                                                          ttCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                         >  2:91363/1‑96 (MQ=255)
                                                                                                          ttCTGATGGAACGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGc                                                         <  1:91363/96‑1 (MQ=255)
                                                                                                                   aaCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                      <  2:1205813/90‑1 (MQ=255)
                                                                                                                   aaCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAAc                                                      >  1:1205813/1‑90 (MQ=255)
                                                                                                                    aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                   >  1:304865/1‑76 (MQ=255)
                                                                                                                    aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATtg                                                                   <  2:304865/76‑1 (MQ=255)
                                                                                                                    aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    <  1:286051/91‑1 (MQ=255)
                                                                                                                    aCGCGGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    >  2:286051/1‑91 (MQ=255)
                                                                                                                      gcgGTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGtt  <  2:511158/139‑1 (MQ=255)
                                                                                                                       cgGTTCGGCTTAAATTTACGCTGTCTTTTGATC‑AATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGa                                                            <  1:1505530/79‑1 (MQ=255)
                                                                                                                       cgGTTCGGCTTAAATTTACGCTGTCTTTTGATC‑AATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGa                                                            >  2:1505530/1‑79 (MQ=255)
                                                                                                                         gTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    >  1:991655/1‑86 (MQ=255)
                                                                                                                         gTTCGGCTTAAATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAg                                                    <  2:991655/86‑1 (MQ=255)
                                                                                                                                   |                                                                                                                             
ACTATCCACCAACCGAAAGATGAACTGGGGGAGCTGGCGATTGATGTACTCATCCATCGGATAACCCAGCCGACCCTTCAGCAACAACGATTACAACTTACTCCGATTCTGATGGAACGCGGTTCGGCTTAGATTTACGCTGTCTTTTGATCAAATTATTACCATCGGTTGTTTTCAGAAGCATGAACATTGCTGCTGAAGCAACAGTAATAATGCCCATTGTGATAAACGTATAGTGGAATTGTTCGACAGTCGTT  >  NC_000913/3939088‑3939344

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