BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,753,749	A→G	I10V (ATC→GTC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,753,749	0	A	G	99.4%	554.5 / ‑4.6	171	I10V (ATC→GTC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (91/79); total (91/80)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 174 total aligned reads displayed.

TCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGAC  >  NC_000913/1753659‑1753838

tCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTc                                                                                         >  1:2446701/1‑93 (MQ=255)

tCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTc                                                                                         <  2:4103488/93‑1 (MQ=255)

  cGGCAGTCCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCgg                                                                                       <  1:4293475/93‑1 (MQ=255)

   ggCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGa                                                                                      >  1:121573/1‑93 (MQ=255)

      aGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCg                                                                                   <  1:771244/93‑1 (MQ=255)

      aGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCg                                                                                   >  1:1650871/1‑93 (MQ=255)

       gTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGa                                                                                  >  2:4400980/1‑93 (MQ=255)

           gcCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAAc                                                                              >  2:1930320/1‑93 (MQ=255)

            cccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAAcc                                                                             >  2:228808/1‑93 (MQ=255)

            cccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAAcc                                                                             >  1:1905962/1‑93 (MQ=255)

            cccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAAcc                                                                             <  1:4460606/93‑1 (MQ=255)

             ccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAAcca                                                                            >  1:4076567/1‑92 (MQ=255)

             ccAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCg                                                                            >  1:1651965/1‑93 (MQ=255)

              cAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGa                                                                           <  1:318216/93‑1 (MQ=255)

               ataaaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          <  1:3362580/91‑1 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          >  2:1691657/1‑93 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          <  2:3167991/93‑1 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          <  1:1272216/93‑1 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          <  1:1181603/93‑1 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          <  1:3660210/93‑1 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGaa                                                                          <  1:3807564/93‑1 (MQ=255)

               aGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTCGCACCGTCGGACCGAAAACCGaa                                                                          >  2:2593078/1‑93 (MQ=255)

                gAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAAt                                                                         <  2:2247903/93‑1 (MQ=255)

                  aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCt                                                                       <  2:1124346/93‑1 (MQ=255)

                  aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCt                                                                       >  1:2582469/1‑93 (MQ=255)

                  aaGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCt                                                                       <  2:1142456/93‑1 (MQ=255)

                   aGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTg                                                                      >  1:751625/1‑93 (MQ=255)

                    gCAAGTTTCTCCCATCCTTCACAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGa                                                                     >  1:4501324/1‑93 (MQ=255)

                     cAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGag                                                                    >  1:1460244/1‑92 (MQ=255)

                     cAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGaa                                                                    <  2:3548391/93‑1 (MQ=255)

                      aaGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAag                                                                   >  1:2105150/1‑93 (MQ=255)

                       aGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAaga                                                                  >  2:1356974/1‑93 (MQ=255)

                         tttCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAgaga                                                                >  2:806668/1‑93 (MQ=255)

                          ttCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGAt                                                               <  2:1754542/93‑1 (MQ=255)

                          ttCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGAt                                                               <  1:364708/93‑1 (MQ=255)

                           tctcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATg                                                              >  1:2667546/1‑93 (MQ=255)

                            ctcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGt                                                             <  1:558037/93‑1 (MQ=255)

                             tcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGtt                                                            >  1:2422615/1‑93 (MQ=255)

                             tcCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGtt                                                            >  2:1219320/1‑93 (MQ=255)

                              cccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTa                                                           <  1:1815113/93‑1 (MQ=255)

                               ccATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAg                                                          >  2:3006671/1‑93 (MQ=255)

                                cATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGTGCGAAAACCGAATCTGAAGAGATGTTAGc                                                         >  1:3235421/1‑93 (MQ=255)

                                cATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGc                                                         >  1:1547821/1‑93 (MQ=255)

                                  tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTa                                                       <  1:2000011/93‑1 (MQ=255)

                                  tCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTa                                                       >  2:928364/1‑93 (MQ=255)

                                        tcAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGc                                                 >  2:2062721/1‑93 (MQ=255)

                                         cAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCt                                                >  1:1279041/1‑93 (MQ=255)

                                           aCTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTgg                                              >  1:4683818/1‑93 (MQ=255)

                                            cTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGa                                             >  2:3254592/1‑93 (MQ=255)

                                              tAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACg                                           >  2:2875003/1‑93 (MQ=255)

                                              tAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACg                                           >  1:2356268/1‑93 (MQ=255)

                                               aaaGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGc                                          <  1:1171857/93‑1 (MQ=255)

                                                aagactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGGGATGTTAGCTAAAATGCTGGACGCt                                         <  1:413798/93‑1 (MQ=255)

                                                aagactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCt                                         >  2:2688407/1‑93 (MQ=255)

                                                aagactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCt                                         >  2:2557276/1‑93 (MQ=255)

                                                 agactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                                        >  1:4137179/1‑93 (MQ=255)

                                                 agactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTg                                        <  1:1696196/93‑1 (MQ=255)

                                                  gactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTgg                                       <  2:1645442/93‑1 (MQ=255)

                                                  gactaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTgg                                       >  1:3733758/1‑93 (MQ=255)

                                                   actaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGc                                      <  2:3776674/93‑1 (MQ=255)

                                                   actaagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGc                                      >  2:317425/1‑93 (MQ=255)

                                                     taagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACAGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCAt                                    >  2:4183257/1‑93 (MQ=255)

                                                      aagactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATg                                   >  1:2902441/1‑93 (MQ=255)

                                                       agactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGa                                  >  1:458446/1‑93 (MQ=255)

                                                        gactGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGaa                                 >  1:4093840/1‑93 (MQ=255)

                                                         actGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAAc                                >  1:3122269/1‑93 (MQ=255)

                                                          ctGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACg                               >  1:4656742/1‑93 (MQ=255)

                                                           tGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGt                              <  1:378303/93‑1 (MQ=255)

                                                           tGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGt                              >  2:2919235/1‑93 (MQ=255)

                                                           tGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGt                              >  1:1067815/1‑93 (MQ=255)

                                                            ggcATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGtt                             <  2:2472488/91‑1 (MQ=255)

                                                            gTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGtt                             >  2:1416682/1‑93 (MQ=255)

                                                             tCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTa                            <  2:3918043/93‑1 (MQ=255)

                                                             tCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTa                            >  2:4308762/1‑93 (MQ=255)

                                                             tCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTa                            <  1:2290971/93‑1 (MQ=255)

                                                             tCATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTa                            >  1:1157024/1‑93 (MQ=255)

                                                              cATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt                           <  1:203445/93‑1 (MQ=255)

                                                              cATGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTAt                           <  1:2968864/93‑1 (MQ=255)

                                                               aTGAAAAAGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATg                          >  1:550746/1‑93 (MQ=255)

                                                                   aaaaGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTc                      <  2:2883702/93‑1 (MQ=255)

                                                                   aaaaGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTc                      >  2:3791040/1‑93 (MQ=255)

                                                                    aaaGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCt                     >  2:4496328/1‑93 (MQ=255)

                                                                     aaGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTg                    <  1:1085918/93‑1 (MQ=255)

                                                                      aGACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGa                   <  1:4514687/93‑1 (MQ=255)

                                                                       gACCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGaa                  >  2:4705890/1‑93 (MQ=255)

                                                                        aCCAAAATTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAAc                 <  1:1673523/93‑1 (MQ=255)

                                                                           aaaaTTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtc              <  2:881000/93‑1 (MQ=255)

                                                                           aaaaTTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtc              >  2:3833379/1‑93 (MQ=255)

                                                                            aaaTTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtct             >  2:818098/1‑93 (MQ=255)

                                                                             aaTTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtctc            <  2:417234/93‑1 (MQ=255)

                                                                             aaTTGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtctc            >  2:4405963/1‑93 (MQ=255)

                                                                               ttGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTtctctc          >  1:959317/1‑93 (MQ=255)

                                                                                tGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCa         <  1:2372032/93‑1 (MQ=255)

                                                                                tGTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCa         <  1:1131683/93‑1 (MQ=255)

                                                                                 gTTTGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCAt        <  2:2430822/93‑1 (MQ=255)

                                                                                  tttGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATg       <  2:1510993/93‑1 (MQ=255)

                                                                                   ttGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAATATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATgg      <  1:403642/93‑1 (MQ=255)

                                                                                   ttGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATgg      <  2:1556216/93‑1 (MQ=255)

                                                                                    tGCACCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGt     <  1:1744450/93‑1 (MQ=255)

                                                                                       aCCGTCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGAc  >  2:143475/1‑93 (MQ=255)

TCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAAAACCGAATCTGAAGAGATGTTAGCTAAAATGCTGGACGCTGGCATGAACGTTATGCGTCTGAACTTCTCTCATGGTGAC  >  NC_000913/1753659‑1753838

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