Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,985,741 T→C V313V (GTA→GTG hemX ← uroporphyrinogen III methylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,985,7410TC99.8% 1300.4 / ‑5.1 407V313V (GTA→GTGhemXuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (223/183);  total (224/183)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

Only 100 of 429 total aligned reads displayed.

AGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985652‑3985830
                                                                                         |                                                                                         
aGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAcc                                                                                        >  2:5393125/1‑93 (MQ=255)
 gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       >  2:8862453/1‑93 (MQ=255)
  aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCa                                                                                      >  1:3909353/1‑93 (MQ=255)
   tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     >  1:6617275/1‑93 (MQ=255)
   tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     <  1:7172564/93‑1 (MQ=255)
    ccATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAgg                                                                                    >  1:4590066/1‑93 (MQ=255)
     cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   <  2:1545843/93‑1 (MQ=255)
     cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   >  2:5045799/1‑93 (MQ=255)
     cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   <  2:9715504/93‑1 (MQ=255)
      aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  <  2:8532203/93‑1 (MQ=255)
      aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  1:5186109/1‑93 (MQ=255)
      aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGGGGTGGCATCATCAGGATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  1:7943401/1‑93 (MQ=255)
       tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTgg                                                                                 >  1:7523501/1‑93 (MQ=255)
       tCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTgg                                                                                 >  2:8207477/1‑93 (MQ=255)
        cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTggg                                                                                >  1:7773386/1‑92 (MQ=255)
        cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                                <  1:9775407/93‑1 (MQ=255)
        cGAGATATTTTGCTGGCTGAACCGGTCCCCCTCGTCGAGGAACCCTTTGGTGGTGGCATCAGCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                                <  1:6443171/93‑1 (MQ=255)
          agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              <  2:560384/93‑1 (MQ=255)
           gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             >  1:8454310/1‑93 (MQ=255)
            atatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGAACCCAGGGGGAGACg                                                                            >  2:403788/1‑93 (MQ=255)
              atTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGtt                                                                          <  1:93823/93‑1 (MQ=255)
               ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                         >  2:3745174/1‑93 (MQ=255)
                  tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                      <  1:2043377/93‑1 (MQ=255)
                  tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                      <  1:8557673/93‑1 (MQ=255)
                   gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     <  2:7310253/93‑1 (MQ=255)
                     tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                   <  1:1346223/93‑1 (MQ=255)
                      ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgcg                                                                  >  2:44242/1‑93 (MQ=255)
                       gcctAACTGGTCCCCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgcgc                                                                 <  1:2710381/90‑1 (MQ=255)
                        cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                                <  1:8276720/93‑1 (MQ=255)
                        cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                                <  2:4291673/93‑1 (MQ=255)
                         ttAACTGGTCCCCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               <  1:3102907/93‑1 (MQ=255)
                          tAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTg                                                              >  2:6619219/1‑93 (MQ=255)
                          tAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTg                                                              >  1:3347670/1‑93 (MQ=255)
                           aaCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTgg                                                             >  1:8159475/1‑93 (MQ=255)
                            aCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGc                                                            >  2:7880420/1‑93 (MQ=255)
                               ggTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                         >  2:1979440/1‑93 (MQ=255)
                               ggTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                         <  1:5329976/93‑1 (MQ=255)
                                gtgcACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                        <  1:4760283/90‑1 (MQ=255)
                                gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                        <  1:6849641/93‑1 (MQ=255)
                                gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                        <  1:1273986/93‑1 (MQ=255)
                                gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                        <  1:2005821/93‑1 (MQ=255)
                                 tCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATaa                                                       >  1:432284/1‑93 (MQ=255)
                                   cACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGt                                                     >  2:2284492/1‑93 (MQ=255)
                                   cACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGt                                                     >  1:8891667/1‑93 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGGCGTTCTCCAGCGCCTGCCGATAAGtc                                                    >  1:2281173/1‑93 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    >  1:5723513/1‑93 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    >  2:1247833/1‑93 (MQ=255)
                                     ccTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtct                                                   >  1:1517952/1‑93 (MQ=255)
                                      cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                  >  2:9049642/1‑93 (MQ=255)
                                       ccgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  1:5781284/92‑1 (MQ=255)
                                       ccgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  2:2888401/92‑1 (MQ=255)
                                        cgtcgAGGAACGCTTTGGTGGTGGCATCATCAGAATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGCCTCtt                                                >  2:5183285/1‑93 (MQ=255)
                                           cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             >  2:7410965/1‑93 (MQ=255)
                                            gAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                            >  2:8353793/1‑93 (MQ=255)
                                             aGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           >  2:3528032/1‑93 (MQ=255)
                                                aaCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGa                                        <  2:3637514/93‑1 (MQ=255)
                                                aaCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGa                                        <  2:9946172/93‑1 (MQ=255)
                                                 aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                       >  1:2441965/1‑93 (MQ=255)
                                                 aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                       >  1:9144967/1‑93 (MQ=255)
                                                  cGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcg                                      >  2:9296072/1‑93 (MQ=255)
                                                  cGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcg                                      >  2:3294178/1‑93 (MQ=255)
                                                     tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTGCAGCGCCTGGCGATAAGTATCTTCCTGGTGACGCgg                                   >  1:10034804/1‑93 (MQ=255)
                                                      ttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGt                                  >  2:6735907/1‑93 (MQ=255)
                                                      ttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGt                                  >  2:4629190/1‑93 (MQ=255)
                                                      ttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTATCTTCCTGGTGACGCGGt                                  >  2:6905427/1‑93 (MQ=255)
                                                       tggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTa                                 <  1:2857011/93‑1 (MQ=255)
                                                       tggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGAGACGCGGTa                                 >  1:2559276/1‑93 (MQ=255)
                                                        ggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTAc                                <  2:2461615/93‑1 (MQ=255)
                                                         gtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACa                               <  1:6565710/93‑1 (MQ=255)
                                                          tggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAg                              >  2:9584662/1‑93 (MQ=255)
                                                             tggCATCATCAGTATCGTAGTAAGCACGCTCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                           <  1:604394/93‑1 (MQ=255)
                                                             tggCATCATCAGTATCGTAGTAAGCACGCCCCCAGGTCGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                           <  1:1887943/93‑1 (MQ=255)
                                                              ggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                          >  1:5602305/1‑93 (MQ=255)
                                                               gCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                         >  2:1066279/1‑93 (MQ=255)
                                                                catcatCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgtg                        >  2:737625/1‑93 (MQ=255)
                                                                 atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       >  2:4814085/1‑93 (MQ=255)
                                                                 atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       >  1:3719021/1‑93 (MQ=255)
                                                                 atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       >  2:710868/1‑93 (MQ=255)
                                                                        gTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAcc                >  2:4886703/1‑93 (MQ=255)
                                                                         tATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAccc               >  1:8734260/1‑92 (MQ=255)
                                                                          aTCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcag              >  1:7002448/1‑93 (MQ=255)
                                                                          aTCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcag              <  2:4185268/93‑1 (MQ=255)
                                                                          aTCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcag              <  1:2544598/93‑1 (MQ=255)
                                                                           tCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAAGCGGTACAGCTTGTGCTGCGACcagc             >  2:3425519/1‑93 (MQ=255)
                                                                           tCGTAGTAAGCACGCACCCAGGTGGAGACGTTCGCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagc             <  1:5872627/93‑1 (MQ=255)
                                                                             gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag           <  2:4578214/93‑1 (MQ=255)
                                                                             gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag           >  1:9721898/1‑93 (MQ=255)
                                                                              tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          >  1:442915/1‑93 (MQ=255)
                                                                              tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          >  2:620437/1‑93 (MQ=255)
                                                                              tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          <  2:5781707/93‑1 (MQ=255)
                                                                              tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg          >  1:7361433/1‑93 (MQ=255)
                                                                              aagtaAGCCCGCACCCAGGAGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGCGACGCGGTACAGCTTGTGCTGCGAGCAGCAgg          >  2:3053291/2‑93 (MQ=255)
                                                                               agtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc         <  1:4921582/93‑1 (MQ=255)
                                                                               agtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGc         >  1:6253309/1‑93 (MQ=255)
                                                                                 taAGCATGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       <  1:7612411/93‑1 (MQ=255)
                                                                                 taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       <  1:4087306/93‑1 (MQ=255)
                                                                                  aaGCCCGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      <  2:1859167/93‑1 (MQ=255)
                                                                                  aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      <  2:1465367/93‑1 (MQ=255)
                                                                                  aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      >  2:1735723/1‑93 (MQ=255)
                                                                                   aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGATGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga     <  2:6098439/93‑1 (MQ=255)
                                                                                         |                                                                                         
AGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985652‑3985830

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