BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	198,418	A→G	D164G (GAC→GGC)	yaeT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	198,418	0	A	G	99.7%	1083.7 / NA	344	D164G (GAC→GGC)	yaeT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (159/184); minor base T (0/1); total (159/185)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 350 total aligned reads displayed.

AGAAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGACCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGA  >  NC_000913/198328‑198507

agaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGcc                                                                                         <  1:6117931/93‑1 (MQ=255)

 gaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCt                                                                                        >  1:4945880/1‑93 (MQ=255)

 gaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCt                                                                                        >  2:1365620/1‑93 (MQ=255)

 gaAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCt                                                                                        >  2:3772979/1‑93 (MQ=255)

  aaaGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTa                                                                                       <  1:3977808/93‑1 (MQ=255)

    aGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTaaa                                                                                     >  2:2320978/1‑93 (MQ=255)

    aGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTaaa                                                                                     >  1:6294649/1‑93 (MQ=255)

     ggTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTaaaa                                                                                    >  1:3628861/1‑93 (MQ=255)

      gTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAAc                                                                                   >  2:7486089/1‑93 (MQ=255)

      gTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAAc                                                                                   <  2:3967425/93‑1 (MQ=255)

           gAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtg                                                                              >  1:1023225/1‑93 (MQ=255)

            aaGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGgtgt                                                                             >  2:2016101/1‑93 (MQ=255)

             aGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGtt                                                                            >  2:5420332/1‑93 (MQ=255)

              gACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTc                                                                           <  1:689522/93‑1 (MQ=255)

               aCTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTcc                                                                          <  1:5015907/93‑1 (MQ=255)

                cTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCa                                                                         <  1:7486648/93‑1 (MQ=255)

                 ttCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAg                                                                        >  2:2905908/1‑93 (MQ=255)

                 ttCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAg                                                                        >  2:4916618/1‑93 (MQ=255)

                   ttactacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGa                                                                      <  1:434968/92‑1 (MQ=255)

                     actacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAg                                                                    <  1:5847982/93‑1 (MQ=255)

                     actacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAg                                                                    <  1:7920095/93‑1 (MQ=255)

                       tacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAggg                                                                  >  1:1987546/1‑92 (MQ=255)

                       tacAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgt                                                                  <  2:1592427/93‑1 (MQ=255)

                        acAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtg                                                                 >  2:6701192/1‑93 (MQ=255)

                         cAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtgt                                                                <  1:974987/93‑1 (MQ=255)

                          aGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtgtg                                                               <  2:2118729/93‑1 (MQ=255)

                           gCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGgtgtgt                                                              >  2:5282300/1‑93 (MQ=255)

                            cGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTc                                                             <  1:5219061/93‑1 (MQ=255)

                            cGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTc                                                             <  2:2564075/93‑1 (MQ=255)

                             gTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCa                                                            <  1:6760858/93‑1 (MQ=255)

                             gTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCa                                                            <  1:549113/93‑1 (MQ=255)

                              tCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAg                                                           <  2:481657/93‑1 (MQ=255)

                              tCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAg                                                           <  1:2226480/93‑1 (MQ=255)

                                 gTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTg                                                        <  2:4118938/93‑1 (MQ=255)

                                 gTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTg                                                        >  1:2056865/1‑93 (MQ=255)

                                   aaaTATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGaa                                                      <  1:5097759/93‑1 (MQ=255)

                                    aaTATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGaaa                                                     <  1:3330280/93‑1 (MQ=255)

                                    aaTATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGaaa                                                     >  2:2481435/1‑93 (MQ=255)

                                      tataGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATc                                                   <  2:4992925/93‑1 (MQ=255)

                                        taGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCa                                                 >  2:1193517/1‑93 (MQ=255)

                                        taGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCa                                                 <  2:865975/93‑1 (MQ=255)

                                        taGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCa                                                 >  1:4424413/1‑93 (MQ=255)

                                          gcgcCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGc                                               <  1:3515799/93‑1 (MQ=255)

                                          gcgcCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGc                                               >  1:3219118/1‑93 (MQ=255)

                                          gcgcCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGc                                               >  1:532857/1‑93 (MQ=255)

                                          gcgcCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGc                                               >  2:3330317/1‑93 (MQ=255)

                                            gcCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCaa                                             >  2:4201273/1‑93 (MQ=255)

                                               aGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAAtt                                          >  2:2704516/1‑93 (MQ=255)

                                               aGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGCGTCAGCTGAAATCCAGCAAAtt                                          >  1:7303095/1‑93 (MQ=255)

                                                gCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCACATTa                                         >  1:4247341/1‑93 (MQ=255)

                                                gCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTa                                         >  1:1070118/1‑93 (MQ=255)

                                                  gTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAAc                                       <  2:1966069/93‑1 (MQ=255)

                                                  gTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAAc                                       >  1:1844644/1‑93 (MQ=255)

                                                   tAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACa                                      >  1:6128525/1‑93 (MQ=255)

                                                   tAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACa                                      >  1:6993938/1‑93 (MQ=255)

                                                   tAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACa                                      <  1:1516483/93‑1 (MQ=255)

                                                    aaaaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAt                                     <  1:6655571/93‑1 (MQ=255)

                                                     aaaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAtt                                    <  1:5786505/93‑1 (MQ=255)

                                                     aaaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAtt                                    >  1:4789135/1‑93 (MQ=255)

                                                      aaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttt                                   >  2:153500/1‑92 (MQ=255)

                                                      aaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttg                                   <  1:1572394/93‑1 (MQ=255)

                                                      aaGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttg                                   <  2:6164838/93‑1 (MQ=255)

                                                       aGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgt                                  <  1:1201772/93‑1 (MQ=255)

                                                       aGCGGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgt                                  <  1:4507646/93‑1 (MQ=255)

                                                        gCTGTCGTGACCCCGCTGCCGCGCAACCGCGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACAttgtt                                 >  2:4579720/1‑93 (MQ=255)

                                                          tGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTgg                               <  1:3764430/93‑1 (MQ=255)

                                                          tGTCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTgg                               >  1:1231312/1‑93 (MQ=255)

                                                            tCGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGGGTTCCAGGAAGGTGGGTCAGCTGAAATCCAGCAAATTAACATGGTTGGTa                             >  1:2687479/1‑93 (MQ=255)

                                                             cGTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTaa                            <  1:2737778/93‑1 (MQ=255)

                                                              gTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAAc                           >  2:3575732/1‑93 (MQ=255)

                                                              gTGACCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAAc                           <  2:5059942/93‑1 (MQ=255)

                                                                ggccccccTGGCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCa                         <  2:5764060/91‑1 (MQ=255)

                                                                 aCCCCGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCAt                        <  1:7179770/93‑1 (MQ=255)

                                                                  ccccGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATg                       >  2:6094938/1‑93 (MQ=255)

                                                                   cccGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGc                      <  2:3260996/93‑1 (MQ=255)

                                                                   cccGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGc                      <  2:2852334/93‑1 (MQ=255)

                                                                   cccGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGc                      <  1:4826173/93‑1 (MQ=255)

                                                                   cccGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGc                      <  2:7801897/93‑1 (MQ=255)

                                                                     cGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCtt                    <  1:2201648/93‑1 (MQ=255)

                                                                     cGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCtt                    >  1:5628509/1‑93 (MQ=255)

                                                                     cGCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCtt                    >  2:2218099/1‑93 (MQ=255)

                                                                      gCTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCttt                   >  1:7889854/1‑93 (MQ=255)

                                                                       cTGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTc                  <  2:5980339/93‑1 (MQ=255)

                                                                        tGCCGCGCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTca                 >  1:5552868/1‑93 (MQ=255)

                                                                           cgcgCAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTcacca              <  2:3090063/93‑1 (MQ=255)

                                                                            gcgcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTcaccac             <  1:5982852/93‑1 (MQ=255)

                                                                             cgcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCaccacc            <  2:7320725/93‑1 (MQ=255)

                                                                             cgcAAACGTGTCGACCTAAAACTGGTGTTCCAGGAAGGTGTGTCCGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCaccacc            <  1:4691368/93‑1 (MQ=255)

                                                                              gcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcg           >  2:5632627/1‑93 (MQ=255)

                                                                              gcAACCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGCGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcg           <  2:4676547/93‑1 (MQ=255)

                                                                                aaCCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcgac         <  2:6445247/93‑1 (MQ=255)

                                                                                 aCCGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACcgacg        >  2:6290463/1‑93 (MQ=255)

                                                                                   cGTGTTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGaa      >  1:653037/1‑93 (MQ=255)

                                                                                    gtgtTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAAc     <  1:2522979/93‑1 (MQ=255)

                                                                                    gtgtTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAAc     <  1:872988/93‑1 (MQ=255)

                                                                                     tgttgtCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACt    <  2:699058/93‑1 (MQ=255)

                                                                                     tgtTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACt    >  1:1816976/1‑93 (MQ=255)

                                                                                     tgtTGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACt    >  2:5588736/1‑93 (MQ=255)

                                                                                       ttGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGa  <  1:2615327/93‑1 (MQ=255)

                                                                                       ttGGCCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGa  <  1:7694837/93‑1 (MQ=255)

AGAAAGGTCTGGAAGACTTCTACTACAGCGTCGGTAAATATAGCGCCAGCGTAAAAGCTGTCGTGACCCCGCTGCCGCGCAACCGTGTTGACCTAAAACTGGTGTTCCAGGAAGGTGTGTCAGCTGAAATCCAGCAAATTAACATTGTTGGTAACCATGCTTTCACCACCGACGAACTGA  >  NC_000913/198328‑198507

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