Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,807 A→G I198I (ATT→ATC hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,8070AG98.1% 495.8 / ‑2.3 158I198I (ATT→ATChemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (80/75);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 174 total aligned reads displayed.

GTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC  >  NC_000913/3987717‑3987896
                                                                                          |                                                                                         
gtgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAt                                                                                         <  2:1931597/93‑1 (MQ=255)
 tgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtc                                                                                        <  1:2591766/93‑1 (MQ=255)
  gATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTTAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtct                                                                                       >  1:4194937/1‑93 (MQ=255)
  gATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtct                                                                                       >  2:2637049/1‑93 (MQ=255)
  gATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtct                                                                                       >  2:4007342/1‑93 (MQ=255)
   aTTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtctc                                                                                      <  1:1805286/93‑1 (MQ=255)
   aTTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtctc                                                                                      >  1:2791849/1‑93 (MQ=255)
   aTTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtctc                                                                                      <  1:601598/93‑1 (MQ=255)
   aTTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGAtctc                                                                                      <  2:838723/93‑1 (MQ=255)
      cAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCggg                                                                                   <  2:1970428/93‑1 (MQ=255)
        gcgcGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGCGTg                                                                                 >  1:2109661/1‑93 (MQ=255)
             gCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaac                                                                            >  2:4107030/1‑92 (MQ=255)
               aaGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATGCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCaccgc                                                                          <  2:3559691/93‑5 (MQ=255)
                agcagcTCGCGCGTGCGTGAATCATCAAGGCGCCATTCAATCCCCACCGCTCCTTGTCCTCCCGCCGGAAGAGAGATCTCGGGTGGCaacgcg                                                                         <  2:2098397/93‑5 (MQ=255)
                   agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGGGGAAACgcggcg                                                                      >  2:2365155/1‑93 (MQ=255)
                    gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGgcgc                                                                     >  1:700657/1‑93 (MQ=255)
                     cTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGcgcg                                                                    >  1:4322961/1‑93 (MQ=255)
                      tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGa                                                                   >  2:2248375/1‑93 (MQ=255)
                      tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGa                                                                   >  1:4745663/1‑93 (MQ=255)
                       cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                  >  1:1056794/1‑93 (MQ=255)
                       cgcgAGGGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGaa                                                                  >  1:1305934/1‑93 (MQ=255)
                        gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAAt                                                                 >  1:3052074/1‑93 (MQ=255)
                         cgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATa                                                                >  2:3378300/1‑93 (MQ=255)
                          gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATAc                                                               >  2:3333098/1‑93 (MQ=255)
                            gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGt                                                             <  2:3061361/93‑1 (MQ=255)
                              gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAATAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGa                                                           <  2:4348034/93‑1 (MQ=255)
                               cGTGAATCATCAAGGCGGCAGTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                          <  2:2845725/93‑1 (MQ=255)
                                gTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGCAATATCGGGTGGCAACGCGGCGCGAATACGTGACt                                                         >  1:923477/1‑93 (MQ=255)
                                  gAATCATCAAGGTGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCACCGCGGCGCGAATACGTGACTcc                                                       >  2:224302/1‑93 (MQ=255)
                                  gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCCGCGCGAATAAGGGACTcc                                                       >  2:2006683/1‑93 (MQ=255)
                                   aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                      <  1:1655821/93‑1 (MQ=255)
                                    atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                     >  1:2755047/1‑93 (MQ=255)
                                    atcatcAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                     >  1:238907/1‑93 (MQ=255)
                                     tcatcaAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                    <  1:1042199/93‑1 (MQ=255)
                                     tcatcaAGGCGGCATTAAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                    <  1:1442309/93‑1 (MQ=255)
                                       atcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAcc                                                  <  2:2853475/93‑1 (MQ=255)
                                          aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                               >  1:2874470/1‑93 (MQ=255)
                                           aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                              >  1:157771/1‑93 (MQ=255)
                                           aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                              >  1:3829495/1‑93 (MQ=255)
                                            agcggaATTCAATACCCACCGCAACTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAc                                             <  2:349341/92‑1 (MQ=255)
                                              cggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGt                                           <  2:639696/93‑1 (MQ=255)
                                               ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                          >  1:2972251/1‑93 (MQ=255)
                                               ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                          >  1:1926082/1‑93 (MQ=255)
                                                 cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                        >  2:2235312/1‑93 (MQ=255)
                                                 cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                        >  2:1322403/1‑93 (MQ=255)
                                                 cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                        >  2:1322402/1‑93 (MQ=255)
                                                   ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGGCCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAGACGTTTTAg                                      >  2:3701376/1‑93 (MQ=255)
                                                   ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGTCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                      >  2:170774/1‑93 (MQ=255)
                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    <  2:3661652/93‑1 (MQ=255)
                                                     cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                    >  1:2266310/1‑93 (MQ=255)
                                                      aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                   >  1:1027903/1‑93 (MQ=255)
                                                      aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                   <  1:1989106/93‑1 (MQ=255)
                                                       aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCg                                  <  2:119460/93‑1 (MQ=255)
                                                         aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                <  1:2968671/93‑1 (MQ=255)
                                                          cccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                               <  1:4358057/93‑1 (MQ=255)
                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCtc                              >  2:4280777/1‑92 (MQ=255)
                                                           ccACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                              >  1:426561/1‑93 (MQ=255)
                                                            cACCGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                             <  2:4061571/93‑1 (MQ=255)
                                                             cccGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                            <  2:1615226/92‑1 (MQ=255)
                                                             aCCGCACCTTGTCCTACCGCCAGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGGTTTAGTACGGCTACg                            >  2:3403036/1‑93 (MQ=255)
                                                              ccGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                           >  2:1393886/1‑93 (MQ=255)
                                                               cGCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                          <  2:3871299/93‑1 (MQ=255)
                                                               cGCACCGTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                          >  2:4606930/1‑93 (MQ=255)
                                                                gCACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                         >  1:951918/1‑93 (MQ=255)
                                                                 cACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCaa                        >  2:1872579/1‑93 (MQ=255)
                                                                 cACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCaa                        >  2:1482905/1‑93 (MQ=255)
                                                                 cACCTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCaa                        >  2:3204694/1‑93 (MQ=255)
                                                                   ccTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                      <  1:3747302/93‑1 (MQ=255)
                                                                    cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                     >  1:3472491/1‑93 (MQ=255)
                                                                    cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                     <  1:4002470/93‑1 (MQ=255)
                                                                    cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                     <  1:3487350/93‑1 (MQ=255)
                                                                    cTTGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                     <  2:2808463/93‑1 (MQ=255)
                                                                     ttGTCGTACCGCCGGAAGAGAGATCTCGGGTGGCAACGGGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                    <  2:289159/93‑1 (MQ=255)
                                                                     ttGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                    <  2:4338341/93‑1 (MQ=255)
                                                                      tGTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                   >  1:1226790/1‑93 (MQ=255)
                                                                       gTCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                  <  2:3093789/93‑1 (MQ=255)
                                                                        tcatACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                 <  2:2531915/90‑1 (MQ=255)
                                                                        tCCTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                 <  1:2143085/93‑1 (MQ=255)
                                                                         ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                >  2:1651172/1‑93 (MQ=255)
                                                                         ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                <  2:1742532/93‑1 (MQ=255)
                                                                         ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                <  2:1732561/93‑1 (MQ=255)
                                                                         ccTACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                <  2:3589566/93‑1 (MQ=255)
                                                                          cTACCGCCGGAAGAGAGATCTCGGGTGGCAACGAGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg               >  1:4314415/1‑93 (MQ=255)
                                                                           tACCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc              >  2:4084913/1‑93 (MQ=255)
                                                                            aactccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCa             <  1:2571318/89‑1 (MQ=255)
                                                                            aCCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCa             >  1:82913/1‑93 (MQ=255)
                                                                            aCCGCCGGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCa             >  2:2502659/1‑93 (MQ=255)
                                                                              cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc           <  1:4651312/93‑1 (MQ=255)
                                                                              cgccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc           <  2:2550930/93‑1 (MQ=255)
                                                                               gccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg          <  1:2895429/93‑1 (MQ=255)
                                                                               gccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg          >  2:739746/1‑93 (MQ=255)
                                                                                ccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt         >  2:4111970/1‑93 (MQ=255)
                                                                                ccgGAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt         <  1:476283/93‑1 (MQ=255)
                                                                                  ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt       >  1:3180852/1‑93 (MQ=255)
                                                                                  ggAAGAGAGATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGCCTACGGCAGGAATGATGGCATCGTAt       >  1:4027849/1‑93 (MQ=255)
                                                                                     agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAGGTTTTAGTCCGGCTACGTCACGAATAAAGGCATCGTATCCg    >  2:357394/1‑93 (MQ=255)
                                                                                     agagagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg    >  2:1677165/1‑93 (MQ=255)
                                                                                      gagagaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc   >  2:4390999/1‑93 (MQ=255)
                                                                                       agagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc  >  2:579734/1‑93 (MQ=255)
                                                                                       agagATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc  <  1:3373805/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC  >  NC_000913/3987717‑3987896

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