Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 631,934 A→G T681A (ACG→GCG)  cstA → carbon starvation protein involved in peptide utilization, APC peptide transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913631,9340AG100.0% 243.0 / NA 65T681A (ACG→GCG) cstAcarbon starvation protein involved in peptide utilization, APC peptide transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/31);  total (34/31)

GTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACC  >  NC_000913/631808‑632066
                                                                                                                              |                                                                                                                                    
gtgtTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAAc                                                                                                                          <  2:1412576/139‑1 (MQ=255)
gtgtTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAAc                                                                                                                          <  1:166878/139‑1 (MQ=255)
  gtTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCa                                                                                                                        >  2:1391268/1‑139 (MQ=255)
  gtTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCa                                                                                                                        <  2:306401/139‑1 (MQ=255)
          acCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTTGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCTGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGa                                                                                                                >  1:1122744/1‑139 (MQ=255)
             gTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGaaaa                                                                                                             >  2:1499148/1‑139 (MQ=255)
                tGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGt                                                                                                          >  1:44252/1‑139 (MQ=255)
                  gATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCg                                                                                                        >  1:1326047/1‑139 (MQ=255)
                      ccGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCgagga                                                                                                    >  1:1110313/1‑139 (MQ=255)
                            tAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGt                                                                                              >  2:1280543/1‑139 (MQ=255)
                                cATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGgcg                                                                                          <  1:1470337/139‑1 (MQ=255)
                                        ttATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCaaa                                                                                  >  2:1143078/1‑139 (MQ=255)
                                            ggtggtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGa                                                                                                                >  2:714499/1‑105 (MQ=255)
                                            ggtggtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGa                                                                                                                <  1:714499/105‑1 (MQ=255)
                                            ggtggtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgc                                                                              <  1:1475244/139‑1 (MQ=255)
                                              tggtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcgc                                                                            <  2:44252/139‑1 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGcc                                                                                                                   >  1:1489712/1‑98 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGcc                                                                                                                   <  2:1489712/98‑1 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          >  2:112460/1‑139 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          >  1:119510/1‑139 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          >  2:33086/1‑139 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          >  2:386467/1‑139 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          >  1:216591/1‑139 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          <  2:1318728/139‑1 (MQ=255)
                                                gtCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGcac                                                                          >  1:1318728/1‑139 (MQ=255)
                                                 tCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcgc                                                                            >  2:1226181/1‑136 (MQ=255)
                                                 tCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcgc                                                                            <  1:1226181/136‑1 (MQ=255)
                                                   gtggtgGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACt                                                                       <  1:1280543/139‑1 (MQ=255)
                                                        ggtTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGt                                                                  <  2:359772/139‑1 (MQ=255)
                                                               gCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCg                                                                                                                                    >  1:1344928/1‑66 (MQ=255)
                                                                    gTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAgggg                                                      >  1:1195177/1‑139 (MQ=255)
                                                                     ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCag                                                               <  2:1472673/129‑1 (MQ=255)
                                                                     ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCag                                                               >  1:1472673/1‑129 (MQ=255)
                                                                     ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggg                                                     <  2:385495/139‑1 (MQ=255)
                                                                     ttCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAggggg                                                     >  2:768973/1‑139 (MQ=255)
                                                                               aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATg                                                                                                           >  2:36264/1‑75 (MQ=255)
                                                                               aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATg                                                                                                           <  1:36264/75‑1 (MQ=255)
                                                                               aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcgc                                                                            <  1:137992/106‑1 (MQ=255)
                                                                               aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcgc                                                                            >  2:137992/1‑106 (MQ=255)
                                                                               aGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGc                                           >  1:705322/1‑139 (MQ=255)
                                                                                     cACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcg                                                                             >  2:338153/1‑99 (MQ=255)
                                                                                     cACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcg                                                                             <  1:338153/99‑1 (MQ=255)
                                                                                       cTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTg                                   <  1:768973/139‑1 (MQ=255)
                                                                                        ttGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcg                                                                             <  2:408634/96‑1 (MQ=255)
                                                                                        ttGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGgcg                                                                             >  1:408634/1‑96 (MQ=255)
                                                                                          gCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCt                                <  1:386467/139‑1 (MQ=255)
                                                                                            ggCATTGAAAGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTa                              <  1:1315607/139‑1 (MQ=255)
                                                                                                   aaaGATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCa                       <  1:1499148/139‑1 (MQ=255)
                                                                                                      gATCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAggg                    <  1:33086/139‑1 (MQ=255)
                                                                                                       aTCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                   <  2:119510/139‑1 (MQ=255)
                                                                                                        tCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAAt                                       <  1:1076117/118‑1 (MQ=255)
                                                                                                        tCCGAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAAt                                       >  2:1076117/1‑118 (MQ=255)
                                                                                                           gAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTa                              <  2:138667/124‑1 (MQ=255)
                                                                                                           gAAGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTa                              >  1:138667/1‑124 (MQ=255)
                                                                                                             aGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc             >  1:897471/1‑139 (MQ=255)
                                                                                                             aGCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctc             >  1:1139495/1‑139 (MQ=255)
                                                                                                              gCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTgg                                            <  2:1372893/107‑1 (MQ=255)
                                                                                                              gCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTgg                                            >  1:1372893/1‑107 (MQ=255)
                                                                                                              gCCAACGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCctct            <  2:1122744/139‑1 (MQ=255)
                                                                                                                 aaCGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAAt                                                                                                            >  1:158585/1‑40 (MQ=255)
                                                                                                                 aaCGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAAt                                                                                                            <  2:158585/40‑1 (MQ=255)
                                                                                                                   cGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCa                       >  1:560495/1‑123 (MQ=255)
                                                                                                                   cGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCa                       <  2:560495/123‑1 (MQ=255)
                                                                                                                   cGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTAGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                   >  1:682848/1‑127 (MQ=255)
                                                                                                                   cGGCGAAAGAAGCGCCGTATGAACCAATGCCGGAAAATGTAGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGt                   <  2:682848/127‑1 (MQ=255)
                                                                                                                       gaaagaagCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAc   <  2:1139495/139‑1 (MQ=255)
                                                                                                                        aaagaagCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAAcc  <  1:428975/139‑1 (MQ=255)
                                                                                                                         aagaagCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTg                                                           <  1:748895/81‑1 (MQ=255)
                                                                                                                         aagaagCGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTg                                                           >  2:748895/1‑81 (MQ=255)
                                                                                                                              |                                                                                                                                    
GTGTTCAACAACCGTCTGGATGCCGGGTTAACCATCTTCTTTATGGTGGTCGTGGTGGTTCTGGCACTGTTCTCGATTAAGACGGCACTTGCGGCATTGAAAGATCCGAAGCCAACGGCGAAAGAAACGCCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACC  >  NC_000913/631808‑632066

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: