Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,170 +C intergenic (+172/+279) hemN → / ← glnG O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,1701.C99.5% 1311.9 / ‑1.4 432intergenic (+172/+279)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (234/196);  total (235/197)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 497 total aligned reads displayed.

CGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAG  >  NC_000913/4051074‑4051259
                                                                                                  |                                                                                         
  cTGGGCAACAAAATCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTTCTCCGCATCAGGTATCTGAGCATGATGCCCGATgcggcgct                                                                                          >  1:5586805/1‑94 (MQ=25)
   tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctg                                                                                         <  1:8142600/96‑4 (MQ=255)
    gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctgg                                                                                        <  2:1697761/96‑5 (MQ=255)
    gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctgg                                                                                        <  1:2124994/96‑5 (MQ=255)
      gCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                                                      <  2:7484651/96‑1 (MQ=25)
       caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCAGATGCGGCGCTGGTcg                                                                                     >  1:12806736/1‑96 (MQ=25)
        aacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                    >  2:7163069/1‑96 (MQ=25)
         accaaaCCTGTTGAAATTATTCCTTTTGTAGGCAGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                                                   <  2:4569965/93‑1 (MQ=18)
         acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                                                   >  2:1156133/1‑96 (MQ=25)
           aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                                                 >  2:4796935/1‑96 (MQ=25)
           aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                                                 >  2:8835907/1‑96 (MQ=25)
            aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                                                >  1:10952969/1‑96 (MQ=25)
             aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                               <  1:11261046/96‑1 (MQ=21)
             aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                               >  1:6199698/1‑96 (MQ=21)
             aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGACtt                                                                               <  2:10097625/96‑1 (MQ=21)
              aCCTGTTGGAATTATTGCTTTTGTAGGCCCGATAAAGCGGTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                                              >  2:12802600/1‑96 (MQ=21)
              aCCTGTTGAAATTATTGCTATTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                                              <  1:12671741/96‑1 (MQ=21)
                ccgTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                                            >  1:11787317/3‑96 (MQ=21)
                 tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                           <  1:8664789/96‑1 (MQ=21)
                 tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGGCTTATCa                                                                           >  2:9397403/1‑96 (MQ=21)
                  gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGGCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCACGTCTTATCAg                                                                          >  2:2437722/1‑96 (MQ=18)
                  gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGTGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                                          <  2:10169378/96‑1 (MQ=18)
                  gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                                          <  2:4384064/96‑1 (MQ=18)
                     gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                                       >  1:9942932/1‑96 (MQ=18)
                       aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                     >  2:9074167/1‑96 (MQ=18)
                       aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                     <  1:4397947/96‑1 (MQ=18)
                         ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                                                   <  1:13314449/96‑1 (MQ=18)
                         ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                                                   <  1:6682916/96‑1 (MQ=18)
                          tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                                  >  1:14292650/1‑96 (MQ=18)
                          tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                                  >  1:8431209/1‑96 (MQ=18)
                            ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                                                >  2:3538999/1‑96 (MQ=18)
                              gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCCCTGGTCGCGTCTTATCAGGCCTACGACttt                                                              >  2:14088940/1‑96 (MQ=18)
                               cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                                             >  1:2835339/1‑96 (MQ=18)
                                 tttGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGACTTATCAGGCCTACGACTTtctc                                                           >  2:600223/1‑96 (MQ=18)
                                  ttGTAGGTCAGATAAGGCGTTCACGCAGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                          <  2:515326/96‑1 (MQ=17)
                                  ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                          <  2:11672164/96‑1 (MQ=18)
                                  ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                          <  2:851961/96‑1 (MQ=18)
                                  ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                          <  2:13170458/96‑1 (MQ=18)
                                  tggTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                          >  1:10054210/3‑96 (MQ=18)
                                   tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                                         <  2:11274579/96‑1 (MQ=18)
                                   tGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                                                         >  1:7849562/1‑96 (MQ=18)
                                     tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                                       >  1:3780284/1‑96 (MQ=18)
                                     tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                                       <  2:5788783/96‑1 (MQ=18)
                                      aGGGCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                                      <  2:6435939/96‑1 (MQ=18)
                                      aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATTCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                                      >  2:10407579/1‑96 (MQ=18)
                                         ccGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                   <  2:11113895/96‑1 (MQ=18)
                                         ccGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                   >  2:9875307/1‑96 (MQ=18)
                                          ctgATAAGGCGTTTCCGCCGCATCAGGCATCTGAGCATGATGCGTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTCCTCCGCATCCgg                                                  <  1:12329607/94‑1 (MQ=18)
                                          cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                                                  >  1:3082599/1‑96 (MQ=18)
                                          cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                                                  >  1:10428584/1‑96 (MQ=18)
                                          cGGATAAGGCGTTTACGCCGCATCAGGCAGCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCATGCCTACGACTTTCTCCGCATCCgg                                                  >  1:12053958/1‑96 (MQ=18)
                                           ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                                 >  2:3736126/1‑96 (MQ=18)
                                           ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                                 <  2:4947944/96‑1 (MQ=18)
                                             aTAAGGCGTTTCCGCCGCATCAGGCATCTGAGCATGATCCCTGATGCGGCGCTGGTCGAGTCTTACCAGGCCTACGACTTCCTCCGCATCCGGCAt                                               <  2:6245505/96‑1 (MQ=255)
                                             aTAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                                               >  2:11949203/1‑96 (MQ=21)
                                                 ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                           <  2:8255335/96‑1 (MQ=21)
                                                 ggCGTTTACGCCGCATCAGGCATCGGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                           <  2:13796885/96‑1 (MQ=21)
                                                  gCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                          >  2:4163708/1‑96 (MQ=21)
                                                  gCGGTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                          >  1:11554322/1‑96 (MQ=21)
                                                   cGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                                         <  1:5103545/96‑1 (MQ=21)
                                                    ggttACGCCGCATCAGGCCTTCGATCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTCTGACTTACTCCGTATCCGGCATAaacacc                                        <  1:6801266/94‑3 (MQ=18)
                                                    gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                                        <  2:10567132/96‑1 (MQ=21)
                                                    gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                                        >  1:467221/1‑96 (MQ=25)
                                                    gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                                        >  1:9793055/1‑96 (MQ=25)
                                                    gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                                        >  1:1563265/1‑96 (MQ=21)
                                                     tttACGTCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                                       >  1:1773399/1‑96 (MQ=25)
                                                     tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg                                       >  1:8975102/1‑96 (MQ=25)
                                                      ttccgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa                                      <  2:12091138/93‑1 (MQ=25)
                                                       tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg                                     >  2:13541516/1‑96 (MQ=25)
                                                        aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                    <  2:11748707/96‑1 (MQ=25)
                                                        aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                    >  2:2973071/1‑96 (MQ=25)
                                                           ccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                                 >  1:5501287/1‑96 (MQ=25)
                                                            cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                                <  1:2505444/96‑1 (MQ=37)
                                                            cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                                >  2:7844199/1‑96 (MQ=37)
                                                             gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACATCGAGCTTCAg                               >  2:10716575/1‑96 (MQ=25)
                                                              cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                              <  1:8329822/96‑1 (MQ=37)
                                                               aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                             <  2:1261696/96‑1 (MQ=37)
                                                               aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                             >  1:13584704/1‑96 (MQ=37)
                                                                tCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGACCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                            <  2:12900562/96‑1 (MQ=38)
                                                                  caggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaa                           >  1:6399911/5‑96 (MQ=38)
                                                                   aggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTaaa                          <  1:3410852/93‑1 (MQ=255)
                                                                    ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                         <  2:10879538/94‑1 (MQ=255)
                                                                      cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAActg                       <  1:10545781/96‑3 (MQ=255)
                                                                      cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCACTTAAACAg                       >  2:6892525/1‑96 (MQ=37)
                                                                           gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa                  >  2:7761875/1‑96 (MQ=255)
                                                                            aGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt                 >  2:8114994/1‑96 (MQ=255)
                                                                             gaatgatgCCTGATGCGGGGCTGGTCGCGTCTTATCAGGTCTACGACTTTCTCCGCATCCGGCATAAACAACTAGCTTCAGCTAAACAGCCCAATc                >  1:10149402/3‑96 (MQ=38)
                                                                               atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAACCCcatcat              >  1:14088509/1‑96 (MQ=38)
                                                                                tgatgCCTGCTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAAGAGCTTCAGCTAAACAGCCCAATCAtt             <  1:3816914/96‑1 (MQ=38)
                                                                                tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCTGCTAAACAGCCCAATCAtt             <  2:12356870/96‑1 (MQ=38)
                                                                                tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAtt             >  2:11379189/1‑96 (MQ=255)
                                                                                  atgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc           >  1:4091096/1‑96 (MQ=255)
                                                                                   tgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg          <  2:5640148/96‑1 (MQ=255)
                                                                                     ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGc        <  1:7217052/96‑1 (MQ=255)
                                                                                     ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGc        >  2:3407176/1‑96 (MQ=255)
                                                                                     ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGc        <  1:1305780/96‑1 (MQ=255)
                                                                                       tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac      <  1:9552471/96‑1 (MQ=255)
                                                                                       tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac      <  2:2102586/96‑1 (MQ=255)
                                                                                       tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac      >  1:13874040/1‑96 (MQ=255)
                                                                                          tGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcacaca   <  2:14292650/96‑1 (MQ=255)
                                                                                                  |                                                                                         
CGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAG  >  NC_000913/4051074‑4051259

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