Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,257,838 +G coding (2222/2241 nt) tdcE ← putative formate acetyltransferase 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,257,8371.G99.7% 1030.1 / ‑4.6 328A741A (GCG→GCCtdcEputative formate acetyltransferase 3
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base G (165/161);  minor base T (0/1);  total (166/162)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 350 total aligned reads displayed.

TTTTTCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTG  >  NC_000913/3257745‑3257924
                                                                                             |                                                                                       
 ttttCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAcg                                                                                        <  1:5558665/93‑1 (MQ=255)
 ttttCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAcg                                                                                        <  1:2332544/93‑1 (MQ=255)
    tCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCAcggcg                                                                                     >  1:6565303/1‑90 (MQ=255)
       tACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAg                                                                                   >  1:3868225/1‑93 (MQ=255)
       tACACCCTCCGTCGTGAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAg                                                                                   <  2:4471604/93‑1 (MQ=255)
        aCATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGc                                                                                  <  1:4582494/93‑1 (MQ=255)
          aTCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCa                                                                                >  1:1655761/1‑93 (MQ=255)
            cctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCaga                                                                              >  2:2881599/1‑93 (MQ=255)
            cctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCaga                                                                              <  1:5266598/93‑1 (MQ=255)
            cctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCaga                                                                              >  2:4890918/1‑93 (MQ=255)
             ctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCagag                                                                             >  2:457442/1‑93 (MQ=255)
             ctccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCagag                                                                             >  1:6236374/1‑93 (MQ=255)
              tccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAgaga                                                                            <  1:5427188/93‑1 (MQ=255)
              tccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAgaga                                                                            <  1:5892802/93‑1 (MQ=255)
              tccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAgaga                                                                            >  1:7500665/1‑93 (MQ=255)
              cccGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAgaga                                                                            <  2:3553367/92‑1 (MQ=255)
                cGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGaca                                                                          <  1:465435/93‑1 (MQ=255)
                cGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGaca                                                                          <  2:5092909/93‑1 (MQ=255)
                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                        >  1:2909689/1‑93 (MQ=255)
                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                        >  1:561192/1‑93 (MQ=255)
                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                        >  1:2381964/1‑93 (MQ=255)
                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                        >  1:244952/1‑93 (MQ=255)
                  gCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg                                                                        >  1:3146537/1‑93 (MQ=255)
                    gTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGa                                                                      <  2:5635828/93‑1 (MQ=255)
                    gTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGa                                                                      <  2:1325155/93‑1 (MQ=255)
                     tCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGAt                                                                     >  2:5294990/1‑93 (MQ=255)
                     tCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGAt                                                                     >  1:1444425/1‑93 (MQ=255)
                         agCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTc                                                                 <  1:2290424/93‑1 (MQ=255)
                         agCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTc                                                                 <  2:2306492/93‑1 (MQ=255)
                          gcgcCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                                <  2:1933275/93‑1 (MQ=255)
                          gcgcCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                                <  2:3804782/93‑1 (MQ=255)
                          gcgcCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                                <  2:3245485/93‑1 (MQ=255)
                          gcgcCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCa                                                                <  1:368431/93‑1 (MQ=255)
                             ccTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGt                                                             >  2:1662159/1‑93 (MQ=255)
                              cTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGtt                                                            >  2:1889573/1‑93 (MQ=255)
                               tGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTa                                                           >  1:3284401/1‑93 (MQ=255)
                               tGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTa                                                           <  2:4711420/93‑1 (MQ=255)
                                gggTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAg                                                          <  2:7208378/93‑1 (MQ=255)
                                 ggTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAgg                                                         <  1:5103411/93‑1 (MQ=255)
                                 ggTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAgg                                                         <  1:6388210/93‑1 (MQ=255)
                                  gTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGa                                                        <  2:6794303/93‑1 (MQ=255)
                                   tAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGat                                                       >  1:5079673/1‑93 (MQ=255)
                                   tAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGat                                                       <  2:2203599/93‑1 (MQ=255)
                                   tAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGat                                                       >  2:3965664/1‑93 (MQ=255)
                                   tAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGat                                                       >  2:5869844/1‑93 (MQ=255)
                                    taaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCTTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGata                                                      <  1:7853224/92‑1 (MQ=255)
                                     aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGatat                                                     >  2:7304876/1‑93 (MQ=255)
                                     aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGatat                                                     >  2:7436072/1‑93 (MQ=255)
                                      aGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAtt                                                    >  1:7690421/1‑93 (MQ=255)
                                      aGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAtt                                                    <  2:6034479/93‑1 (MQ=255)
                                       ggTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttt                                                   >  1:4530464/1‑93 (MQ=255)
                                       ggTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttt                                                   >  2:7742532/1‑93 (MQ=255)
                                       agTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttt                                                   <  1:4149379/92‑1 (MQ=255)
                                         tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttttt                                                 >  2:6661699/1‑93 (MQ=255)
                                         tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATAttttt                                                 >  1:1071065/1‑93 (MQ=255)
                                          aCGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTc                                                <  1:2708616/93‑1 (MQ=255)
                                           cGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                                               >  2:7544297/1‑93 (MQ=255)
                                           cGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTcc                                               >  2:5332075/1‑93 (MQ=255)
                                            gTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCg                                              <  1:7165471/93‑1 (MQ=255)
                                            gTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCg                                              <  2:2082769/93‑1 (MQ=255)
                                              gAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCggg                                            >  1:4724935/1‑93 (MQ=255)
                                              gAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCggg                                            <  2:4026025/93‑1 (MQ=255)
                                              gAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCggg                                            <  2:902591/93‑1 (MQ=255)
                                               aaaTAGCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCGCGGCGTAGCCAGAGACAAGGATTGTCGGGTTAGGATATTTTTCCGGGt                                           <  2:2751411/93‑1 (MQ=37)
                                               aaaTAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGt                                           <  1:1207174/93‑1 (MQ=255)
                                                aaTAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTg                                          <  1:993976/93‑1 (MQ=255)
                                                 aTAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGc                                         >  2:420872/1‑93 (MQ=255)
                                                  tAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCt                                        >  2:7919030/1‑93 (MQ=255)
                                                  tAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCt                                        >  1:7295650/1‑93 (MQ=255)
                                                  tAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCt                                        <  1:7808610/93‑1 (MQ=255)
                                                    aCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCg                                      >  1:5800578/1‑93 (MQ=255)
                                                    aCATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCg                                      >  2:6506062/1‑93 (MQ=255)
                                                     cATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGa                                     <  1:1923697/93‑1 (MQ=255)
                                                      aTCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGAt                                    <  1:3468428/93‑1 (MQ=255)
                                                          tGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCa                                >  2:5754183/1‑93 (MQ=255)
                                                          tGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCa                                <  2:1191743/93‑1 (MQ=255)
                                                           gTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCAt                               <  2:1496154/93‑1 (MQ=255)
                                                           gTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCAt                               >  1:7055414/1‑93 (MQ=255)
                                                              gCTGTTCACGGGTCAGTGCGTTGAAGCGCACTGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCca                            <  1:4909459/93‑1 (MQ=255)
                                                               cTGTTCACGTGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcag                           <  2:3174607/93‑1 (MQ=255)
                                                                tGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagc                          >  2:6187051/1‑93 (MQ=255)
                                                                tGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagc                          >  1:1508986/1‑93 (MQ=255)
                                                                tGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagc                          >  1:6756412/1‑93 (MQ=255)
                                                                 gTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCcagca                         <  2:7884176/93‑1 (MQ=255)
                                                                    cACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAgcagca                      <  1:907377/93‑1 (MQ=255)
                                                                       gggTCAGTGCGTTGACGCGCACGGCGTAGCCAGAGACACGGATAGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAatt                   <  2:7027837/93‑1 (MQ=37)
                                                                       gggTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAttt                   >  2:1994592/1‑93 (MQ=255)
                                                                        ggTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTc                  <  1:6880760/93‑1 (MQ=255)
                                                                         gTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTcc                 <  1:4845171/93‑1 (MQ=255)
                                                                          tCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTccc                >  1:1997065/1‑93 (MQ=255)
                                                                          tCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTccc                >  2:514633/1‑93 (MQ=255)
                                                                          gcAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACGCGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTccc                <  1:1795772/92‑1 (MQ=37)
                                                                                cGTTGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCa          >  1:4046733/1‑93 (MQ=255)
                                                                                  ttGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCAtt        >  2:2461441/1‑93 (MQ=255)
                                                                                  ttGAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCAtt        <  2:2614848/93‑1 (MQ=255)
                                                                                    gAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTac      <  2:6688143/93‑1 (MQ=255)
                                                                                    gAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTac      <  2:4386783/93‑1 (MQ=255)
                                                                                    gAAGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTac      <  2:6281144/93‑1 (MQ=255)
                                                                                     aaGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTacg     <  1:3578450/93‑1 (MQ=255)
                                                                                      aGCGCACGGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTacgt    >  2:6557910/1‑93 (MQ=255)
                                                                                             |                                                                                       
TTTTTCATACATCCTCCGGCGTCAGAGCGCCTGGGTAAAGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTG  >  NC_000913/3257745‑3257924

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: