Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,144 Δ1 bp intergenic (+146/+305) 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,1430C.99.8% 1549.6 / ‑3.5 454intergenic (+145/+306)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (231/222);  total (231/223)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 488 total aligned reads displayed.

TGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTA  >  NC_000913/4051048‑4051233
                                                                                               |                                                                                           
     tGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                       >  1:13787877/1‑96 (MQ=255)
      gagaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGa                                                                                      >  1:515326/1‑96 (MQ=255)
       agaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAg                                                                                     >  2:6138455/1‑96 (MQ=255)
         aaCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCa                                                                                   <  2:7724013/96‑1 (MQ=255)
          aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCat                                                                                  <  1:7078258/96‑1 (MQ=255)
            aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGTCGCATCAGGCATCTGAGCatga                                                                                >  1:14240181/1‑96 (MQ=255)
            aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                                                <  2:5407721/96‑1 (MQ=255)
              gATTGCCTGATGCTCTTGGCAACAATACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGTCGTTTACGCCGCATCAGGCATCTGAGCatgatg                                                                              <  1:8386510/96‑1 (MQ=25)
                 tGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                           >  2:3787020/1‑96 (MQ=255)
                  gCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTg                                                                          <  2:5076034/96‑1 (MQ=255)
                   ccTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                         <  2:1065724/96‑1 (MQ=255)
                     tGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATg                                                                       >  1:3269301/1‑96 (MQ=255)
                     tGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATg                                                                       >  1:6526309/1‑96 (MQ=255)
                         gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                                                   >  2:13053365/1‑96 (MQ=255)
                         gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                                                   <  2:2201653/96‑1 (MQ=255)
                         gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                                                   >  1:5416927/1‑96 (MQ=255)
                         gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                                                   <  1:1557230/96‑1 (MQ=255)
                         gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCAGCTGAGCATGATGCCTGATgcggc                                                                   >  1:12525548/1‑96 (MQ=255)
                         gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGGTTACGACGCATCTGGCATCGGAGCATGATGCCTGCTgcggc                                                                   >  1:41548/1‑96 (MQ=25)
                          cgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                                  <  1:8191304/96‑1 (MQ=255)
                           ggtGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgc                                                                 >  1:12082463/3‑95 (MQ=255)
                            cTGGGCAACAAAATCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTTCTCCGCATCAGGTATCTGAGCATGATGCCCGATgcggcgct                                                                >  1:5586805/1‑94 (MQ=25)
                            cTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGtt                                                               <  1:9195822/96‑2 (MQ=255)
                             tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctg                                                               >  2:4763362/1‑93 (MQ=255)
                             tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctg                                                               <  2:8879598/96‑4 (MQ=255)
                                 caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                           <  1:6201055/96‑1 (MQ=25)
                                  aacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                          <  2:8501277/96‑1 (MQ=25)
                                      aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                      >  2:8172776/1‑96 (MQ=21)
                                       aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                     <  2:11030828/96‑1 (MQ=21)
                                        aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                    <  2:56829/96‑1 (MQ=21)
                                          ccgTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                  >  1:11787317/3‑96 (MQ=21)
                                          cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                  <  2:14274272/96‑1 (MQ=21)
                                          cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGACCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                  <  2:5909081/96‑1 (MQ=21)
                                           tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                 >  2:7255010/1‑96 (MQ=21)
                                           tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGGCTTATCa                                                 >  2:9397403/1‑96 (MQ=21)
                                            gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTCTTAGCAg                                                >  1:9913918/1‑96 (MQ=18)
                                                aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                            <  2:176272/96‑1 (MQ=18)
                                                 aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                           <  2:12229348/96‑1 (MQ=18)
                                                 aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                           <  1:1011988/96‑1 (MQ=18)
                                                  attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                          >  2:9846636/1‑96 (MQ=18)
                                                   ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                         <  1:13314449/96‑1 (MQ=18)
                                                    tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                        >  2:4371246/1‑96 (MQ=18)
                                                       tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                     <  2:10078535/96‑1 (MQ=18)
                                                        gCTTTTGTAGGTCGTATTAGTCGTTTACGCCGAATCAGGCATCTGAGCATTATGCCTGATGCGGCGCTGGTCGTGTCTTATCAGGTCTCCGACttt                                    <  1:6787202/96‑1 (MQ=255)
                                                        gCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                    >  1:10862618/1‑96 (MQ=18)
                                                         cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtc                                   >  1:2835339/1‑96 (MQ=18)
                                                         cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGCCTTATCAGGCCTACGCCTTtc                                   >  2:3575929/1‑96 (MQ=17)
                                                            ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                <  2:4509319/96‑1 (MQ=18)
                                                            ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                <  1:6041286/96‑1 (MQ=18)
                                                            ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTCTTATCAGGCCTACGACTTTCTcc                                >  2:68910/1‑96 (MQ=18)
                                                            tggTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                >  1:10054210/3‑96 (MQ=18)
                                                              gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                              >  2:148865/1‑96 (MQ=18)
                                                               tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCGTCTGAGCATGATGCCTGATGCGGCGCTGGGCACGTCTTATCAGGCCTACGACTTTCTCCGCa                             >  2:14088328/1‑96 (MQ=18)
                                                               tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                             >  2:11349996/1‑96 (MQ=18)
                                                                aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATTCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                            >  2:10407579/1‑96 (MQ=18)
                                                                aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCGGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                            >  2:10298649/1‑96 (MQ=18)
                                                                 ggCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                           >  1:13982810/1‑96 (MQ=18)
                                                                  gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                          >  1:2391132/1‑96 (MQ=18)
                                                                   ccGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                         <  2:11196084/96‑1 (MQ=18)
                                                                    ctgATAAGGCGTTTCCGCCGCATCAGGCATCTGAGCATGATGCGTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTCCTCCGCATCCgg                        <  1:12329607/94‑1 (MQ=18)
                                                                    cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                        <  2:115799/96‑1 (MQ=18)
                                                                    cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                        >  1:3082599/1‑96 (MQ=18)
                                                                     ggATAAGGCTTTTACGCCGCATCGGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                       <  1:9585795/96‑1 (MQ=18)
                                                                     ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                       <  2:4947944/96‑1 (MQ=18)
                                                                       aTAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                     <  2:12798188/96‑1 (MQ=18)
                                                                       aTAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                     >  2:11949203/1‑96 (MQ=21)
                                                                        tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                    >  2:2791594/1‑96 (MQ=21)
                                                                         aaGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                   >  2:55218/1‑96 (MQ=21)
                                                                         aaGGCGTTGACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                   >  1:12992447/1‑96 (MQ=21)
                                                                           ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                 >  1:7572390/1‑96 (MQ=25)
                                                                           ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                 >  2:11697717/1‑96 (MQ=25)
                                                                           ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCAGCATCCGGCATAaac                 >  1:2118306/1‑96 (MQ=21)
                                                                           ggCGTTTACGCCGCATCAGGCATCGGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                 <  2:13796885/96‑1 (MQ=21)
                                                                             cGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa               >  1:346137/1‑96 (MQ=25)
                                                                             cGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa               >  1:4440933/1‑96 (MQ=21)
                                                                             agTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa               >  1:11363059/2‑96 (MQ=25)
                                                                              gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac              >  2:9149307/1‑96 (MQ=25)
                                                                              gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac              <  2:10567132/96‑1 (MQ=21)
                                                                              gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac              >  1:467221/1‑96 (MQ=25)
                                                                              gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac              >  1:8847493/1‑96 (MQ=21)
                                                                               tttACGTCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg             >  1:1773399/1‑96 (MQ=25)
                                                                               tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg             >  2:6391065/1‑96 (MQ=25)
                                                                               tatACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACg             >  1:10507051/3‑96 (MQ=25)
                                                                                ttccgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa            <  2:13361303/93‑1 (MQ=25)
                                                                                  aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc          <  2:10147899/96‑1 (MQ=25)
                                                                                  aCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc          <  2:3239646/96‑1 (MQ=25)
                                                                                   cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCt         >  1:3791105/1‑96 (MQ=25)
                                                                                    gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt        >  2:6462864/1‑96 (MQ=37)
                                                                                    gacgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtg        >  2:8736157/3‑95 (MQ=25)
                                                                                     ccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc       >  2:3970905/1‑96 (MQ=37)
                                                                                      cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa      >  1:8645282/1‑96 (MQ=37)
                                                                                      cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa      >  2:7844199/1‑96 (MQ=37)
                                                                                      cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa      <  1:11179736/96‑1 (MQ=37)
                                                                                        cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATTAACAACGAGCTTCAGc    >  1:10333492/1‑96 (MQ=37)
                                                                                        cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATCAACAACGAGCTTCata    >  1:12783683/1‑94 (MQ=25)
                                                                                        cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc    >  2:1178824/1‑96 (MQ=37)
                                                                                        cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc    <  1:8329822/96‑1 (MQ=37)
                                                                                        cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc    >  1:13517443/1‑96 (MQ=37)
                                                                                         aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt   >  1:8072235/1‑96 (MQ=37)
                                                                                         aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt   <  2:1261696/96‑1 (MQ=37)
                                                                                               |                                                                                           
TGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTA  >  NC_000913/4051048‑4051233

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