Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,032,327 A→G Q44R (CAG→CGG)  hemG → protoporphyrin oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,032,3270AG100.0% 748.3 / NA 231Q44R (CAG→CGG) hemGprotoporphyrin oxidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (114/117);  total (114/117)

Only 100 of 237 total aligned reads displayed.

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416
                                                                                          |                                                                                         
cgcgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGt                                                                                         <  1:3636534/93‑1 (MQ=255)
  cgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTgg                                                                                       >  1:5809976/1‑93 (MQ=255)
      agaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaaa                                                                                   >  1:1978905/1‑93 (MQ=255)
        aTTGCCTCCTACCTGGCTTCGGAAATGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAAc                                                                                 >  2:937716/1‑93 (MQ=255)
         ttGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGCAAACt                                                                                >  2:4819064/1‑93 (MQ=255)
           gCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTAt                                                                              >  2:2185330/1‑93 (MQ=255)
           gCCTCCTACCAGGCTTCGGAACTGAACGTACAGGGGAACCAGACAGATGTCACCAATGTGCACCGCATTAAAGAACCACGGTGGGACCACTAt                                                                              <  1:6297817/93‑1 (MQ=255)
            ccccctACCTGGCTTCGGAACTGAAAGAACTGGTGATCCAGGCGGATGTCGCCAATGTGCACAGCATTGAAGAACCACGGTGGGAAAACTATg                                                                             <  1:2138638/90‑1 (MQ=255)
              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           <  2:3062977/93‑1 (MQ=255)
              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           <  2:567052/93‑1 (MQ=255)
               cctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAcc                                                                          <  1:1293982/93‑1 (MQ=255)
                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  1:4760294/1‑93 (MQ=255)
                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  2:1418640/1‑93 (MQ=255)
                 tACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgt                                                                        >  1:3371558/1‑93 (MQ=255)
                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGCAAACTATGACCgtgt                                                                      >  1:3012653/1‑93 (MQ=255)
                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      >  2:2676566/1‑93 (MQ=255)
                   ccTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtgt                                                                      >  2:4751157/1‑93 (MQ=255)
                     tGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTgg                                                                    <  2:4586709/93‑1 (MQ=255)
                     tGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTgg                                                                    >  2:3304025/1‑93 (MQ=255)
                     tGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTgg                                                                    >  2:5312413/1‑93 (MQ=255)
                      ggCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGt                                                                   <  1:6153832/93‑1 (MQ=255)
                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  >  2:4980246/1‑93 (MQ=255)
                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  <  1:3488198/93‑1 (MQ=255)
                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  >  2:6153876/1‑93 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                <  1:875011/93‑1 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  1:514967/1‑93 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                <  1:3337860/93‑1 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  1:2754810/1‑93 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  1:2506971/1‑93 (MQ=255)
                          tCGGAACTGAAAGATCTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAGCCACGGTGGGAAACCTATGACCGTGTGGTCAtt                                                               <  1:960003/93‑1 (MQ=255)
                           cGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTg                                                              <  2:2538074/93‑1 (MQ=255)
                           cGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTg                                                              <  1:5903358/93‑1 (MQ=255)
                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          <  1:1075031/93‑1 (MQ=255)
                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          >  1:5476047/1‑93 (MQ=255)
                                cTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCt                                                         <  2:2746890/93‑1 (MQ=255)
                                 tGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCtt                                                        >  2:4483410/1‑93 (MQ=255)
                                  gAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTc                                                       <  1:1694299/93‑1 (MQ=255)
                                  gAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTc                                                       <  1:1555455/93‑1 (MQ=255)
                                      gAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTTGTGCTTCTAtt                                                   <  2:2963104/93‑1 (MQ=255)
                                       aaCTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTc                                                  <  2:2766375/93‑1 (MQ=255)
                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  1:5186975/93‑1 (MQ=255)
                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  2:2238516/93‑1 (MQ=255)
                                            gggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTAt                                             <  2:1321165/93‑1 (MQ=255)
                                            gggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTAt                                             <  1:44734/93‑1 (MQ=255)
                                              gATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATgg                                           >  2:5410593/1‑93 (MQ=255)
                                              gATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATgg                                           <  1:2032265/93‑1 (MQ=255)
                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          <  2:1386875/93‑1 (MQ=255)
                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          >  2:112998/1‑93 (MQ=255)
                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          <  1:4614928/93‑1 (MQ=255)
                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          <  1:3829361/93‑1 (MQ=255)
                                                tCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTc                                         >  1:5543440/1‑93 (MQ=255)
                                                 ccAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCa                                        <  2:2070856/93‑1 (MQ=255)
                                                 ccAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCa                                        <  1:3541063/93‑1 (MQ=255)
                                                  cAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCAc                                       >  1:747858/1‑93 (MQ=255)
                                                   aGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACt                                      >  2:4401456/1‑93 (MQ=255)
                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  1:1002836/1‑93 (MQ=255)
                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  1:482600/1‑93 (MQ=255)
                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  2:1088780/1‑93 (MQ=255)
                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  2:574010/1‑93 (MQ=255)
                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 >  1:4180804/1‑93 (MQ=255)
                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 >  2:1726957/1‑93 (MQ=255)
                                                          tGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTc                               >  1:1597976/1‑93 (MQ=255)
                                                           gTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCa                              >  1:23374/1‑93 (MQ=255)
                                                           gTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCa                              <  2:316446/93‑1 (MQ=255)
                                                             cGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgt                            <  2:1051128/93‑2 (MQ=255)
                                                             cGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgt                            <  1:1349489/93‑2 (MQ=255)
                                                              gCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtg                           >  1:4381495/1‑91 (MQ=255)
                                                              gCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtg                           <  2:2381684/93‑3 (MQ=255)
                                                               ccAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          <  1:4891900/93‑1 (MQ=255)
                                                               ccAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          <  2:6015906/93‑1 (MQ=255)
                                                                cAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGtt                         >  1:4012019/1‑93 (MQ=255)
                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  1:724692/1‑93 (MQ=255)
                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  1:633491/1‑93 (MQ=255)
                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        <  1:82023/93‑1 (MQ=255)
                                                                  aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                       >  2:3779473/1‑93 (MQ=255)
                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      >  1:1364826/1‑93 (MQ=255)
                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      <  1:5992833/93‑1 (MQ=255)
                                                                    gtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAg                     >  1:1142720/1‑93 (MQ=255)
                                                                      gCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGa                   <  1:5768378/93‑1 (MQ=255)
                                                                      gCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGa                   <  1:206371/93‑1 (MQ=255)
                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 <  2:5905640/93‑1 (MQ=255)
                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 >  2:3756224/1‑93 (MQ=255)
                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 >  2:4229558/1‑93 (MQ=255)
                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  1:6195518/1‑93 (MQ=255)
                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  2:3000622/1‑93 (MQ=255)
                                                                            cATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGt             >  1:5085359/1‑93 (MQ=255)
                                                                             aTTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTc            <  1:1829423/93‑1 (MQ=255)
                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  1:2766455/1‑93 (MQ=255)
                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  2:3210624/1‑93 (MQ=255)
                                                                                 aagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGGATTTGTCaaaa        >  2:1556552/1‑93 (MQ=255)
                                                                                 aagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaa        >  1:5282461/1‑93 (MQ=255)
                                                                                  agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaa       <  2:5812587/93‑1 (MQ=255)
                                                                                   gaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaaa      <  1:610691/93‑1 (MQ=255)
                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTTCCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     <  1:608895/93‑1 (MQ=255)
                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     <  1:2716310/93‑1 (MQ=255)
                                                                                    aaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     <  2:4041270/93‑1 (MQ=255)
                                                                                    aaCCACGGTGCGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAAc     <  2:59421/93‑1 (MQ=255)
                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  1:4864396/93‑1 (MQ=255)
                                                                                      ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt   <  2:366172/93‑1 (MQ=255)
                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  <  1:3449765/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

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