Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,986,029 A→G R217R (CGT→CGC hemX ← uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,986,0290AG99.6% 879.6 / NA 253R217R (CGT→CGChemXuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (129/123);  minor base T (0/1);  total (129/124)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

TGCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGAT  >  NC_000913/3985940‑3986118
                                                                                         |                                                                                         
tGCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGgcgc                                                                                        >  2:6140875/1‑93 (MQ=255)
 gCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCa                                                                                       <  2:2470865/93‑1 (MQ=255)
  cAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAg                                                                                      >  2:481514/1‑93 (MQ=255)
  cAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAg                                                                                      <  1:3614605/93‑1 (MQ=255)
  cAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAg                                                                                      <  1:6247691/93‑1 (MQ=255)
  cAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAg                                                                                      <  1:727891/93‑1 (MQ=255)
  cAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAg                                                                                      <  1:1268190/93‑1 (MQ=255)
   aGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAgg                                                                                     >  2:794634/1‑93 (MQ=255)
   aGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAgg                                                                                     >  1:984047/1‑93 (MQ=255)
    gATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGt                                                                                    >  1:1708322/1‑93 (MQ=255)
    gATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGt                                                                                    >  1:4170960/1‑93 (MQ=255)
     aTTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCCGGtt                                                                                   >  2:5612112/1‑93 (MQ=255)
      ttGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTa                                                                                  >  2:3978535/1‑93 (MQ=255)
       tGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTAt                                                                                 >  1:4887891/1‑93 (MQ=255)
        gATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATc                                                                                >  2:6196741/1‑93 (MQ=255)
        gATACGCCATTCGCTGATGGAACTGGAAAGCTCTACACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATc                                                                                <  2:520903/93‑1 (MQ=255)
         aTACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCt                                                                               >  1:5985779/1‑93 (MQ=255)
             gCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCt                                                                           <  2:2625821/93‑1 (MQ=255)
             gCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCt                                                                           >  2:5137557/1‑93 (MQ=255)
              ccATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTg                                                                          >  2:5965439/1‑93 (MQ=255)
              ccATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTg                                                                          >  2:6412120/1‑93 (MQ=255)
               cATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGa                                                                         >  1:2494524/1‑93 (MQ=255)
               cATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGa                                                                         <  2:4402668/93‑1 (MQ=255)
                aTTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGAt                                                                        >  2:3610883/1‑93 (MQ=255)
                aTTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGAt                                                                        >  1:608973/1‑93 (MQ=255)
                 ttCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGAtt                                                                       >  1:4652325/1‑93 (MQ=255)
                   cGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTg                                                                     >  2:1349115/1‑93 (MQ=255)
                   cGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTg                                                                     >  2:5870658/1‑93 (MQ=255)
                    gCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGa                                                                    >  2:3520830/1‑93 (MQ=255)
                    gCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGa                                                                    >  2:6300597/1‑93 (MQ=255)
                    gCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGa                                                                    <  2:2036151/93‑1 (MQ=255)
                     cTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGaa                                                                   >  2:5184325/1‑93 (MQ=255)
                       gATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAg                                                                 >  1:3452568/1‑93 (MQ=255)
                       gATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAg                                                                 >  2:341682/1‑93 (MQ=255)
                       gATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAg                                                                 >  1:6285284/1‑93 (MQ=255)
                        aTGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGc                                                                <  1:4965663/93‑1 (MQ=255)
                          ggAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTg                                                              >  1:2960423/1‑93 (MQ=255)
                           gAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGa                                                             <  1:590675/93‑1 (MQ=255)
                            aaCTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGAt                                                            <  1:1361870/93‑1 (MQ=255)
                            aaCTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGAt                                                            <  1:4361372/93‑1 (MQ=255)
                               tGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTaa                                                         >  2:3266162/1‑93 (MQ=255)
                               tGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTaa                                                         >  1:5713137/1‑93 (MQ=255)
                               tGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTaa                                                         <  2:4713890/93‑1 (MQ=255)
                                  aaaGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCt                                                      <  2:752882/93‑1 (MQ=255)
                                   aaGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCtt                                                     >  2:2968228/1‑93 (MQ=255)
                                   aaGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCtt                                                     >  1:1833258/1‑93 (MQ=255)
                                   aaGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCtt                                                     >  1:3823703/1‑93 (MQ=255)
                                    aGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTa                                                    <  1:241464/93‑1 (MQ=255)
                                    aGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTa                                                    >  1:3882757/1‑93 (MQ=255)
                                    aGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTa                                                    >  2:5574390/1‑93 (MQ=255)
                                    aGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTa                                                    >  1:1000722/1‑93 (MQ=255)
                                     gCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTaa                                                   <  2:846945/93‑1 (MQ=255)
                                       tctTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAgg                                                 <  2:1182529/93‑1 (MQ=255)
                                       tctTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAgg                                                 >  2:5261407/1‑93 (MQ=255)
                                        ctTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGga                                                <  2:2048861/93‑1 (MQ=255)
                                        ctTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGga                                                >  1:2765275/1‑93 (MQ=255)
                                         ttCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGgat                                               <  1:4760607/93‑1 (MQ=255)
                                          tCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGgatg                                              <  2:4040771/93‑1 (MQ=255)
                                           cACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGgatga                                             <  2:4957988/93‑1 (MQ=255)
                                            aCCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGgatgat                                            <  1:3356652/93‑1 (MQ=255)
                                            aCCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGgatgat                                            <  1:2628500/93‑1 (MQ=255)
                                             ccGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGatgatg                                           >  2:5962761/1‑93 (MQ=255)
                                               gTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGcc                                         >  2:5851938/1‑93 (MQ=255)
                                                tCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCg                                        >  1:5465158/1‑93 (MQ=255)
                                                 cTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGt                                       >  1:6403592/1‑93 (MQ=255)
                                                  tGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTc                                      <  2:5676394/93‑1 (MQ=255)
                                                     aTCCATCGGCGAACCATCGCTGTCATTATCGGGCAGGCGCGCGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATa                                   >  1:873586/1‑93 (MQ=255)
                                                       ccATCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAAt                                 <  1:6074151/93‑1 (MQ=255)
                                                         aTCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTAAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATcc                               <  1:2158362/93‑1 (MQ=255)
                                                          tCGGCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCa                              <  1:2832345/93‑1 (MQ=255)
                                                            ggCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCAcc                            >  2:4586633/1‑93 (MQ=255)
                                                            ggCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCAcc                            <  2:1902869/93‑1 (MQ=255)
                                                            ggCGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCAcc                            >  2:5994833/1‑93 (MQ=255)
                                                              cGAACCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTg                          >  2:2351517/1‑93 (MQ=255)
                                                                aaCCATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCg                        >  1:2826357/1‑93 (MQ=255)
                                                                 aCCATCGCTGTCATTATCGGCCAGGCTCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGa                       >  2:5728217/1‑93 (MQ=255)
                                                                   cATCGCTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATa                     <  1:153410/93‑1 (MQ=255)
                                                                      cGCTGTCATTATCGGCCAGTCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTg                  <  1:4455038/93‑1 (MQ=255)
                                                                        cTGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCa                >  1:3160468/1‑93 (MQ=255)
                                                                         tGTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAg               >  2:3314532/1‑93 (MQ=255)
                                                                          gTCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGa              <  2:3147104/93‑1 (MQ=255)
                                                                           tCATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATATTGCAGaa             >  1:1110054/1‑93 (MQ=255)
                                                                            cATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGaaa            >  2:5505281/1‑93 (MQ=255)
                                                                            cATTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGaaa            >  2:2825070/1‑93 (MQ=255)
                                                                             aTTATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAg           <  2:223327/93‑1 (MQ=255)
                                                                              ttATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAgg          >  2:1836998/1‑93 (MQ=255)
                                                                              ttATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAgg          <  1:2320267/93‑1 (MQ=255)
                                                                              ttATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTAATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAgg          >  1:5244423/1‑93 (MQ=255)
                                                                               tATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGc         >  2:5312321/1‑93 (MQ=255)
                                                                               tATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGc         <  1:3780516/93‑1 (MQ=255)
                                                                               tATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGc         <  2:5040647/93‑1 (MQ=255)
                                                                               tATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGc         >  1:1308059/1‑93 (MQ=255)
                                                                               tATCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGc         <  1:192732/93‑1 (MQ=255)
                                                                                aTCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCt        <  2:1040324/93‑1 (MQ=255)
                                                                                 tCGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTg       <  1:3947195/93‑1 (MQ=255)
                                                                                  cGGCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTgg      <  1:3326764/93‑1 (MQ=255)
                                                                                   ggCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGc     <  1:5686061/93‑1 (MQ=255)
                                                                                   ggCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGc     >  1:4305185/1‑93 (MQ=255)
                                                                                    gCCAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCg    >  1:4814637/1‑93 (MQ=255)
                                                                                      cAGGCGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGat  <  1:5322754/93‑1 (MQ=255)
                                                                                         |                                                                                         
TGCAGATTGATACGCCATTCGCTGATGGAACTGGAAAGCTCTTCACCGTCTGAATCCATCGGCGAACCATCGCTGTCATTATCGGCCAGACGCAGGTTATCTACCTGATTTGAAAGCTGATTAAGCTTAAGGATGATGCCGTCATAATCCACCTGCGATACTGCAGAAAGGCTGGCGAT  >  NC_000913/3985940‑3986118

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