BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,628,644	A→T	C99S (TGT→AGT)	yfgF ←	cyclic‑di‑GMP phosphodiesterase, anaerobic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,644	0	A	T	98.9%	592.5 / ‑4.2	175	C99S (TGT→AGT)	yfgF	cyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand): ref base A (0/2); new base T (103/69); total (104/71)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 178 total aligned reads displayed.

ACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCC  >  NC_000913/2628548‑2628739

acacAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACtga                                                                                               >  1:100523/1‑99 (MQ=255)

  acaATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTgaga                                                                                             >  1:3919274/1‑99 (MQ=255)

   caATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaa                                                                                            <  1:1206898/99‑1 (MQ=255)

    aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaat                                                                                           <  2:2815472/99‑1 (MQ=255)

    aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaat                                                                                           <  1:75996/99‑1 (MQ=255)

    aaTCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaat                                                                                           <  1:1740289/99‑1 (MQ=255)

      tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatga                                                                                         <  1:863294/99‑1 (MQ=255)

      tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatga                                                                                         >  2:4013295/1‑99 (MQ=255)

      tCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACACCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatga                                                                                         <  2:2145808/99‑1 (MQ=255)

       cTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgat                                                                                        <  1:1443619/99‑1 (MQ=255)

        tGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgata                                                                                       <  1:1182333/99‑1 (MQ=255)

         gCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataa                                                                                      >  2:4643654/1‑99 (MQ=255)

         gCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataa                                                                                      <  1:1970080/99‑1 (MQ=255)

         gCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataa                                                                                      <  1:2729000/99‑1 (MQ=255)

          ccAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataat                                                                                     >  2:2304099/1‑99 (MQ=255)

          ccAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataat                                                                                     <  1:4206784/99‑1 (MQ=255)

          ccAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataat                                                                                     >  2:458470/1‑99 (MQ=255)

           cAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataatg                                                                                    >  2:2581511/1‑99 (MQ=255)

            aGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataatga                                                                                   >  2:2355319/1‑99 (MQ=255)

            aGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGaatgataatga                                                                                   <  1:2058218/99‑1 (MQ=255)

                tataCGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATgataatgataat                                                                               <  1:3314573/99‑1 (MQ=255)

                    cGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATaaaa                                                                           <  1:2343215/99‑1 (MQ=255)

                     gCTGTGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAAt                                                                          <  2:309119/99‑1 (MQ=255)

                     gCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAAt                                                                          >  1:2504575/1‑99 (MQ=255)

                     gCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAAt                                                                          <  2:1777951/99‑1 (MQ=255)

                     gCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAAt                                                                          <  1:905/99‑1 (MQ=255)

                        gggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAg                                                                       <  1:1365041/99‑1 (MQ=255)

                        gggAAATGAAACGGGTATCACCATGTGAAACGTTGCTACGCCGGGGAGCAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAg                                                                       >  2:5299785/1‑99 (MQ=255)

                         ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGc                                                                      >  1:5256366/1‑99 (MQ=255)

                         ggAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGc                                                                      <  1:4530498/99‑1 (MQ=255)

                         ggAAATTAAACGGGGACCACCATATGGAACTTTGTCAGGCCGGGGGGTAAAGAACCTGTAACCGCCCCAACTGAAAAGATTAAGGATAATTACATTAGc                                                                      <  2:5503544/99‑1 (MQ=255)

                             aTTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGtt                                                                  <  1:4469175/99‑1 (MQ=255)

                              ttAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                 >  1:3734267/1‑99 (MQ=255)

                              ttAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTa                                                                 <  2:1768276/99‑1 (MQ=255)

                               tAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTaa                                                                >  1:9707/1‑99 (MQ=255)

                                aaaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAAt                                                               <  2:1195700/99‑1 (MQ=255)

                                aaaCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAAt                                                               <  1:10198/99‑1 (MQ=255)

                                  aggggTATCACCATCTAAAACGTTGTTACGCCAGGAAGTAAAGACCCTGTAACCTCCCCAACTTAGAATGATAATGATAATAAAATTAGCCGTTAATGa                                                             <  2:362162/97‑1 (MQ=255)

                                  aCGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGa                                                             >  2:857832/1‑99 (MQ=255)

                                    gggTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATa                                                           >  1:2933846/1‑99 (MQ=255)

                                          cACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGttt                                                     >  1:3521094/1‑99 (MQ=255)

                                           aCCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTc                                                    <  2:2957866/99‑1 (MQ=255)

                                            ccATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCg                                                   >  2:1632941/1‑99 (MQ=255)

                                             cATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                  >  2:4850901/1‑99 (MQ=255)

                                              aTGTGAAACGTTGTTACTCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCGGTTAATGATAGCGTTTCGTa                                                 <  2:4218178/99‑1 (MQ=255)

                                                gtgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGAAAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaa                                               >  1:2803860/1‑99 (MQ=255)

                                                 tgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa                                              >  1:1867594/1‑99 (MQ=255)

                                                 tgAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaa                                              >  1:5044456/1‑99 (MQ=255)

                                                  gAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTaaaaa                                             <  2:3823251/99‑1 (MQ=255)

                                                   aaaCGTTGTTACGCGGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAc                                            >  2:4194170/1‑99 (MQ=255)

                                                   aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAc                                            >  1:3247321/1‑99 (MQ=255)

                                                   aaaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAc                                            >  2:1521919/1‑99 (MQ=255)

                                                    aaCGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAcc                                           >  1:409142/1‑99 (MQ=255)

                                                       gttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                        <  1:4739225/99‑1 (MQ=255)

                                                        ttgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGc                                       >  2:649074/1‑99 (MQ=255)

                                                         tgttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGct                                      <  2:418768/99‑1 (MQ=255)

                                                           ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGCGTGATAATGATAATAAAAACAGCAGTTAACGATAGCGTTTCGTAAACACCCAGcgct                                    >  2:4495554/1‑99 (MQ=255)

                                                           ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                    <  2:3365047/99‑1 (MQ=255)

                                                           ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                    >  1:4638785/1‑99 (MQ=255)

                                                           ttACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGctct                                    >  1:111060/1‑99 (MQ=255)

                                                             aCGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                  >  1:123835/1‑99 (MQ=255)

                                                             aCGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGc                                  <  2:2977313/99‑1 (MQ=255)

                                                             aCGCCAGGTAGTAAAGACCCTGTAACCGCGCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGAGAGCGTTTCGTAAAAACCCAGCTCTGc                                  <  2:4317553/99‑1 (MQ=255)

                                                              cGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAACATTAGCCGTTAATGATACCGTTTCGTAAGAACCCAGCTCTGCa                                 >  2:4828329/1‑99 (MQ=255)

                                                               gCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCat                                >  1:615780/1‑99 (MQ=255)

                                                               gCCGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCat                                >  1:936527/1‑99 (MQ=255)

                                                                ccGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCata                               <  2:137634/99‑1 (MQ=255)

                                                                 cGGGGAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCatat                              >  1:4131897/1‑99 (MQ=255)

                                                                   gggAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATAttt                            >  1:2865964/1‑99 (MQ=255)

                                                                     gAGTAAAGACCCTGTAACCGCCCCAACTGCGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                          <  2:4153031/99‑1 (MQ=255)

                                                                     gAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                          >  2:798230/1‑99 (MQ=255)

                                                                     gAGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTACCCGATAATGATAGCGTTTCGAAAAAACCCAGCTCTGCATATTTGc                          >  2:5160521/1‑99 (MQ=255)

                                                                     gAGTAAAGACCCTCTAACCGCCCCAACTGAGCATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                          <  2:4325609/99‑1 (MQ=255)

                                                                      aGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                         <  2:1141522/99‑1 (MQ=255)

                                                                      aGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                         >  1:569656/1‑99 (MQ=255)

                                                                      aGTAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCg                         >  1:4423728/1‑99 (MQ=255)

                                                                        tAAAGACCCTGTAAGCGCCCCAACTGAGACTGATAATGATAATAGAATTAGCCGTTGATGATAGCGTTTCGTAAAGACCCAGCTCTGCATATTTGCGGa                       >  1:4903350/1‑99 (MQ=255)

                                                                        tAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                       >  2:4699803/1‑99 (MQ=255)

                                                                        tAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                       >  1:3946419/1‑99 (MQ=255)

                                                                        tAAAGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGa                       <  2:4544646/99‑1 (MQ=255)

                                                                          aaGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACa                     >  2:3297220/1‑99 (MQ=255)

                                                                          aaGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACa                     >  1:3632241/1‑99 (MQ=255)

                                                                          aaGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACa                     >  2:222726/1‑99 (MQ=255)

                                                                           aGACCCTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaa                    >  2:2393096/1‑99 (MQ=255)

                                                                              cccTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaaaaa                 >  1:526117/1‑99 (MQ=255)

                                                                               ccTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                <  2:4425337/99‑1 (MQ=255)

                                                                               ccTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                >  2:4723616/1‑99 (MQ=255)

                                                                               ccTGTAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAAt                <  1:3337468/99‑1 (MQ=255)

                                                                                   tAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt            >  2:2555456/1‑99 (MQ=255)

                                                                                   tAACCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt            >  2:539100/1‑99 (MQ=255)

                                                                                    aaCCGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTa           >  1:5202208/1‑99 (MQ=255)

                                                                                      ccGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAg         >  1:1541206/1‑99 (MQ=255)

                                                                                      ccGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAg         >  1:677106/1‑99 (MQ=255)

                                                                                      ccGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAg         >  1:5481932/1‑99 (MQ=255)

                                                                                       cGCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAgg        <  1:3063051/99‑1 (MQ=255)

                                                                                        gCCCCAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGAAAAAAATCCCTAAGGc       >  2:588572/1‑99 (MQ=255)

                                                                                          cccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCTAGCTCTGCATATTTACGGACAAAAATCCCTAAGGCGa     >  1:5127330/1‑99 (MQ=255)

                                                                                          cccAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGa     >  1:4188108/1‑99 (MQ=255)

                                                                                            cAACTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATc   >  1:2160990/1‑99 (MQ=255)

                                                                                             aaCTGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATcc  <  2:3187887/99‑1 (MQ=255)

ACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGGTATCACCATGTGAAACGTTGTTACGCCGGGGAGTAAAGACCCTGTAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCC  >  NC_000913/2628548‑2628739

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