Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,822 2 bp→CG coding (578‑579/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,8220GC100.0% 1940.4 / NA 564G193A (GGC→GCGhemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (275/289);  total (275/289)
*NC_0009133,987,8230CG99.8% 1934.0 / NA 561G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (271/289);  total (271/290)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 583 total aligned reads displayed.

CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG  >  NC_000913/3987735‑3987909
                                                                                        ||                                                                                      
cagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgc                                                                                    <  1:6914188/93‑1 (MQ=255)
cagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgc                                                                                    <  1:2336102/93‑1 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  <  1:8238509/93‑1 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  >  1:3574590/1‑93 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  <  1:5577305/93‑1 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  >  2:4988751/1‑93 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  <  1:1515411/93‑1 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  <  1:1080804/93‑1 (MQ=255)
   cTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGcgcg                                                                                 <  2:6733215/93‑1 (MQ=255)
    tCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGa                                                                                >  2:7023477/1‑93 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGCAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            >  1:5878834/1‑93 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            >  2:5422934/1‑93 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATAGCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGAGGCAACGAGGCGCGCATACg                                                                           >  2:1364173/1‑93 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACg                                                                           >  1:4548213/1‑93 (MQ=255)
            gCGTGAATCATCAAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGa                                                                        <  1:3253253/93‑1 (MQ=255)
             cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                                       <  1:8072271/93‑1 (MQ=255)
               tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTc                                                                     >  2:914625/1‑93 (MQ=255)
                gAATCATCAAGGCGGCATTCGATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                                    >  1:3120569/1‑93 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                   <  1:3951458/93‑1 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                   >  2:6819509/1‑93 (MQ=255)
                     atcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAcc                                                               >  1:7144102/1‑93 (MQ=255)
                      tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                              >  1:6045887/1‑93 (MQ=255)
                        aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                            >  1:4729485/1‑93 (MQ=255)
                        aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                            <  2:4038977/93‑1 (MQ=255)
                         aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                           <  2:7158849/93‑1 (MQ=255)
                          ggcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAc                                                          <  1:6108631/93‑1 (MQ=255)
                           gcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                         >  1:6817139/1‑93 (MQ=255)
                            cggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGt                                                        <  1:6539762/93‑1 (MQ=255)
                             ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                       <  2:1758014/93‑1 (MQ=255)
                             ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                       <  1:2816974/93‑1 (MQ=255)
                               cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                     <  1:1001666/93‑1 (MQ=255)
                                atgcAACACCCACCGCCCCCTGTCCTACGGCCGGAAGAGAAATCTCGAGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  2:3498685/90‑1 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                   >  1:2001964/1‑93 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                   >  1:7439409/1‑93 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 >  1:1876543/1‑93 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  2:5812656/93‑1 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  2:6071435/93‑1 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  2:8091292/93‑1 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  1:958548/93‑1 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 >  1:8357304/1‑93 (MQ=255)
                                     aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCg                                               <  2:4709174/93‑1 (MQ=255)
                                        cccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                            >  2:2321425/1‑93 (MQ=255)
                                        cccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                            <  2:1194963/93‑1 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           <  2:854032/93‑1 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           <  1:8717730/93‑1 (MQ=255)
                                          cACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                          <  1:4209344/93‑1 (MQ=255)
                                          cACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                          <  2:5140625/93‑1 (MQ=255)
                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACg                                         >  1:5085457/1‑93 (MQ=255)
                                            ccGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                                        >  2:659194/1‑93 (MQ=255)
                                            ccGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                                        >  1:7260093/1‑93 (MQ=255)
                                             cGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                                       <  2:737696/93‑1 (MQ=255)
                                              gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                                      >  2:7745608/1‑93 (MQ=255)
                                              gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                                      >  1:5965393/1‑93 (MQ=255)
                                               cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCaa                                     >  2:8721915/1‑93 (MQ=255)
                                                aCCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAg                                    >  1:6235370/1‑93 (MQ=255)
                                                 ccTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACATGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                   >  1:2664577/1‑93 (MQ=255)
                                                  cTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                                  <  1:4337537/93‑1 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGTAAGAatg                                >  2:6212852/1‑93 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                <  1:7684728/93‑1 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                <  2:6410967/93‑1 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                >  2:8240380/1‑93 (MQ=255)
                                                     gTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                               >  1:5460417/1‑93 (MQ=255)
                                                     gTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                               <  2:4524813/93‑1 (MQ=255)
                                                      tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                              <  2:1097960/93‑1 (MQ=255)
                                                      tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                              <  2:314559/93‑1 (MQ=255)
                                                       ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                             <  2:2825745/93‑1 (MQ=255)
                                                        cTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg                            <  2:2471247/93‑1 (MQ=255)
                                                         tACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                           <  2:3860064/93‑1 (MQ=255)
                                                         tACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                           >  2:1866423/1‑93 (MQ=255)
                                                         tACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                           >  2:4817285/1‑93 (MQ=255)
                                                            cgccgGAAGAGAAATCTCGGGTGGCAACGTGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                        <  1:1396204/93‑1 (MQ=255)
                                                            cgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                        >  2:8367515/1‑93 (MQ=255)
                                                             gccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg                       <  2:8003445/93‑1 (MQ=255)
                                                              ccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt                      >  2:6887411/1‑93 (MQ=255)
                                                                ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                    <  1:3812115/93‑1 (MQ=255)
                                                                ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                    <  1:6598637/93‑1 (MQ=255)
                                                                ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                    >  1:4973927/1‑93 (MQ=255)
                                                                ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                    <  2:2996422/93‑1 (MQ=255)
                                                                ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAt                    <  2:1987074/93‑1 (MQ=255)
                                                                 gAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                   <  1:2215180/93‑1 (MQ=255)
                                                                 gAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                   <  1:1181578/93‑1 (MQ=255)
                                                                   agagAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                 >  2:5991116/1‑93 (MQ=255)
                                                                   agagAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                 >  1:84303/1‑93 (MQ=255)
                                                                    gagaAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc                >  2:4326823/1‑93 (MQ=255)
                                                                    aagaAATCTCGGGGTGGCAAGCGGCCCCGATTCCGTGCTCCCAACCTAAACGTTTTTGGTCGGGTTCCGCCAGAATGATGGCATCGTATTctc                <  2:4222011/92‑3 (MQ=255)
                                                                      gaAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCg              >  1:259629/1‑93 (MQ=255)
                                                                       aaaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGt             <  1:2976357/93‑1 (MQ=255)
                                                                         aTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAACCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTa           <  2:1616134/93‑1 (MQ=255)
                                                                         aTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTa           <  2:2678831/93‑1 (MQ=255)
                                                                         aTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTa           >  2:3382133/1‑93 (MQ=255)
                                                                           ctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATc         <  2:7552935/93‑1 (MQ=255)
                                                                           ctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATc         <  2:6330651/93‑1 (MQ=255)
                                                                            tcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATTGCATCGTATTCGCCGTTATcc        <  2:5539346/93‑1 (MQ=255)
                                                                            tcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATcc        >  1:538120/1‑93 (MQ=255)
                                                                             tggggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATcc        <  1:7972730/91‑1 (MQ=255)
                                                                             cGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa       <  2:2135120/93‑1 (MQ=255)
                                                                               gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg      >  1:635929/1‑93 (MQ=255)
                                                                                ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGt     >  1:5810347/1‑93 (MQ=255)
                                                                                 gTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGtt    <  2:3652450/93‑1 (MQ=255)
                                                                                   ggCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTg  >  1:1816667/1‑93 (MQ=255)
                                                                                        ||                                                                                      
CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG  >  NC_000913/3987735‑3987909

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: