Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,586 2 bp→GC coding (814‑815/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,5860CG100.0% 438.6 / NA 137R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (76/61);  total (76/61)
*NC_0009133,987,5870GC100.0% 437.8 / NA 137R272A (CGG→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (76/61);  total (76/61)

Only 100 of 147 total aligned reads displayed.

TCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGG  >  NC_000913/3987494‑3987672
                                                                                            ||                                                                                     
tCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCTTCTTGCGGCGCACCGCGCCGTTCACTTCGAATAATCTCGGAACCGTccggcc                                                                                        >  2:2942650/1‑93 (MQ=255)
      gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  >  1:4256209/1‑93 (MQ=255)
      gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  <  2:4704345/93‑1 (MQ=255)
      gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                  <  1:1854602/93‑1 (MQ=255)
       ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAc                                                                                 <  2:2362565/93‑1 (MQ=255)
        tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                >  2:4002470/1‑93 (MQ=255)
        tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                >  2:1712591/1‑93 (MQ=255)
        tcttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                >  1:2845725/1‑93 (MQ=255)
         cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCTCCGACCa                                                                               <  1:73206/93‑1 (MQ=255)
             tGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCCAACCGTCCGGCACGCCGACCAgcgc                                                                           >  2:713176/1‑93 (MQ=255)
              gcaaGCGAAACCCCCATTTGTTCGGCATCTTGCGGCGCCCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgca                                                                          <  2:1754062/90‑1 (MQ=255)
              gCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgca                                                                          >  2:167699/1‑93 (MQ=255)
                cAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacg                                                                        >  1:2968391/1‑93 (MQ=255)
                 aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgc                                                                       >  1:3836471/1‑93 (MQ=255)
                 aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgc                                                                       <  1:3463512/93‑1 (MQ=255)
                 aGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcacgc                                                                       >  1:2390996/1‑93 (MQ=255)
                   cGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAg                                                                     <  1:1954369/93‑1 (MQ=255)
                    gAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                    <  2:4252044/93‑1 (MQ=255)
                     aaaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACgccgcc                                                                   >  2:2675159/1‑93 (MQ=255)
                     aaaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGcc                                                                   <  1:1636940/93‑1 (MQ=255)
                      actCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  <  2:2844740/91‑1 (MQ=255)
                      aacccccTTTTGTTCGTCTTCTTGCGGCGCCCCGCGGCGTTCACCGCGAATAATCTCCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  <  2:3440833/90‑1 (MQ=255)
                      aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  >  1:4362960/1‑93 (MQ=255)
                      aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  >  2:962310/1‑93 (MQ=255)
                      aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCa                                                                  >  1:2550930/1‑93 (MQ=255)
                         ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAt                                                               >  2:4127182/1‑93 (MQ=255)
                         ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGAGCACGCAGCCAgct                                                               >  2:2747661/1‑91 (MQ=255)
                          cccATTTGTTCGGCATCTTGCGGCGCCCCGCGGCGTTCCCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACTCAGCCAGAtt                                                              <  2:1521952/93‑1 (MQ=255)
                          cccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAtt                                                              <  1:1217781/93‑1 (MQ=255)
                           ccATTTGTTCGGCATCTTGCGGAGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                             >  2:2991273/1‑93 (MQ=255)
                            cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  1:3418247/1‑93 (MQ=255)
                            cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            <  2:2105169/93‑1 (MQ=255)
                            cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  2:1922197/1‑93 (MQ=255)
                            cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  2:3487350/1‑93 (MQ=255)
                            cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  2:3632868/1‑93 (MQ=255)
                            cATTTGTTCGACATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTc                                                            >  1:2914249/1‑93 (MQ=255)
                                tGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCa                                                        >  1:200857/1‑93 (MQ=255)
                                tGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCa                                                        <  2:2338684/93‑1 (MQ=255)
                                tGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCa                                                        >  1:2063588/1‑93 (MQ=255)
                                 gTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCAt                                                       <  2:4313038/93‑1 (MQ=255)
                                  ttCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                      >  1:4109009/1‑93 (MQ=255)
                                  ttCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                      <  2:4018434/93‑1 (MQ=255)
                                  ttCGGCATCTGGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATc                                                      <  2:2562941/93‑1 (MQ=255)
                                    cGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                    <  2:3887034/93‑1 (MQ=255)
                                      gCATCTTGCGGCGCGCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  <  1:4007153/93‑1 (MQ=255)
                                      gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  <  1:718620/93‑1 (MQ=255)
                                      gCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                  <  1:1348042/93‑1 (MQ=255)
                                        aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAg                                                <  2:1897577/93‑1 (MQ=255)
                                          cTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                              >  2:2726660/1‑93 (MQ=255)
                                          cTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                              <  2:1289044/93‑1 (MQ=255)
                                           ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             <  1:1130733/93‑1 (MQ=255)
                                           ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             <  2:253199/93‑1 (MQ=255)
                                           ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             >  2:1012738/1‑93 (MQ=255)
                                           ttGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTc                                             >  1:620746/1‑93 (MQ=255)
                                              cggcgCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                          >  2:1455550/1‑93 (MQ=255)
                                                gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                        <  1:4720129/93‑1 (MQ=255)
                                                gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                        >  2:3236753/1‑93 (MQ=255)
                                                 cgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCTCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                       <  1:2622627/93‑1 (MQ=255)
                                                    aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  1:1732561/1‑93 (MQ=255)
                                                    aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCt                                    >  1:1732562/1‑93 (MQ=255)
                                                      cgcgGCGTTCCCCGCGAATAATCTGGGAACCGTCCGGCGCGGCGACCAGCCCACGCAGCCAGATTTCGCCATCAATAAGCTCTGCGTAGcgcc                                  >  2:3995818/1‑90 (MQ=255)
                                                       gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCACCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTAcc                                 >  2:1805286/1‑93 (MQ=255)
                                                         ggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCaa                               <  2:4517487/93‑1 (MQ=255)
                                                          gcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAt                              >  2:649461/1‑93 (MQ=255)
                                                          gcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAt                              >  1:879780/1‑93 (MQ=255)
                                                           cgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAtt                             <  2:569083/93‑1 (MQ=255)
                                                           cgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAtt                             <  2:974954/93‑1 (MQ=255)
                                                           cgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAAtt                             <  1:948203/93‑1 (MQ=255)
                                                            gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                            <  2:4255167/93‑1 (MQ=255)
                                                            gTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTg                            >  1:3475398/1‑93 (MQ=255)
                                                             ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATACGCTCGGCGTAGCTACCAATTgg                           >  1:284702/1‑93 (MQ=255)
                                                             ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                           >  1:3190794/1‑93 (MQ=255)
                                                             ttCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTgg                           >  1:554848/1‑93 (MQ=255)
                                                               cACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTTGCa                         >  2:157434/1‑93 (MQ=255)
                                                               cACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                         >  1:1615226/1‑93 (MQ=255)
                                                                cccGCGAATAATCTGCGAACCGTCCGGCGCTCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                        <  1:4537685/92‑1 (MQ=255)
                                                                aCCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                        >  1:1514195/1‑93 (MQ=255)
                                                                 ccGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAcc                       <  1:604440/93‑1 (MQ=255)
                                                                   gcgAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                     <  1:4385275/93‑1 (MQ=255)
                                                                     gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAc                   <  2:3388494/93‑1 (MQ=255)
                                                                     gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGAc                   >  1:3589566/1‑93 (MQ=255)
                                                                     gAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTCGCACCTGAc                   >  1:3589567/1‑93 (MQ=255)
                                                                      aataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACa                  >  2:1432307/1‑93 (MQ=255)
                                                                      aataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACa                  <  2:593882/93‑1 (MQ=255)
                                                                       ataatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACAt                 >  1:3220265/1‑93 (MQ=255)
                                                                        taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                <  2:3338591/93‑1 (MQ=255)
                                                                        taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                <  2:4299038/93‑1 (MQ=255)
                                                                        taatCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATc                >  1:2668524/1‑93 (MQ=255)
                                                                          atCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGATACCAATTGGCACCTGACATCCg              >  2:3227664/1‑93 (MQ=255)
                                                                            cTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGcc            >  2:2150859/1‑93 (MQ=255)
                                                                               cGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc         <  1:3318285/93‑1 (MQ=255)
                                                                                gAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCg        <  1:2020977/93‑1 (MQ=255)
                                                                                 aaCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCga       <  1:1076164/93‑1 (MQ=255)
                                                                                  aCCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag      >  2:1046558/1‑93 (MQ=255)
                                                                                   ccGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     <  2:287476/93‑1 (MQ=255)
                                                                                   ccGTCCGGCGCACCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga     <  2:875686/93‑1 (MQ=255)
                                                                                    cGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGAc    >  1:3871299/1‑93 (MQ=255)
                                                                                     gTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg   >  2:4651312/1‑93 (MQ=255)
                                                                                      tCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg  <  1:4041289/93‑1 (MQ=255)
                                                                                      gccggcgcgCCGACCACCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACgg  <  2:3943296/92‑1 (MQ=255)
                                                                                            ||                                                                                     
TCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCCGGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGG  >  NC_000913/3987494‑3987672

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