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% 1968.0 / NA 503G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (251/252);  total (251/252)
*NC_0009133,987,8230CG99.8% 1909.7 / NA 498G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (248/249);  minor base T (0/1);  total (248/250)
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 527 total aligned reads displayed.

AAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTC  >  NC_000913/3987732‑3987912
                                                                                          ||                                                                                         
aaGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgc                                                                                       >  2:1190616/1‑96 (MQ=255)
 agcagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACgcggcg                                                                                      >  2:343172/1‑96 (MQ=255)
  gcagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                     >  1:393947/1‑96 (MQ=255)
  gcagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                     <  1:45944/96‑1 (MQ=255)
    agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGa                                                                                   >  2:11488085/1‑96 (MQ=255)
      cTCGCGAGTGCGTGAATCATCAAGGCGGCATTCGATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAAt                                                                                 >  2:14082447/1‑96 (MQ=255)
      cTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAAt                                                                                 <  1:1821112/96‑1 (MQ=255)
        cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                               <  1:8216648/96‑1 (MQ=255)
        cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                               >  1:2784643/1‑96 (MQ=255)
        cgcgAGTGCGTGAATCATCAAGGAGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                               <  1:12477699/96‑1 (MQ=255)
         gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACg                                                                              <  2:3987021/96‑1 (MQ=255)
          cgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGCAATCTCGGGTGGCAACGCGGCGCGAATACGt                                                                             >  2:11214138/1‑96 (MQ=255)
           gAGTGCGTGAATCATCAAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTg                                                                            <  2:6037032/96‑1 (MQ=255)
           gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTg                                                                            <  1:12592400/96‑1 (MQ=255)
           gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTg                                                                            >  2:10998447/1‑96 (MQ=255)
           gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTg                                                                            >  2:6982905/1‑96 (MQ=255)
             gTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGCATACGTGAc                                                                          >  1:14358116/1‑96 (MQ=255)
              tGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACt                                                                         >  1:1189964/1‑96 (MQ=255)
              tGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACt                                                                         <  1:5098837/96‑1 (MQ=255)
               gCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCACCGCGGCGCGAATACGTGACTc                                                                        <  1:3620510/96‑1 (MQ=255)
               gCGTGAATCATCAAGGCGGCATTCAATACCCACCCCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGAGGCAACGCGGCGCGAATACGTGACTc                                                                        >  1:9289173/1‑96 (MQ=255)
                cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                                       <  2:574423/96‑1 (MQ=255)
                 gTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                      >  1:8521053/1‑96 (MQ=255)
                 gTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                      >  1:12201786/1‑96 (MQ=255)
                  tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcccg                                                                     <  1:11101043/96‑3 (MQ=255)
                  tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                                     >  2:951045/1‑96 (MQ=255)
                  tGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAg                                                                     <  2:2300143/96‑1 (MQ=255)
                   gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGa                                                                    <  1:4321638/96‑1 (MQ=255)
                      tcatcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                                 >  2:2554166/1‑96 (MQ=255)
                       catcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCTCTAATACGTGACTGCAGACCTa                                                                >  1:5354832/1‑96 (MQ=255)
                       catcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTa                                                                <  1:641897/96‑1 (MQ=255)
                        ttcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                               >  1:13819310/2‑96 (MQ=255)
                        atcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                               >  2:13307048/1‑96 (MQ=255)
                        atcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                               >  2:14358566/1‑96 (MQ=255)
                         tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                              >  2:10013796/1‑96 (MQ=255)
                         tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                              >  2:4356788/1‑96 (MQ=255)
                          caAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGCCCTAAAc                                                             >  2:6791304/1‑96 (MQ=255)
                          caAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaacc                                                             <  1:10263217/96‑3 (MQ=255)
                          caAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAc                                                             >  2:379169/1‑96 (MQ=255)
                           aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                            >  1:12978330/1‑96 (MQ=255)
                           aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                            >  2:9341002/1‑96 (MQ=255)
                           aaGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACg                                                            >  2:3323823/1‑96 (MQ=255)
                              gcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGttt                                                         >  2:5107233/1‑96 (MQ=255)
                               cgccATTCAATAGCCACCGCACCTTGTCCTACCTCCGTAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACTTGACTCCAAACCTAAACGtttt                                                        <  1:13450364/93‑1 (MQ=255)
                                ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                       >  2:7570869/1‑96 (MQ=255)
                                 gcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                      <  1:7975946/96‑1 (MQ=255)
                                  cATTGAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                     >  2:6383257/1‑96 (MQ=255)
                                  cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTact                                                     <  1:10621428/96‑3 (MQ=255)
                                  cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                     <  1:7017784/96‑1 (MQ=255)
                                  cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                     >  2:10346835/1‑96 (MQ=255)
                                   aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGACTACGTGACTCCAGACCTAAACGTTTTAGTc                                                    >  1:6768212/1‑96 (MQ=255)
                                   aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                    >  1:12774004/1‑96 (MQ=255)
                                   aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                    >  2:12524257/1‑96 (MQ=255)
                                      cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCCACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                                 <  2:13959270/96‑1 (MQ=255)
                                      cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcccg                                                 <  1:11668813/96‑3 (MQ=255)
                                      cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                                 >  2:293344/1‑96 (MQ=255)
                                      cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCgg                                                 >  1:5546511/1‑96 (MQ=255)
                                       aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                                <  1:3463812/96‑1 (MQ=255)
                                       aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGc                                                >  1:2664066/1‑96 (MQ=255)
                                        aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                               <  1:1321510/96‑1 (MQ=255)
                                        aTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                               <  1:5734359/96‑1 (MQ=255)
                                          aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                             <  1:3784312/96‑1 (MQ=255)
                                          aCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                             >  2:10699931/1‑96 (MQ=255)
                                            ccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                                           <  2:1797133/96‑1 (MQ=255)
                                             cACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGc                                          <  1:5852806/96‑1 (MQ=255)
                                                 gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                      >  2:13012468/1‑96 (MQ=255)
                                                 gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                      >  2:12739246/1‑96 (MQ=255)
                                                 gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                      <  1:8795339/96‑1 (MQ=255)
                                                 gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                      >  1:10011494/1‑96 (MQ=255)
                                                  cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                                     <  2:863435/96‑1 (MQ=255)
                                                   aCCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                    >  2:7273126/1‑96 (MQ=255)
                                                     cTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                                  <  2:3532747/96‑1 (MQ=255)
                                                       tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                                <  1:7328886/96‑1 (MQ=255)
                                                         tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                              <  1:11194648/96‑1 (MQ=255)
                                                         tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                              >  2:12250652/1‑96 (MQ=255)
                                                          ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGATTGATGGCa                             >  2:8338011/1‑96 (MQ=255)
                                                          ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCa                             >  2:9135124/1‑96 (MQ=255)
                                                           cTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCAt                            <  1:9490209/96‑1 (MQ=255)
                                                            tACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAACCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                           <  1:14082447/96‑1 (MQ=255)
                                                             aCCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg                          <  2:4790696/96‑1 (MQ=255)
                                                               cgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa                        >  1:597435/1‑96 (MQ=255)
                                                                 ccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                      >  2:886867/1‑96 (MQ=255)
                                                                  cgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTc                     >  2:5353984/1‑96 (MQ=255)
                                                                  cgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTc                     <  2:3094839/96‑1 (MQ=255)
                                                                   ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                    <  2:13763897/96‑1 (MQ=255)
                                                                   ggAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                    <  2:10145695/96‑1 (MQ=255)
                                                                    gAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc                   >  1:8359804/1‑96 (MQ=255)
                                                                       gagaAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAGGATGGCATCGGATTCGCCgg                >  1:4203444/1‑95 (MQ=255)
                                                                          aaaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAt             <  2:1352271/96‑1 (MQ=255)
                                                                           aaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATc            <  1:3011873/96‑1 (MQ=255)
                                                                           aaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATc            <  2:7802144/96‑1 (MQ=255)
                                                                            aTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATcc           <  2:8254671/96‑1 (MQ=255)
                                                                             tctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa          >  1:2269110/1‑96 (MQ=255)
                                                                             tctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa          <  1:10356637/96‑1 (MQ=255)
                                                                              cgcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg         >  1:10785899/3‑96 (MQ=255)
                                                                               tcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGt        <  1:4485994/96‑1 (MQ=255)
                                                                                 gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCCCCGTTATCCAGttt      <  1:1546111/96‑1 (MQ=255)
                                                                                  ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTg     >  1:2147240/1‑96 (MQ=255)
                                                                                  ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTg     >  1:9581161/1‑96 (MQ=255)
                                                                                     ggCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTc  >  2:1420733/1‑96 (MQ=255)
                                                                                          ||                                                                                         
AAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTC  >  NC_000913/3987732‑3987912

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