Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,263,021 G→A D29N (GAT→AAT)  hemA → HemA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,263,0210GA100.0% 940.4 / NA 251D29N (GAT→AAT) hemAHemA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (131/120);  total (131/120)

Only 100 of 262 total aligned reads displayed.

GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110
                                                                                          |                                                                                         
gCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAGCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAAt                                                                                         >  2:4817135/1‑93 (MQ=255)
 cAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAATCTCAATc                                                                                        <  1:5095370/93‑1 (MQ=255)
 cAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc                                                                                        <  1:2768631/93‑1 (MQ=255)
  aGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  1:499484/93‑1 (MQ=255)
  aGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  2:5152188/93‑1 (MQ=255)
   gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  1:5739957/1‑93 (MQ=255)
   gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  1:3432485/1‑93 (MQ=255)
   gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  1:2709836/1‑93 (MQ=255)
   gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      <  1:145250/93‑1 (MQ=255)
    aCATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAgg                                                                                     <  1:4340473/93‑1 (MQ=255)
     cATGACCGTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    <  1:217009/93‑1 (MQ=255)
     cATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    <  2:2678103/93‑1 (MQ=255)
     cATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgc                                                                                    <  1:1844273/93‑1 (MQ=255)
      aTGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgcg                                                                                   >  1:4662760/1‑93 (MQ=255)
       tGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgcgc                                                                                  <  2:3485794/93‑1 (MQ=255)
         aCCCTTTTAGCACTCGGTATCAACCATAAAGCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtt                                                                                <  2:2866721/93‑1 (MQ=255)
         aCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtt                                                                                >  1:1085629/1‑93 (MQ=255)
          cccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTg                                                                               >  2:2924279/1‑93 (MQ=255)
           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGa                                                                              >  2:3435188/1‑93 (MQ=255)
           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGa                                                                              >  1:746471/1‑93 (MQ=255)
           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGa                                                                              >  2:6575973/1‑93 (MQ=255)
           ccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGa                                                                              <  1:2424565/93‑1 (MQ=255)
             ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGATGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            >  1:5952299/1‑93 (MQ=255)
               ttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGc                                                                          <  2:4142123/93‑1 (MQ=255)
               ttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGc                                                                          <  2:6480234/93‑1 (MQ=255)
                tAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGcc                                                                         <  1:5623857/93‑1 (MQ=255)
                tAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGcc                                                                         <  2:4657258/93‑1 (MQ=255)
                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        <  1:3087373/93‑1 (MQ=255)
                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        <  2:1642288/93‑1 (MQ=255)
                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        >  2:3384855/1‑93 (MQ=255)
                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        <  1:3740334/93‑1 (MQ=255)
                    ccTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     <  2:1276340/92‑1 (MQ=255)
                    aCTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCt                                                                     >  2:3413156/1‑93 (MQ=255)
                      tCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGATCAATCAGGCGCCTGACAGCCTGCTTg                                                                   >  1:2124340/1‑93 (MQ=255)
                       cGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgc                                                                  >  2:6054655/1‑93 (MQ=255)
                         gTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgcgc                                                                >  1:3896974/1‑93 (MQ=255)
                         gTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgcgc                                                                >  2:4098794/1‑93 (MQ=255)
                             cAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGcc                                                            >  1:785330/1‑93 (MQ=255)
                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGAt                                                         <  2:786089/93‑1 (MQ=255)
                                 cATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTTGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATg                                                        >  1:5185038/1‑93 (MQ=255)
                                 cATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATg                                                        >  1:2664884/1‑93 (MQ=255)
                                  aTAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATgg                                                       >  2:373030/1‑93 (MQ=255)
                                   tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      >  1:253251/1‑93 (MQ=255)
                                   tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      >  1:3759451/1‑93 (MQ=255)
                                   tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACACCCTGCTTGCGCAGCCGATggg                                                      >  1:4414224/1‑92 (MQ=255)
                                    aaaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTg                                                     >  2:507458/1‑93 (MQ=255)
                                     aaaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGc                                                    <  2:1157295/93‑1 (MQ=255)
                                      aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   >  1:6497692/1‑93 (MQ=255)
                                      aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   >  2:3680723/1‑93 (MQ=255)
                                       aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  >  2:1771889/1‑93 (MQ=255)
                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 >  2:3575551/1‑93 (MQ=255)
                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 >  2:1397243/1‑93 (MQ=255)
                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 >  2:1975442/1‑93 (MQ=255)
                                        cGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAgg                                                 <  2:2557831/93‑1 (MQ=255)
                                          gCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggc                                               <  2:5825355/93‑1 (MQ=255)
                                          gCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggc                                               <  1:3327700/93‑1 (MQ=255)
                                            aCCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcgg                                             >  1:6080510/1‑93 (MQ=255)
                                            aCCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcgg                                             >  1:1783261/1‑93 (MQ=255)
                                             ccTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcggc                                            >  1:6436929/1‑93 (MQ=255)
                                               tGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgt                                          <  2:4220126/93‑1 (MQ=255)
                                               tGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgt                                          <  1:5655338/93‑1 (MQ=255)
                                                gTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtg                                         >  2:1856010/1‑93 (MQ=255)
                                                 tATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtgg                                        >  2:2270861/1‑93 (MQ=255)
                                                  aTCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtggt                                       <  1:5358060/93‑1 (MQ=255)
                                                  aTCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtggt                                       >  1:4545337/1‑93 (MQ=255)
                                                   tCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtggtg                                      >  1:3485886/1‑93 (MQ=255)
                                                    cGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGc                                     >  2:566411/1‑93 (MQ=255)
                                                     gCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCt                                    <  2:2064011/93‑1 (MQ=255)
                                                     gCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCt                                    >  2:2408185/1‑93 (MQ=255)
                                                      cTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTg                                   <  1:1754675/93‑1 (MQ=255)
                                                              aCGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTg                           >  1:6684162/1‑93 (MQ=255)
                                                                  gtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       <  1:3681664/93‑1 (MQ=255)
                                                                       gTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCAc                  >  2:2042970/1‑93 (MQ=255)
                                                                       gTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCAc                  >  1:580575/1‑93 (MQ=255)
                                                                       gTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCAc                  >  1:378593/1‑93 (MQ=255)
                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  2:3954082/1‑93 (MQ=255)
                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  2:6411243/1‑93 (MQ=255)
                                                                         tttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACgg                >  2:3264753/1‑93 (MQ=255)
                                                                          ttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGa               <  2:2483467/93‑1 (MQ=255)
                                                                           tCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGaa              <  1:4025238/93‑1 (MQ=255)
                                                                            cGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAAc             <  1:6385762/93‑1 (MQ=255)
                                                                            cGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAAc             >  2:446740/1‑93 (MQ=255)
                                                                             gCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACt            <  1:4922166/93‑1 (MQ=255)
                                                                              cgggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACtt           <  1:4227804/91‑1 (MQ=255)
                                                                               cGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACttt          >  1:4785811/1‑93 (MQ=255)
                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         >  2:6290577/1‑93 (MQ=255)
                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         >  2:3090285/1‑93 (MQ=255)
                                                                                  aTAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATc       <  1:2877457/93‑1 (MQ=255)
                                                                                  aTAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATc       <  1:1577703/93‑1 (MQ=255)
                                                                                  aTAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATc       <  1:1224086/93‑1 (MQ=255)
                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCt      >  2:5982299/1‑93 (MQ=255)
                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCt      <  2:6003324/93‑1 (MQ=255)
                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCt      >  1:6215662/1‑93 (MQ=255)
                                                                                    aaGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCtt     <  1:5230689/93‑1 (MQ=255)
                                                                                    aaGCCCAATCAGGCGCTGGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCtt     <  2:1629150/93‑1 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    >  1:2620637/1‑93 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    >  2:5296878/1‑93 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    >  1:5493589/1‑93 (MQ=255)
                                                                                      gCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAg   <  2:3896580/93‑1 (MQ=255)
                                                                                       cTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGc  <  2:5432348/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110

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