BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,649	A→G	V201A (GTG→GCG)	yagY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,649	0	A	G	97.0%	513.4 / 3.2	168	V201A (GTG→GCG)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/4); new base G (83/80); total (84/84)
Fisher's exact test for biased strand distribution p-value = 3.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 178 total aligned reads displayed.

TCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAG  >  NC_000913/308559‑308738

 cGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCgcgc                                                                                        <  1:2169653/93‑1 (MQ=255)

 cGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCgcgc                                                                                        >  1:1351544/1‑93 (MQ=255)

  gTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCcgcgcg                                                                                       >  1:2005296/1‑93 (MQ=255)

   ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAATCCGTCCTTTGTTATCCGCAGTGTCCGCGCgg                                                                                      >  1:1711261/1‑93 (MQ=255)

   ttATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCgg                                                                                      <  1:4517643/93‑1 (MQ=255)

    tATCCATCTGCACAGCTATCTATTTCCCGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCggg                                                                                     >  1:584510/1‑93 (MQ=255)

     aTCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGt                                                                                    >  2:4129005/1‑93 (MQ=255)

      tCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTa                                                                                   >  1:3415380/1‑93 (MQ=255)

       ccATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaa                                                                                  <  2:239827/93‑1 (MQ=255)

       ccATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTaa                                                                                  <  2:2010423/93‑1 (MQ=255)

        cATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTaaa                                                                                 >  1:2752215/1‑93 (MQ=255)

          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               <  2:3517780/93‑1 (MQ=255)

          tCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               <  2:3517779/93‑1 (MQ=255)

          acTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAac                                                                               <  2:3920152/92‑1 (MQ=255)

            tGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacga                                                                             <  1:1280767/93‑1 (MQ=255)

             gCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAacgac                                                                            <  2:1934179/93‑1 (MQ=255)

               acaGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAcgc                                                                          <  1:907019/93‑1 (MQ=255)

               acaGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAcgc                                                                          <  2:4629703/93‑1 (MQ=255)

                caGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGGGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAcgcg                                                                         >  1:1128422/1‑93 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  1:844348/93‑1 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       >  1:2036447/1‑93 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  2:4587415/93‑1 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       >  1:3760076/1‑93 (MQ=255)

                  gCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGAc                                                                       <  2:4080805/93‑1 (MQ=255)

                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCACCCCGTACTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     <  2:4456153/93‑1 (MQ=255)

                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     <  1:1011112/93‑1 (MQ=255)

                    tatctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACtt                                                                     >  2:1417229/1‑93 (MQ=255)

                     atctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTg                                                                    >  1:602758/1‑93 (MQ=255)

                     atctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTg                                                                    >  2:3372541/1‑93 (MQ=255)

                       ctatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGcc                                                                  <  2:1898547/93‑1 (MQ=255)

                        tatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGccc                                                                 <  2:2581275/93‑1 (MQ=255)

                        tatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGccc                                                                 <  2:3119832/93‑1 (MQ=255)

                        tatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGccc                                                                 >  1:2731483/1‑93 (MQ=255)

                        tatTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTAccc                                                                 >  2:1526913/1‑93 (MQ=255)

                         atTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCg                                                                >  1:3058582/1‑93 (MQ=255)

                         atTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCg                                                                >  1:4644723/1‑93 (MQ=255)

                              aCGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATc                                                           >  2:1387815/1‑93 (MQ=255)

                               cGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTACCCGGCATCa                                                          >  1:4720614/1‑93 (MQ=255)

                                gggAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAg                                                         >  1:2412301/1‑93 (MQ=255)

                                 ggAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                                                        <  1:3019461/93‑1 (MQ=255)

                                   aaTGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCCTCCggaa                                                      >  2:3831715/1‑91 (MQ=255)

                                   aaTGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgta                                                      >  1:4434929/1‑93 (MQ=255)

                                    aTGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtag                                                     >  2:145950/1‑93 (MQ=255)

                                    aTGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtag                                                     <  1:2250566/93‑1 (MQ=255)

                                     tGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGAGGGNACTTGCCCGGCATCAGgtagn                                                    >  1:4085088/1‑92 (MQ=255)

                                     tGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagt                                                    <  2:2696043/93‑1 (MQ=255)

                                     tGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagt                                                    <  1:2792245/93‑1 (MQ=255)

                                     tGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagt                                                    >  1:3104090/1‑93 (MQ=255)

                                     tGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagt                                                    >  1:1180244/1‑93 (MQ=255)

                                      gAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGgtagta                                                   <  1:4234670/93‑1 (MQ=255)

                                          ttATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtc                                               >  1:2623046/1‑93 (MQ=255)

                                            aTCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATtctc                                             <  1:3507732/93‑1 (MQ=255)

                                             tCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTACGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTctct                                            >  2:916801/1‑93 (MQ=255)

                                               aCCCTGCCAAAGTGCAGCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCttt                                          <  2:2191845/93‑1 (MQ=255)

                                               aCCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAATGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCttt                                          >  2:424539/1‑93 (MQ=255)

                                                  cTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                                       >  2:4532068/1‑93 (MQ=255)

                                                  cTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                                       <  1:2568734/93‑1 (MQ=255)

                                                  cTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTaca                                       >  2:1818552/1‑93 (MQ=255)

                                                   tGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTGCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac                                      >  2:2340191/1‑93 (MQ=255)

                                                   tGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac                                      <  2:4546079/93‑1 (MQ=255)

                                                      cAAAGTGCAGCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTcc                                   <  2:119979/93‑1 (MQ=255)

                                                      cAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTcc                                   >  1:4596426/1‑93 (MQ=255)

                                                      cAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTcc                                   >  1:4749014/1‑93 (MQ=255)

                                                       aaaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCt                                  <  1:2203598/93‑1 (MQ=255)

                                                        aaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt                                 >  2:2035253/1‑93 (MQ=255)

                                                        aaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt                                 >  1:553021/1‑93 (MQ=255)

                                                        aaGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCtt                                 <  2:2839026/93‑1 (MQ=255)

                                                            gCAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCg                             >  1:368382/1‑93 (MQ=255)

                                                             cAACCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGt                            <  1:445150/93‑1 (MQ=255)

                                                              caCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           <  1:410823/92‑1 (MQ=255)

                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           <  2:1325841/93‑1 (MQ=255)

                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           >  2:2885652/1‑93 (MQ=255)

                                                              aaCCCGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGtt                           >  2:2721582/1‑93 (MQ=255)

                                                                cccGTCCTTTGTTATCCGCAGTGTCCGCGCGGTTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTgg                         <  1:1054034/93‑1 (MQ=255)

                                                                cccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTgg                         <  2:1825796/93‑1 (MQ=255)

                                                                 ccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggc                        <  2:1675833/93‑1 (MQ=255)

                                                                 ccGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggc                        >  2:3968445/1‑93 (MQ=255)

                                                                  cGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcg                       >  2:3484885/1‑93 (MQ=255)

                                                                  cGTCCTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTACCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcg                       <  1:3196994/93‑1 (MQ=255)

                                                                    tCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggcggc                     <  2:1586391/93‑1 (MQ=255)

                                                                      cTTTGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCtt                   <  2:2139413/93‑1 (MQ=255)

                                                                       tttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCttt                  <  2:2926873/93‑1 (MQ=255)

                                                                       tttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCttt                  <  2:105198/93‑1 (MQ=255)

                                                                        ttGTTATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTc                 <  2:619359/93‑1 (MQ=255)

                                                                           gtATCCGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAgg              <  1:1640138/92‑1 (MQ=255)

                                                                               ccGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAg          <  2:1875448/93‑1 (MQ=255)

                                                                                ctcAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAgg         >  1:198331/3‑93 (MQ=255)

                                                                                cGCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAgg         <  2:3696353/93‑1 (MQ=255)

                                                                                 gCAGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGt        >  1:141768/1‑93 (MQ=255)

                                                                                   aGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTcc      <  2:3630373/93‑1 (MQ=255)

                                                                                   aGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTcc      <  1:2915037/93‑1 (MQ=255)

                                                                                   aGTGTCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTcc      <  1:3845214/93‑1 (MQ=255)

                                                                                   aGTGTCCGCGCGGGTAAAACGACGCGAATTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTcc      <  2:4260348/93‑1 (MQ=255)

                                                                                    gtgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCg     >  1:3641508/1‑93 (MQ=255)

                                                                                     tgtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGt    >  2:3455786/1‑93 (MQ=255)

                                                                                     tgtCCGCGCGGGTAAAACGACGCGACTTACCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGt    <  2:1209181/93‑1 (MQ=255)

                                                                                      gtCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTa   >  2:1127183/1‑93 (MQ=255)

                                                                                       tCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAg  <  1:826301/93‑1 (MQ=255)

                                                                                       tCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAg  <  2:2515403/93‑1 (MQ=255)

                                                                                       tCCGCGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAg  <  2:4721297/93‑1 (MQ=255)

TCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAG  >  NC_000913/308559‑308738

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: `‑`