BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	4,402,594	0	A	G	22.3%	249.6 / 25.1	104	intergenic (+53/+76)	ulaF/yjfY	L‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 family protein YjfY
Reads supporting (aligned to +/- strand): ref base A (62/18); new base G (0/23); total (62/42)
Fisher's exact test for biased strand distribution p-value = 3.78e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 132 total aligned reads displayed.

TTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGTTATACCT  >  NC_012967/4402523‑4402665

ttattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCa                                                                                 <  1:3675130/72‑1 (MQ=255)

ttattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCa                                                                                 <  2:2715541/72‑1 (MQ=255)

gtattcCGAGCAGAATTACAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGTTCTAATTCcaca                                                                                 <  1:2238627/71‑1 (MQ=255)

 tattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAg                                                                                >  2:1933949/1‑72 (MQ=255)

 tattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTACTCGGGAGATAACTGGTCTAATTCCGCAg                                                                                <  2:339384/72‑1 (MQ=255)

  attaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGTGAGATAACTGGTCTAATTCCGCAGc                                                                               >  1:1727065/1‑72 (MQ=255)

  attaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGc                                                                               >  1:2737084/1‑72 (MQ=255)

  attaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGc                                                                               >  2:2149230/1‑72 (MQ=255)

   ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGcc                                                                              <  1:4330115/72‑1 (MQ=255)

   ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGcc                                                                              <  2:6138465/72‑1 (MQ=255)

   ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGcc                                                                              <  1:5486615/72‑1 (MQ=255)

   ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTGCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGcc                                                                              <  1:934948/72‑1 (MQ=255)

   gtgcGGGGAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGcc                                                                              <  1:3067060/69‑1 (MQ=255)

    taCGGGCAGAAGTAGAAGACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCg                                                                             <  2:3856211/72‑1 (MQ=255)

     aCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑t                                                                           <  1:806796/72‑1 (MQ=255)

      cGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCgct                                                                           >  2:3897682/1‑70 (MQ=255)

      cGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑tt                                                                          >  2:6275347/1‑72 (MQ=255)

        ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑tttt                                                                        <  1:2242027/72‑1 (MQ=255)

        ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑tttt                                                                        <  2:5442793/72‑1 (MQ=255)

        ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑tttt                                                                        <  1:3206729/72‑1 (MQ=255)

         gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑ttttt                                                                       >  1:3953389/1‑72 (MQ=255)

         gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑ttttt                                                                       >  2:4123231/1‑72 (MQ=255)

         gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑ttttt                                                                       <  1:6013022/72‑1 (MQ=255)

         gCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑ttttt                                                                       >  1:6372262/1‑72 (MQ=255)

          cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTc                                                                      <  1:1862305/72‑1 (MQ=255)

          cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTc                                                                      <  1:2605362/72‑1 (MQ=255)

           aGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑a                                                                    <  1:2039081/72‑1 (MQ=255)

             aaGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑aaa                                                                  >  1:5023920/1‑72 (MQ=255)

              aGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑aaaa                                                                 <  1:4149105/72‑1 (MQ=255)

              aGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑aaaa                                                                 <  2:4397876/72‑1 (MQ=255)

              aGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑acaa                                                                 >  1:2533359/1‑72 (MQ=255)

              aGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑aaaa                                                                 <  1:5783640/72‑1 (MQ=255)

              aGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑aaaa                                                                 >  1:4384269/1‑72 (MQ=255)

               gTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑aaaaa                                                                >  1:4843766/1‑72 (MQ=255)

                tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTACTTCCGCAGCCG‑TGTTTC‑aaaaaa                                                               >  1:6241156/1‑72 (MQ=255)

                tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑aaaaaa                                                               >  2:1194548/1‑72 (MQ=255)

                  gAACCCGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTANTTCCGCGGCCG‑TTTTTC‑AAAAAAA‑g                                                            <  1:741820/72‑1 (MQ=255)

                  gAACACGCGCTGCGGAAATTTCCTTCNTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑g                                                            <  1:456827/72‑1 (MQ=255)

                  gAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑g                                                            >  2:6311806/1‑72 (MQ=255)

                   accacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑c                                                          <  1:946226/70‑1 (MQ=255)

                   aaCACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑c                                                          <  2:3204883/72‑1 (MQ=255)

                    acacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑cc                                                         <  1:6249101/72‑1 (MQ=255)

                     cacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTTCGCAACCGCTTTTTC‑AAAAAAA‑c‑cc                                                         >  2:5602034/1‑72 (MQ=255)

                      acGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑cccc                                                       <  1:4633675/72‑1 (MQ=255)

                        gcgcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCTGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑t                                                    <  1:3780710/72‑1 (MQ=255)

                        gcgcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTCAAAAAAAA‑G‑ccccc                                                      >  1:4089827/1‑72 (MQ=255)

                        gcgcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑t                                                    >  2:4976949/1‑72 (MQ=255)

                         ngcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑g                                                  <  2:2473961/71‑1 (MQ=255)

                            tGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑gag‑a                                              <  1:2834371/72‑1 (MQ=255)

                            tGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑gag‑a                                              <  1:5266847/72‑1 (MQ=255)

                             gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑aa                                             <  1:5662265/72‑1 (MQ=255)

                             gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑C‑GAG‑aa                                             >  2:6008488/1‑72 (MQ=255)

                             gCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CACCC‑T‑CAG‑aa                                             >  1:2124408/1‑72 (MQ=255)

                              cGGAAATTTCTTGCCTCGGGAGAAAACTGGTCTATTTCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑g                                           <  2:754579/72‑1 (MQ=255)

                               ggAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑gg                                          <  1:2442653/72‑1 (MQ=255)

                               ggAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑C‑GGG‑AA‑gg                                          >  1:5973503/1‑72 (MQ=255)

                                gAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑ggg                                         <  1:6459301/72‑1 (MQ=255)

                                gAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTCC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑ggg                                         >  1:4640645/1‑72 (MQ=255)

                                 aaaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑gggg                                        >  1:1816710/1‑72 (MQ=255)

                                 aaaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑C‑GTG‑AA‑aggg                                        >  2:4012937/1‑72 (MQ=255)

                                  attttCCTTCCTGGGGAGATAACTGGTCTAATCCCGCGGCCG‑TTTTTA‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑ggggg                                       <  1:4552125/70‑1 (MQ=255)

                                  aaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑ggggg                                       >  1:2606724/1‑72 (MQ=255)

                                   aTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGgc                                      >  1:4945114/1‑72 (MQ=255)

                                    tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑ggggggc                                     >  1:5929698/1‑70 (MQ=255)

                                     ttccttccTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑A‑GCGCC‑C‑GAG‑GA‑gngggggg                                    >  2:448276/1‑69 (MQ=255)

                                     tgcctgccgCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGTGgcaa                                    <  2:4499783/70‑1 (MQ=255)

                                      tccttcctCGGGATATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGgcaaa                                   >  1:3232691/1‑72 (MQ=255)

                                      tccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTCAAAAAAAA‑C‑CCCCC‑T‑CGG‑AA‑ggggggan                                    >  1:769317/1‑69 (MQ=255)

                                      tccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGggaaa                                   >  1:4330576/1‑72 (MQ=255)

                                         ttcctCGGGAGATAACTGGTCTAATTCCGCAGCAG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGggaaagca                                >  2:1554393/1‑72 (MQ=255)

                                           cgtCGGGAGATAACTGGTCTATTCCCGCGGCCG‑TTTTTA‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGgcaaagcaaa                              <  2:4586428/70‑1 (MQ=255)

                                            ctCGGGAGATAACTGGTCTAATTCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAAc                             <  2:6612693/72‑1 (MQ=255)

                                            ctCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAGAAA‑GGGGgcaaaccaac                              >  2:5980062/1‑71 (MQ=255)

                                            ctCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑C‑GAG‑AA‑AGGGGCGaaacaaac                             >  2:1962497/1‑72 (MQ=255)

                                              cGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCA‑T‑GAG‑AA‑GGGGGCAAAGCAAACTa                           <  1:2956055/72‑1 (MQ=255)

                                               gggAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCCAACTAt                          >  1:1653522/1‑72 (MQ=255)

                                               gggAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTAt                          >  2:4511743/1‑72 (MQ=255)

                                               gggAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑C‑GAG‑AA‑GGGGGGCAAACAAAatat                          >  2:1971970/1‑72 (MQ=255)

                                                 gagaTAACTGGTCTAATCCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATgg                        <  1:3175014/72‑1 (MQ=255)

                                                 aaaaTAACTGGTCTAATCCCGCGGCC‑‑TTTTTA‑AAAAAAA‑C‑CCCCCGT‑GAG‑AG‑GGGGGCAAAGCAAACTATgg                        <  2:5462913/69‑1 (MQ=255)

                                                   gaTAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑GCCCCCC‑C‑GGA‑AA‑GGGGGGAAAGCAAACTATgga                       >  2:4078725/1‑71 (MQ=255)

                                                   gaTAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GCGGGCAAAGCAAACTATGGCa                      <  2:4433473/72‑1 (MQ=255)

                                                   gaTAACTGCTCCAATTCCCCCGCCT‑TTTTAA‑AAAAAAACC‑CCCCC‑G‑GAG‑AG‑GGGGGCAAAGCAAACTATGGc                       <  1:5967789/72‑1 (MQ=255)

                                                    aTAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑GCCCCCC‑T‑GGG‑AA‑GGGGGGAAAGCAAACTATGGCa                      >  2:3809138/1‑72 (MQ=255)

                                                      aagtGGTTTAATCCCGCAGCCG‑TTTTTA‑AAAAAAA‑G‑CCCCC‑G‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATg                   <  2:1862592/69‑1 (MQ=255)

                                                      aaCTGGTCTAATCCCGCTGCCG‑TTTTTA‑AAAAAAA‑C‑CCCCC‑G‑GAG‑AG‑GGGGGCAAAGCAAACTATGGCAATg                   <  1:3119761/72‑1 (MQ=255)

                                                           ggttAATTCCGGCGCCC‑TTTTTA‑AAAAAAA‑A‑CCCCC‑TGGAG‑AG‑GGGGGCAAAGCAAACTATGGCAATGTTTc               <  2:4935750/69‑1 (MQ=255)

                                                            tCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑A‑CCCCC‑T‑GAG‑GA‑GGGGGGCAAGCAAACTATGGCACCGttttct             >  2:1409690/1‑69 (MQ=255)

                                                            tCTAACCCCGCGGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCtt             <  2:2854080/72‑1 (MQ=255)

                                                               aaTTCCGCAGCCG‑TTTTTC‑AAACAAA‑A‑CCCCC‑C‑GAG‑AAGGGGGGCAAACAAAATAATGGCAATGTTTTTTTg           >  2:508332/1‑72 (MQ=255)

                                                               aaTTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCttttg          >  2:1256761/1‑70 (MQ=255)

                                                                aTTCCGCAGCCG‑TTTTTCAAAAAAAA‑A‑GCCCC‑C‑CAA‑AA‑GGGGGGAAGGCGAAATTTTGTAATGTTTTTTTgt          >  1:972733/1‑71 (MQ=255)

                                                                 ttCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGtt        >  2:6544479/1‑72 (MQ=255)

                                                                  cccGCTGCCG‑TTTTTA‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AG‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGTta       <  2:6372259/71‑1 (MQ=255)

                                                                   ccGCGGCCG‑TTTTAA‑AAAAAAA‑C‑CCCCC‑G‑GAG‑AG‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGTtat      <  1:60281/72‑1 (MQ=255)

                                                                    cGCAGCCG‑TTTTTC‑AAAAAAA‑G‑GCCCC‑C‑CCG‑CA‑GGGGGCCAAGCAAAATACGGCAATGCTTCTTTTGtttaa     >  1:1493527/1‑69 (MQ=255)

                                                                     gCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCCAAGCAAACTCTTGCAATGCTTCTTTGGTTATAc    >  1:1669446/1‑72 (MQ=255)

                                                                     gCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGTTATAc    >  2:5057090/1‑72 (MQ=255)

                                                                     gCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGTTATAc    >  2:3278676/1‑72 (MQ=255)

                                                                       aGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGGAAAACAAAATATGGCAATGTTTCTTTGGTTATACCt  >  1:6092163/1‑72 (MQ=255)

TTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCG‑TTTTTC‑AAAAAAA‑G‑CCCCC‑T‑GAG‑AA‑GGGGGCAAAGCAAACTATGGCAATGTTTCTTTGGTTATACCT  >  NC_012967/4402523‑4402665

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