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,7410TC98.3% 558.3 / ‑2.0 178V313V (GTA→GTGhemXuroporphyrinogen III methylase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (87/88);  total (88/90)
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.99e-01

Only 100 of 186 total aligned reads displayed.

AGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985652‑3985830
                                                                                         |                                                                                         
 gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       <  2:2588252/93‑1 (MQ=255)
 gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       <  1:1391682/93‑1 (MQ=255)
 gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                       <  1:4540097/93‑1 (MQ=255)
  aTCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCa                                                                                      <  1:1771678/93‑1 (MQ=255)
   tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     >  2:3039152/1‑93 (MQ=255)
   tCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAg                                                                                     <  2:3341529/93‑1 (MQ=255)
     cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                   <  1:3026260/93‑1 (MQ=255)
      aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  >  1:4472221/1‑93 (MQ=255)
      aTCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTg                                                                                  <  1:779006/93‑1 (MQ=255)
        cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                                <  2:4133534/93‑1 (MQ=255)
         gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgag                                                                               >  2:1736936/1‑93 (MQ=255)
         gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgag                                                                               <  2:3072202/93‑1 (MQ=255)
         gagaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgag                                                                               >  2:4409341/1‑93 (MQ=255)
          agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgcga                                                                              >  1:4317696/1‑93 (MQ=255)
          agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              >  1:4619356/1‑93 (MQ=255)
          agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                              >  2:350331/1‑93 (MQ=255)
           gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             <  1:4218848/93‑1 (MQ=255)
           gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                             <  1:2095421/93‑1 (MQ=255)
            atatTTTGCTGGCTTAACTTGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACg                                                                            >  1:999877/1‑93 (MQ=255)
             tatTTTGCTGGCTTATCTGGTCCACCTTGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           <  2:811748/93‑1 (MQ=255)
             tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           <  2:1880182/93‑1 (MQ=255)
             tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                           >  2:162008/1‑93 (MQ=255)
              atTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGtt                                                                          <  2:3545638/93‑1 (MQ=255)
               ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                         >  2:3485546/1‑93 (MQ=255)
               ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                         <  1:1122462/93‑1 (MQ=255)
                 ttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtctc                                                                       <  1:1898364/93‑1 (MQ=255)
                  tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                      >  1:2372880/1‑93 (MQ=255)
                  tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                      <  2:1658989/93‑1 (MQ=255)
                   gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     <  2:1836001/93‑1 (MQ=255)
                   gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     <  1:4494209/93‑1 (MQ=255)
                   gCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCa                                                                     >  2:894296/1‑93 (MQ=255)
                    cTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAg                                                                    >  2:4482059/1‑93 (MQ=255)
                    cTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAg                                                                    <  2:1188755/93‑1 (MQ=255)
                      ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgcg                                                                  >  2:3673816/1‑93 (MQ=255)
                        cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCACCAGCGcc                                                                >  2:225346/1‑93 (MQ=255)
                         ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                               <  2:3883523/93‑1 (MQ=255)
                           aaCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTgg                                                             >  1:3327741/1‑93 (MQ=255)
                            aCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGc                                                            >  2:1152870/1‑93 (MQ=255)
                             cTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCg                                                           >  2:3326688/1‑93 (MQ=255)
                             cTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCg                                                           <  1:2409294/93‑1 (MQ=255)
                             cTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCg                                                           <  1:906905/93‑1 (MQ=255)
                              tGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGa                                                          <  1:2062458/93‑1 (MQ=255)
                                 tCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATaa                                                       <  2:1255302/93‑1 (MQ=255)
                                  ccACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAg                                                      >  1:1351148/1‑93 (MQ=255)
                                  ccACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAg                                                      <  1:3836351/93‑1 (MQ=255)
                                    cccTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    <  1:475775/92‑1 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    <  1:2701776/93‑1 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    >  2:4283698/1‑93 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    >  2:4201899/1‑93 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    <  2:4709600/93‑1 (MQ=255)
                                    aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                    <  1:3299138/93‑1 (MQ=255)
                                       tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  1:3967603/93‑1 (MQ=255)
                                       tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 >  2:1139434/1‑93 (MQ=255)
                                       tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                 <  1:3488037/93‑1 (MQ=255)
                                         gtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTc                                               <  2:439131/93‑1 (MQ=255)
                                          tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                              >  2:2783451/1‑93 (MQ=255)
                                          tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                              >  1:2018766/1‑93 (MQ=255)
                                           cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             <  1:1979629/93‑1 (MQ=255)
                                           cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             >  2:4514675/1‑93 (MQ=255)
                                           cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             <  1:2315512/93‑1 (MQ=255)
                                           cgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCt                                             <  2:3550365/93‑1 (MQ=255)
                                             aGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGGCGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTgg                                           >  1:711997/1‑93 (MQ=255)
                                              ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          <  1:2257727/93‑1 (MQ=255)
                                              ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          >  1:225528/1‑93 (MQ=255)
                                              ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          >  1:3515396/1‑93 (MQ=255)
                                              agAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                          <  1:1673229/92‑1 (MQ=255)
                                                  cGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcg                                      >  1:541595/1‑93 (MQ=255)
                                                   gCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcgc                                     <  1:2153528/93‑1 (MQ=255)
                                                   gCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcgc                                     >  1:3586102/1‑93 (MQ=255)
                                                     tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCgg                                   >  2:4081098/1‑93 (MQ=255)
                                                     tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCgg                                   <  1:1307162/93‑1 (MQ=255)
                                                       tggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTa                                 >  2:1997582/1‑93 (MQ=255)
                                                       tggtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTa                                 <  2:2440739/93‑1 (MQ=255)
                                                         gtggtggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACa                               <  2:2150521/93‑1 (MQ=255)
                                                             tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                           <  1:2462408/93‑1 (MQ=255)
                                                             tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                           <  2:150686/93‑1 (MQ=255)
                                                              ggtatcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                          >  2:3089232/4‑93 (MQ=255)
                                                              gatatcatcAGTATCTTAGTAAGCACGCATGCAGGTGGAGACGTTCTCAAGCGCCTGTCGATAAAACTGTGCATGGTGCCGCGGTAGAGCTtg                          <  2:1998213/90‑1 (MQ=255)
                                                                 atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       <  2:4649615/93‑1 (MQ=255)
                                                                 atcatcAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgc                       >  1:2648295/1‑93 (MQ=255)
                                                                  tcatcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgct                      >  1:3993715/1‑93 (MQ=255)
                                                                  tcatcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgct                      >  2:1821725/1‑93 (MQ=255)
                                                                    atcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctgc                    >  1:1190309/1‑93 (MQ=255)
                                                                    atcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctgc                    >  1:3265488/1‑93 (MQ=255)
                                                                       aGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                 >  2:2691026/1‑93 (MQ=255)
                                                                       aGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                 >  2:2691025/1‑93 (MQ=255)
                                                                        gTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAcc                >  2:1444008/1‑93 (MQ=255)
                                                                         tATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACca               <  2:4107699/93‑1 (MQ=255)
                                                                         tATCGTAGTAAGCACGCTCCCAGGTGGAGACGGTCTCCAGCGCCTGGCGATAAGTCGCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACca               <  1:4137500/93‑1 (MQ=255)
                                                                             gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag           >  2:1341742/1‑93 (MQ=255)
                                                                             gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag           >  1:1543273/1‑93 (MQ=255)
                                                                             gtagtaAGCACGCACCCAGGTGGAGACGTCCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag           >  1:2870155/1‑93 (MQ=255)
                                                                                 taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga       >  2:3939651/1‑93 (MQ=255)
                                                                                  aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      >  1:4380862/1‑93 (MQ=255)
                                                                                  aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag      <  1:1082516/93‑1 (MQ=255)
                                                                                   aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga     <  1:3085280/93‑1 (MQ=255)
                                                                                   aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga     <  1:1673111/93‑1 (MQ=255)
                                                                                    gcacgcacCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGagag    <  2:1349351/93‑1 (MQ=255)
                                                                                    gcacgcacCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGagag    <  2:203921/93‑1 (MQ=255)
                                                                                      aCGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCg  >  2:3173050/1‑93 (MQ=255)
                                                                                         |                                                                                         
AGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGCG  >  NC_000913/3985652‑3985830

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