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% 477.9 / NA 133D29N (GAT→AAT) hemAHemA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (81/52);  total (81/52)

Only 100 of 137 total aligned reads displayed.

GCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCGATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGC  >  NC_000913/1262931‑1263110
                                                                                          |                                                                                         
gCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAAt                                                                                         <  2:2688767/93‑1 (MQ=255)
gCAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAAt                                                                                         >  1:2542027/1‑93 (MQ=255)
 cAGACATTACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc                                                                                        <  2:3835955/93‑1 (MQ=255)
 cAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc                                                                                        <  1:764683/93‑1 (MQ=255)
 cAGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATc                                                                                        >  2:1941494/1‑93 (MQ=255)
  aGACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCa                                                                                       <  1:650782/93‑1 (MQ=255)
   gACATGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAg                                                                                      >  2:2003641/1‑93 (MQ=255)
    aCATGACCCTCTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAgg                                                                                     <  2:3464812/93‑1 (MQ=255)
      aTGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGgcg                                                                                   >  1:3268588/1‑93 (MQ=255)
      aTGACCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGccg                                                                                   >  2:50446/1‑93 (MQ=255)
         aCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtt                                                                                >  1:1123352/1‑93 (MQ=255)
         aCCCTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCtc                                                                                >  1:3375276/1‑92 (MQ=255)
          cccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTg                                                                               >  1:3642286/1‑93 (MQ=255)
          cccTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTg                                                                               >  2:3929304/1‑93 (MQ=255)
            cTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGAc                                                                             <  2:3440355/93‑1 (MQ=255)
            cTTTTAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGAc                                                                             >  1:2296025/1‑93 (MQ=255)
             ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            >  1:1808229/1‑93 (MQ=255)
             ttttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACa                                                                            >  2:4531080/1‑93 (MQ=255)
              tttAGCACTCGGTATCAACCATAAACCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAg                                                                           <  1:1244323/93‑1 (MQ=255)
              tttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAg                                                                           <  1:2145140/93‑1 (MQ=255)
              tttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAg                                                                           <  1:1244324/93‑1 (MQ=255)
              tttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAg                                                                           >  1:4428973/1‑93 (MQ=255)
               ttAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGc                                                                          >  1:475966/1‑93 (MQ=255)
                tAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGcc                                                                         >  1:3398387/1‑93 (MQ=255)
                tAGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGcc                                                                         >  2:1310755/1‑93 (MQ=255)
                 aGCACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGAGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCt                                                                        <  1:3541388/93‑1 (MQ=255)
                   cACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGc                                                                      <  1:2552290/93‑1 (MQ=255)
                   cACTCGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGc                                                                      <  1:2918959/93‑1 (MQ=255)
                       cGGTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgc                                                                  <  1:631305/93‑1 (MQ=255)
                        ggTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgcg                                                                 >  1:1747235/1‑93 (MQ=255)
                         gTATCAACCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTgcgc                                                                >  1:3282758/1‑93 (MQ=255)
                              aaCCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCg                                                           >  1:1125114/1‑93 (MQ=255)
                               aCCATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGa                                                          >  1:4076021/1‑93 (MQ=255)
                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGag                                                         >  1:3348630/1‑92 (MQ=255)
                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGAt                                                         >  1:2204577/1‑93 (MQ=255)
                                ccATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGAt                                                         >  1:404347/1‑93 (MQ=255)
                                 cATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATg                                                        >  1:4361312/1‑93 (MQ=255)
                                 cATAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATg                                                        >  1:1817950/1‑93 (MQ=255)
                                   tAAAACGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGt                                                      <  2:3713397/93‑1 (MQ=255)
                                      aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   <  1:178800/93‑1 (MQ=255)
                                      aaCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCa                                                   <  1:5258/93‑1 (MQ=255)
                                       aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  >  1:3392574/1‑93 (MQ=255)
                                       aCGGCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAg                                                  >  2:1738023/1‑93 (MQ=255)
                                          gCACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCGTGCTTGCGCAGCCGATGGTGCAGggc                                               >  2:2829832/1‑93 (MQ=255)
                                           cACCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcg                                              >  1:1559177/1‑93 (MQ=255)
                                            aCCTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGggcgg                                             >  2:824968/1‑93 (MQ=255)
                                              cTGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGgcggcg                                           <  1:997698/93‑1 (MQ=255)
                                               tGTATCGCTGCGAGAGCGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgt                                          <  2:3669205/93‑1 (MQ=255)
                                               tGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGgcggcgg                                          >  1:3789192/1‑92 (MQ=255)
                                               tGTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgt                                          <  2:2319561/93‑1 (MQ=255)
                                                gTATCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGGTGGTGCAGGGCGGCgtg                                         >  2:3907584/1‑93 (MQ=255)
                                                  aTCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtggt                                       <  2:1001757/93‑1 (MQ=255)
                                                   tCGCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCgtggtg                                      >  2:488642/1‑93 (MQ=255)
                                                     gCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCt                                    >  1:4311372/1‑93 (MQ=255)
                                                     gCTGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCt                                    >  2:3653643/1‑93 (MQ=255)
                                                       tGCGAGGTCTTGTATCTTTGTCTCCGGATAAGCTCAATCAGGCTCTGGACTGCCAGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  <  2:3584647/93‑1 (MQ=255)
                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  <  2:2024085/93‑1 (MQ=255)
                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  <  2:2848389/93‑1 (MQ=255)
                                                       tGCGAGAACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGt                                  >  1:2650471/1‑93 (MQ=255)
                                                          gagaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGa                               >  1:1215949/1‑93 (MQ=255)
                                                           agaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGAc                              >  2:2086772/1‑93 (MQ=255)
                                                            gaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACg                             >  1:854357/1‑93 (MQ=255)
                                                            gaACGTGTATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACg                             >  1:4236275/1‑93 (MQ=255)
                                                                gtgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCa                         >  1:2926587/1‑93 (MQ=255)
                                                                gtgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCa                         >  2:623728/1‑93 (MQ=255)
                                                                gtgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCa                         >  1:2698289/1‑93 (MQ=255)
                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGTaa                        <  1:340452/93‑1 (MQ=255)
                                                                 tgtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCaa                        >  1:2662840/1‑93 (MQ=255)
                                                                  gtATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       <  1:4133039/93‑1 (MQ=255)
                                                                  gtATCGTGTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAc                       <  1:3652372/93‑1 (MQ=255)
                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAcc                      >  2:2494267/1‑93 (MQ=255)
                                                                   tATCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAcc                      <  1:225201/93‑1 (MQ=255)
                                                                   gaTCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAAcc                      <  2:1925654/92‑1 (MQ=255)
                                                                    aTCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCg                     >  1:1940119/1‑93 (MQ=255)
                                                                     tCGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGc                    >  2:2232473/1‑93 (MQ=255)
                                                                      cGTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCa                   <  2:1370100/93‑1 (MQ=255)
                                                                      cGTTTTCCCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCa                   <  1:3834581/93‑1 (MQ=255)
                                                                       gTTTTCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCAc                  >  1:3321764/1‑93 (MQ=255)
                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  2:4321516/1‑93 (MQ=255)
                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  2:3770420/1‑93 (MQ=255)
                                                                        ttttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACg                 >  2:2401710/1‑93 (MQ=255)
                                                                          ttCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGa               >  1:402931/1‑93 (MQ=255)
                                                                           tCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGaa              <  1:1827149/93‑1 (MQ=255)
                                                                           tCGCCGGATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGaa              <  2:4674066/93‑1 (MQ=255)
                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         >  2:4216119/1‑93 (MQ=255)
                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         <  1:2056162/93‑1 (MQ=255)
                                                                                ggATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTa         >  1:2169438/1‑93 (MQ=255)
                                                                                 gATAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTAt        >  1:2232792/1‑93 (MQ=255)
                                                                                  aTAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATc       <  1:1152216/93‑1 (MQ=255)
                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCATGGCGGCGTGGTGCTGTCGACGTGCAAGCGCACGGAACTTTATCt      >  2:2493497/1‑93 (MQ=255)
                                                                                   tAAGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCt      >  1:4523514/1‑93 (MQ=255)
                                                                                    aaGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCtt     >  1:80413/1‑93 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTTACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    <  1:3800520/93‑1 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    <  1:1629422/93‑1 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    <  1:2011369/93‑1 (MQ=255)
                                                                                     aGCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTa    <  1:1630339/93‑1 (MQ=255)
                                                                                      gCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGGACGGAACTTTATCTTAg   >  2:3302400/1‑93 (MQ=255)
                                                                                      gCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAg   >  2:319095/1‑93 (MQ=255)
                                                                                      gCTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAg   <  2:3113641/93‑1 (MQ=255)
                                                                                       cTCAATCAGGCGCTTGACAGCCTGCTTGCGCAGCCGATGGTGCAGGGCGGCGTGGTGCTGTCGACGTGCAACCGCACGGAACTTTATCTTAGc  <  1:958378/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 < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: