Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 372,966 A→G T16A (ACC→GCC)  mhpF → acetaldehyde‑CoA dehydrogenase II, NAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913372,9660AG100.0% 284.3 / NA 78T16A (ACC→GCC) mhpFacetaldehyde‑CoA dehydrogenase II, NAD‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/39);  total (39/39)

TTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372848‑373099
                                                                                                                      |                                                                                                                                     
ttGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAAtt                                                                                                                   <  2:639328/139‑1 (MQ=255)
  gAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAAtttt                                                                                                                 <  2:1069823/139‑1 (MQ=255)
        cATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCa                                                                                                           >  1:656976/1‑139 (MQ=255)
         atatTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               >  1:1628629/1‑118 (MQ=255)
         atatTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               <  2:1628629/118‑1 (MQ=255)
         atatTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCAc                                                                                                          >  2:180294/1‑139 (MQ=255)
           atTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACgg                                                                                                        <  2:1247651/139‑1 (MQ=255)
             tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:487219/139‑1 (MQ=255)
             tGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  2:656976/139‑1 (MQ=255)
               aaGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAg                                                                                                    <  1:1328861/139‑1 (MQ=255)
                aGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  2:1615566/1‑136 (MQ=255)
                aGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:1615566/136‑1 (MQ=255)
                aGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGc                                                                                                   >  2:1523486/1‑139 (MQ=255)
                 ggCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCa                                                                                                  >  1:80454/1‑139 (MQ=255)
                   cATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               >  1:52756/1‑108 (MQ=255)
                   cATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               <  2:52756/108‑1 (MQ=255)
                    aTAGGTTCAGTTGCTGCGACATTTTCAAGTGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCt                                                                                               >  2:1644210/1‑139 (MQ=255)
                      aGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTgg                                                                                             >  2:848414/1‑139 (MQ=255)
                           cAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATg                                                                                        <  1:848414/139‑1 (MQ=255)
                                   gcGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGAt                                                                                >  2:312027/1‑139 (MQ=255)
                                         ttttCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTgg                                                                          <  1:1372024/139‑1 (MQ=255)
                                                   cAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               >  1:1152335/1‑76 (MQ=255)
                                                   cAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCtg                                                                                                                               <  2:1152335/76‑1 (MQ=255)
                                                   cAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCt                                                                >  1:1779582/1‑139 (MQ=255)
                                                   cAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCt                                                                <  2:190264/139‑1 (MQ=255)
                                                    aGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACgg                                                                                                        >  2:502111/1‑98 (MQ=255)
                                                    aGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACgg                                                                                                        <  1:502111/98‑1 (MQ=255)
                                                                aaGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTc                                                   >  2:1250195/1‑139 (MQ=255)
                                                                  gTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:23168/86‑1 (MQ=255)
                                                                  gTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  2:23168/1‑86 (MQ=255)
                                                                  gTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTg                                                 >  2:563878/1‑139 (MQ=255)
                                                                    cTGTCATGGGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgc                                               >  1:1714954/1‑139 (MQ=255)
                                                                    cTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgc                                               <  1:1250195/139‑1 (MQ=255)
                                                                     tGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:256449/83‑1 (MQ=255)
                                                                     tGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  2:256449/1‑83 (MQ=255)
                                                                     tGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCa                                                                        >  2:1770989/1‑113 (MQ=255)
                                                                     tGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCa                                                                        <  1:1770989/113‑1 (MQ=255)
                                                                       tCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcgcg                                            <  2:80454/139‑1 (MQ=255)
                                                                        cATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  <  2:1562184/100‑1 (MQ=255)
                                                                        cATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  >  1:1562184/1‑100 (MQ=255)
                                                                          tGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcgcgcgc                                         <  1:1523486/139‑1 (MQ=255)
                                                                            aGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCa                                       <  2:1779582/139‑1 (MQ=255)
                                                                              tAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATg                                                                                        >  1:839834/1‑88 (MQ=255)
                                                                              tAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATg                                                                                        <  2:839834/88‑1 (MQ=255)
                                                                                 gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  2:699770/71‑1 (MQ=255)
                                                                                 gCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  1:699770/1‑71 (MQ=255)
                                                                                   gTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTAt                                >  2:669527/1‑139 (MQ=255)
                                                                                      aaGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACgg                                                                                                        <  1:1581267/64‑1 (MQ=255)
                                                                                      aaGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACgg                                                                                                        >  2:1581267/1‑64 (MQ=255)
                                                                                      aaGTCGCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATggg                             >  2:1690192/1‑139 (MQ=255)
                                                                                           gCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGgag                                                                                           >  2:1603711/1‑72 (MQ=255)
                                                                                           gCCATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGgag                                                                                           <  1:1603711/72‑1 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  2:246235/60‑1 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  1:246235/1‑60 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      <  1:1599897/60‑1 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTc                                                                                                      >  2:1599897/1‑60 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCa                                                                                                     <  1:1650639/61‑1 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCa                                                                                                     >  2:1650639/1‑61 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCa                                                              <  1:1325858/100‑1 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCa                                                              >  2:1325858/1‑100 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCg                                                         <  2:81606/105‑1 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCg                                                         >  1:81606/1‑105 (MQ=255)
                                                                                            ccATTATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGc                       >  1:449973/1‑139 (MQ=255)
                                                                                               ttATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    <  1:669527/139‑1 (MQ=255)
                                                                                               ttATCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGGCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGccac                    <  1:1644210/139‑1 (MQ=255)
                                                                                                  tCGGTTCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCcaccac                 <  1:563878/139‑1 (MQ=255)
                                                                                                       tCTGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATg            >  1:608052/1‑139 (MQ=255)
                                                                                                         tGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGc                       >  2:1405733/1‑126 (MQ=255)
                                                                                                         tGGCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGc                       <  1:1405733/126‑1 (MQ=255)
                                                                                                          ggCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  >  1:1075771/1‑66 (MQ=255)
                                                                                                          ggCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTg                                                                                  <  2:1075771/66‑1 (MQ=255)
                                                                                                           gCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGt                                                                             <  2:1031469/70‑1 (MQ=255)
                                                                                                           gCAACATTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGt                                                                             >  1:1031469/1‑70 (MQ=255)
                                                                                                                aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcg                                              >  2:1253379/1‑96 (MQ=255)
                                                                                                                aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGgcg                                              <  1:1253379/96‑1 (MQ=255)
                                                                                                                aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGTGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGa   >  2:1003666/1‑139 (MQ=255)
                                                                                                                aTTGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGa   >  1:1627204/1‑139 (MQ=255)
                                                                                                                 ttGGTGCCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAt  <  1:312027/139‑1 (MQ=255)
                                                                                                                      |                                                                                                                                     
TTGAAGCCCATATTGAAGGCATAGGTTCAGTTGCTGCGACATTTTCAAGCGCAGCCCCAAAAGGAAGTCTGTCATGAGTAAGCGTAAAGTCGCCATTATCGGTTCTGGCAACATTGGTACCGATCTGATGATTAAAATTTTGCGTCACGGTCAGCATCTGGAGATGGCGGTGATGGTTGGCATTGATCCTCAGTCCGACGGTCTGGCGCGCGCCAGACGTATGGGCGTCGCCACCACCCATGAAGGGGTGAT  >  NC_000913/372848‑373099

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