BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,520	A→G	S531P (TCC→CCC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,520	0	A	G	98.3%	534.2 / ‑3.2	174	S531P (TCC→CCC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (88/83); total (90/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

TACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594431‑1594609

 acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCa                                                                                       <  1:5835716/93‑1 (MQ=255)

   acCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAtt                                                                                     <  2:2891190/93‑1 (MQ=255)

    ccGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAttt                                                                                    <  1:2332904/93‑1 (MQ=255)

     cGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtc                                                                                   >  1:1653289/1‑93 (MQ=255)

     cGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtc                                                                                   <  1:3613843/93‑1 (MQ=255)

         cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtcctcc                                                                               <  1:6602994/93‑1 (MQ=255)

         cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtcctcc                                                                               <  1:4882845/93‑1 (MQ=255)

          tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCg                                                                              >  1:1007768/1‑93 (MQ=255)

          tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCg                                                                              >  2:5071712/1‑93 (MQ=255)

             aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATa                                                                           >  1:993821/1‑93 (MQ=255)

              gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATaa                                                                          >  1:4420799/1‑93 (MQ=255)

               cAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTCCCGGGGGCATTTCCTCCGATaaa                                                                         <  2:4944493/93‑1 (MQ=255)

               cAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATaaa                                                                         >  2:4131670/1‑93 (MQ=255)

               cAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATaaa                                                                         <  2:3458429/93‑1 (MQ=255)

                aaTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAAt                                                                        >  1:2620503/1‑93 (MQ=255)

                  ttGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATcc                                                                      >  2:1156988/1‑93 (MQ=255)

                  ttGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATcc                                                                      <  1:4237109/93‑1 (MQ=255)

                   tgtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCa                                                                     >  1:2613188/1‑93 (MQ=255)

                   tgtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCa                                                                     <  2:3098274/93‑1 (MQ=255)

                    gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCAt                                                                    <  1:5485618/93‑1 (MQ=255)

                    gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCAt                                                                    <  2:1428620/93‑1 (MQ=255)

                    gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCAt                                                                    <  2:6591704/93‑1 (MQ=255)

                    gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCAt                                                                    <  1:1185963/93‑1 (MQ=255)

                     tgCGTTAGTAGCTGATGTCTTTAAATGCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  2:5896513/93‑1 (MQ=255)

                     tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  1:5123933/93‑1 (MQ=255)

                     tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  1:7149745/93‑1 (MQ=255)

                        gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGcc                                                                >  2:7741167/1‑93 (MQ=255)

                        gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGcc                                                                >  2:800146/1‑93 (MQ=255)

                        gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGcc                                                                >  2:5380758/1‑93 (MQ=255)

                          tagtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCg                                                              >  1:7180411/1‑93 (MQ=255)

                           agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGt                                                             <  2:6209228/93‑1 (MQ=255)

                            gtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                            <  1:3782001/93‑1 (MQ=255)

                            gtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                            >  1:6598484/1‑93 (MQ=255)

                            gtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                            >  2:5500819/1‑93 (MQ=255)

                              agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGa                                                          >  2:2486868/1‑93 (MQ=255)

                              agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGa                                                          >  1:3072984/1‑93 (MQ=255)

                               gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  2:4999981/1‑93 (MQ=255)

                                   atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTaa                                                     <  2:1045173/93‑1 (MQ=255)

                                    tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGACTTAAg                                                    >  1:4720467/1‑93 (MQ=255)

                                    tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  2:5783681/1‑93 (MQ=255)

                                    tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  2:6002336/93‑1 (MQ=255)

                                    tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    >  2:3024049/1‑93 (MQ=255)

                                    tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  1:1085685/93‑1 (MQ=255)

                                    tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  1:3235874/93‑1 (MQ=255)

                                     gTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                   >  2:1633928/1‑93 (MQ=255)

                                     gTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                   <  2:3899596/93‑1 (MQ=255)

                                     gTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAgg                                                   >  1:335826/1‑93 (MQ=255)

                                      tCTTTAAATTCCAGAGACTCCCCGTTGGAAGATCAACCATCCCCCTACCGCGGGCATTTCCTCCGACAAAACCAAAGACCATTAAATTAAgaa                                                  >  2:7317069/1‑91 (MQ=255)

                                      gcTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAggg                                                  <  2:2704463/92‑1 (MQ=255)

                                       cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 <  1:7520462/93‑1 (MQ=255)

                                       cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 <  2:6376964/93‑1 (MQ=255)

                                       cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 <  2:3949474/93‑1 (MQ=255)

                                        tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  1:2371503/1‑93 (MQ=255)

                                        tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  2:1822006/1‑93 (MQ=255)

                                        tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                <  2:5076346/93‑1 (MQ=255)

                                        tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  2:4113669/1‑93 (MQ=255)

                                          tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              >  1:1529968/1‑93 (MQ=255)

                                          tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              <  1:6721240/93‑1 (MQ=255)

                                             aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                           <  1:3682875/93‑1 (MQ=255)

                                                 cAGAGACTGTCCATTGAAATATTAACTATCCCCTTACCGGGGACATTACCTCCGATAAATCCATAGCCCGTTGAAATAAGGGTGGCGCCAttc                                       <  2:7035344/93‑2 (MQ=255)

                                                  agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTAc                                      >  2:5744511/1‑93 (MQ=255)

                                                  agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTAc                                      >  1:7634639/1‑93 (MQ=255)

                                                   gagaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                     >  1:3751634/1‑93 (MQ=255)

                                                     gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  1:4638777/1‑93 (MQ=255)

                                                     gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  1:7033598/1‑93 (MQ=255)

                                                      gcTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                  <  1:2426671/92‑1 (MQ=255)

                                                      aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                  >  2:3245713/1‑93 (MQ=255)

                                                       cTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACa                                 <  2:1180225/93‑1 (MQ=255)

                                                        tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                >  1:361825/1‑93 (MQ=255)

                                                         gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               <  2:1335914/93‑1 (MQ=255)

                                                         gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               <  2:4413134/93‑1 (MQ=255)

                                                          tCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAttt                              <  1:2008234/93‑1 (MQ=255)

                                                           ccGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                             <  1:3530793/93‑1 (MQ=255)

                                                           ccGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                             <  2:1832121/93‑1 (MQ=255)

                                                            cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                            <  1:4117189/93‑1 (MQ=255)

                                                            cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                            >  1:4856855/1‑93 (MQ=255)

                                                             gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                           >  1:137039/1‑93 (MQ=255)

                                                              ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          >  1:4080390/1‑93 (MQ=255)

                                                              ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          <  2:217182/93‑1 (MQ=255)

                                                              ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          <  1:6761161/93‑1 (MQ=255)

                                                                gAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                        >  2:5915469/1‑93 (MQ=255)

                                                                  aaTATTACCTATCCCATTACCGGGGGCATTTCCTCCGATAAATCCATCGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATccc                      <  2:3134844/93‑1 (MQ=255)

                                                                   atatTAACTATCCCCTTACCGGGAGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                     <  2:4274029/93‑1 (MQ=255)

                                                                   atatTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGCTACCg                     >  1:1527794/1‑93 (MQ=255)

                                                                    tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                    >  2:1632922/1‑93 (MQ=255)

                                                                     atTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCt                   >  1:5234171/1‑93 (MQ=255)

                                                                     atTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCt                   >  2:7174948/1‑93 (MQ=255)

                                                                        aaCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAcc                <  1:5782002/93‑1 (MQ=255)

                                                                         aCTCTCCCCTTACCGGGGGCATTGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa               <  1:3812221/93‑1 (MQ=255)

                                                                           tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa             >  2:5738421/1‑93 (MQ=255)

                                                                           tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa             >  1:7405782/1‑93 (MQ=255)

                                                                             tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt           >  1:4778475/1‑93 (MQ=255)

                                                                             tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt           >  1:728374/1‑93 (MQ=255)

                                                                              ccccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGACTTAAGGGTGGCACCATTACTGACATTTACGATACCGCTTCCCAAAtt          >  1:3578356/1‑93 (MQ=255)

                                                                               cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa         >  2:238024/1‑93 (MQ=255)

                                                                                ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc        >  2:2407714/1‑93 (MQ=255)

                                                                                ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc        >  2:4628586/1‑93 (MQ=255)

                                                                                    aCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc    <  1:5432587/93‑1 (MQ=255)

                                                                                     ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   >  1:4603665/1‑93 (MQ=255)

                                                                                     ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   <  2:1915204/93‑1 (MQ=255)

TACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594431‑1594609

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