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% 803.5 / NA 244Q44R (CAG→CGG) hemGprotoporphyrin oxidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (127/117);  total (127/117)

Only 100 of 252 total aligned reads displayed.

CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416
                                                                                          |                                                                                         
cgcgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGt                                                                                         >  2:6552456/1‑93 (MQ=255)
 gcgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTg                                                                                        >  2:6253850/1‑93 (MQ=255)
  cgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTgg                                                                                       <  1:1691968/93‑1 (MQ=255)
  cgcgAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTgg                                                                                       <  1:2813205/93‑1 (MQ=255)
     gagaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaa                                                                                    >  2:4660647/1‑93 (MQ=255)
     gagaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaa                                                                                    >  2:5481425/1‑93 (MQ=255)
      agaTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGaaa                                                                                   <  1:3533919/93‑1 (MQ=255)
        aTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAAc                                                                                 <  1:4111307/93‑1 (MQ=255)
        aTTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAAc                                                                                 <  2:5049598/93‑1 (MQ=255)
         ttGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACt                                                                                >  2:2650862/1‑93 (MQ=255)
         ttGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACt                                                                                <  1:2489192/93‑1 (MQ=255)
          tGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTa                                                                               <  1:1239283/93‑1 (MQ=255)
             ctcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGa                                                                            <  2:680106/93‑1 (MQ=255)
              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  2:5616481/1‑93 (MQ=255)
              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  1:514321/1‑93 (MQ=255)
              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           <  2:2728693/93‑1 (MQ=255)
              tcctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGAc                                                                           >  1:1998023/1‑93 (MQ=255)
                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         >  1:2569608/1‑93 (MQ=255)
                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCg                                                                         <  2:3385667/93‑1 (MQ=255)
                ctACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTAGGAAAACTATGACCg                                                                         <  1:637815/93‑1 (MQ=255)
                  aCNTGGNTTCGGAACTGAAAGANCTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtg                                                                       <  2:5877229/93‑1 (MQ=255)
                  aCCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCgtg                                                                       <  1:2979690/93‑1 (MQ=255)
                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGTCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     >  2:2152793/1‑93 (MQ=255)
                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     >  2:130260/1‑93 (MQ=255)
                    cTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGtgtg                                                                     <  1:39389/93‑1 (MQ=255)
                     tGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTgg                                                                    >  1:6075938/1‑93 (MQ=255)
                       gCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTc                                                                  >  2:1670863/1‑93 (MQ=255)
                        cTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCa                                                                 >  2:4194255/1‑93 (MQ=255)
                        cTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCa                                                                 <  2:5354703/93‑1 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  2:1983977/1‑93 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                <  1:4208392/93‑1 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  2:1415839/1‑93 (MQ=255)
                         ttCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAt                                                                >  1:1826551/1‑93 (MQ=255)
                          tCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAtt                                                               <  1:1476514/93‑1 (MQ=255)
                          tCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAtt                                                               <  2:2877076/93‑1 (MQ=255)
                          tCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAtt                                                               <  1:6293414/93‑1 (MQ=255)
                          tCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCAtt                                                               <  2:5321077/93‑1 (MQ=255)
                           cGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTg                                                              <  1:1273580/93‑1 (MQ=255)
                           cGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTg                                                              <  1:240003/93‑1 (MQ=255)
                            ggAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTgg                                                             >  1:2959923/1‑93 (MQ=255)
                            ggAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTgg                                                             >  2:5146651/1‑93 (MQ=255)
                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGc                                                          >  1:3351790/1‑93 (MQ=255)
                               aCTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCACTGGTGc                                                          >  2:1502336/1‑93 (MQ=255)
                                 tGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCtt                                                        >  2:4461438/1‑93 (MQ=255)
                                  gAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTc                                                       >  2:458771/1‑93 (MQ=255)
                                   aaaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCt                                                      <  2:2618429/93‑1 (MQ=255)
                                   aaaGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCt                                                      <  2:4855764/93‑1 (MQ=255)
                                     aGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAt                                                    >  2:5509475/1‑93 (MQ=255)
                                      gAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTAtt                                                   <  1:3344481/93‑1 (MQ=255)
                                         cTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGCACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGc                                                <  2:3811394/93‑1 (MQ=255)
                                         cTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGc                                                <  1:4560355/93‑1 (MQ=255)
                                         cTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGc                                                <  1:3840111/93‑1 (MQ=255)
                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              >  2:4981198/1‑93 (MQ=255)
                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  1:4491341/93‑1 (MQ=255)
                                           ggggATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTa                                              <  2:1855956/93‑1 (MQ=255)
                                              gATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATgg                                           >  2:2799566/1‑93 (MQ=255)
                                               aTCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGt                                          >  1:3651920/1‑93 (MQ=255)
                                                tCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTc                                         >  2:2202863/1‑93 (MQ=255)
                                                  cAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCAc                                       >  1:2316014/1‑93 (MQ=255)
                                                    ggCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTa                                     >  1:5645949/1‑93 (MQ=255)
                                                     gCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTAc                                    <  2:6155514/93‑1 (MQ=255)
                                                       ggATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCa                                  >  1:5201642/1‑93 (MQ=255)
                                                        gATGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCAt                                 <  1:3400842/93‑1 (MQ=255)
                                                          tGTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTc                               >  2:2850637/1‑93 (MQ=255)
                                                           gTCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCa                              <  2:3037016/93‑1 (MQ=255)
                                                            tCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAg                             <  1:2283933/93‑1 (MQ=255)
                                                            tCGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAg                             >  1:4364416/1‑93 (MQ=255)
                                                             cGCCAATGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgt                            <  2:3728047/93‑2 (MQ=255)
                                                               cgaaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                          <  1:844957/91‑1 (MQ=255)
                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        <  1:5936686/93‑1 (MQ=255)
                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  1:5881076/1‑93 (MQ=255)
                                                                 aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                        >  2:541561/1‑93 (MQ=255)
                                                                  aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                       >  1:6700303/1‑93 (MQ=255)
                                                                   tgtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCa                      <  2:3928364/93‑1 (MQ=255)
                                                                    gtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAg                     <  2:1383870/93‑1 (MQ=255)
                                                                    gtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAg                     <  1:2151884/93‑1 (MQ=255)
                                                                        aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                 >  2:172159/1‑93 (MQ=255)
                                                                         ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                >  1:4175942/1‑93 (MQ=255)
                                                                         ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                >  2:578285/1‑93 (MQ=255)
                                                                          cGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAttt               <  2:6077539/93‑1 (MQ=255)
                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  2:2333326/1‑93 (MQ=255)
                                                                           gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg              >  1:3444977/1‑93 (MQ=255)
                                                                            cATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGt             >  1:656387/1‑93 (MQ=255)
                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  1:3758131/1‑93 (MQ=255)
                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  1:4385352/1‑93 (MQ=255)
                                                                              ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa           >  2:1289003/1‑93 (MQ=255)
                                                                                gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa         >  2:113424/1‑93 (MQ=255)
                                                                                gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa         >  2:412432/1‑93 (MQ=255)
                                                                                 aagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaa        >  1:5690364/1‑93 (MQ=255)
                                                                                  agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaa       >  2:5303163/1‑93 (MQ=255)
                                                                                  agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaa       >  2:2358501/1‑93 (MQ=255)
                                                                                  agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTAAGTGTTCCAGGAATTTGTCaaaaa       >  1:1948349/1‑93 (MQ=255)
                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  1:2860029/93‑1 (MQ=255)
                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    <  1:4650334/93‑1 (MQ=255)
                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    >  2:4520200/1‑93 (MQ=255)
                                                                                     aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATACGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa    >  1:6032937/1‑93 (MQ=255)
                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  >  2:5677318/1‑93 (MQ=255)
                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  <  1:4608178/93‑1 (MQ=255)
                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  <  1:4097113/93‑1 (MQ=255)
                                                                                       cacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATg  <  1:1639197/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCGCGAGATTGCCTCCTACCTGGCTTCGGAACTGAAAGAACTGGGGATCCAGGCGGATGTCGCCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATG  >  NC_000913/4032237‑4032416

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