BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	708,421	G→A	G289S (GGC→AGC)	ybfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	708,421	0	G	A	99.0%	1033.1 / ‑2.4	290	G289S (GGC→AGC)	ybfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/1); new base A (157/130); total (159/131)
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.96e-01

Only 100 of 299 total aligned reads displayed.

GATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAC  >  NC_000913/708331‑708510

 atatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCa                                                                                        >  1:5034899/1‑93 (MQ=255)

  tatCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCAc                                                                                       <  2:5749404/93‑1 (MQ=255)

     cGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGc                                                                                    <  2:3063558/93‑1 (MQ=255)

     cGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGc                                                                                    >  2:6526954/1‑93 (MQ=255)

       ccGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGccg                                                                                  >  1:6816497/1‑92 (MQ=255)

       ccGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCt                                                                                  <  1:974933/93‑1 (MQ=255)

       ccGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCt                                                                                  >  2:2110957/1‑93 (MQ=255)

       ccGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCt                                                                                  >  2:6468807/1‑93 (MQ=255)

         ggTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTgg                                                                                >  1:1004937/1‑93 (MQ=255)

          gTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGc                                                                               <  2:3200076/93‑1 (MQ=255)

           tAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCt                                                                              >  1:4567992/1‑93 (MQ=255)

            aGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTg                                                                             >  1:4504390/1‑93 (MQ=255)

              cccGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCa                                                                           >  2:2215770/1‑93 (MQ=255)

              cccGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCa                                                                           >  2:7757423/1‑93 (MQ=255)

               ccGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCaa                                                                          >  1:150103/1‑93 (MQ=255)

                cGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAg                                                                         <  2:4693777/93‑1 (MQ=255)

                cGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAg                                                                         >  1:4382384/1‑93 (MQ=255)

                 gTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGc                                                                        <  2:7389811/93‑1 (MQ=255)

                  ttAACCCCCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCg                                                                       <  2:5366823/93‑1 (MQ=255)

                    aaCCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGtt                                                                     >  1:6080190/1‑93 (MQ=255)

                     accaccAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTg                                                                    >  1:6174441/1‑93 (MQ=255)

                      ccaccaGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTgc                                                                   >  1:292333/1‑92 (MQ=255)

                         ccaGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCt                                                                >  2:6076054/1‑93 (MQ=255)

                           aGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCttt                                                              <  2:6964982/93‑1 (MQ=255)

                           aGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCttt                                                              <  2:1992031/93‑1 (MQ=255)

                             cTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTgg                                                            <  2:3504049/93‑1 (MQ=255)

                              tACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGt                                                           <  1:7905058/93‑1 (MQ=255)

                              tACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGt                                                           >  1:4964156/1‑93 (MQ=255)

                              tACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGt                                                           >  1:5904161/1‑93 (MQ=255)

                                 cAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAc                                                        >  1:7416464/1‑93 (MQ=255)

                                  aGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAcc                                                       >  2:7662322/1‑93 (MQ=255)

                                  aGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTAcc                                                       >  1:8354962/1‑93 (MQ=255)

                                   gTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCg                                                      >  1:6972895/1‑93 (MQ=255)

                                   gTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCg                                                      >  1:2659746/1‑93 (MQ=255)

                                    ttCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCgg                                                     <  2:5506676/93‑1 (MQ=255)

                                     tCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCggg                                                    >  2:5707948/1‑93 (MQ=255)

                                      cTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggc                                                   <  2:6024060/93‑1 (MQ=255)

                                       tACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggcg                                                  >  1:5452504/1‑93 (MQ=255)

                                       tACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggcg                                                  <  1:4942383/93‑1 (MQ=255)

                                         cGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggcggc                                                <  1:7103341/93‑1 (MQ=255)

                                         cGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGggcggc                                                <  1:8632487/93‑1 (MQ=255)

                                          ggTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCt                                               >  1:4847537/1‑93 (MQ=255)

                                          ggTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCt                                               <  2:8335779/93‑1 (MQ=255)

                                           gTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTg                                              >  2:4896102/1‑93 (MQ=255)

                                           gTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTg                                              <  1:2138236/93‑1 (MQ=255)

                                            tACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGa                                             <  2:1367672/93‑1 (MQ=255)

                                            tACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGa                                             >  2:113591/1‑93 (MQ=255)

                                              cccGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACg                                           >  2:4976677/1‑93 (MQ=255)

                                              cccGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACg                                           <  1:8721316/93‑1 (MQ=255)

                                              cccGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTACGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACg                                           <  1:2744691/93‑1 (MQ=255)

                                               ccGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGt                                          >  1:1031168/1‑93 (MQ=255)

                                                cGAGATAAAGTTGACGCTCGCACCGTCAACGAACTTTATGACTGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGATGtg                                         >  2:3913248/1‑92 (MQ=255)

                                                cGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTa                                         <  2:4570879/93‑1 (MQ=255)

                                                 gagaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGCCGTAg                                        >  1:7565332/1‑93 (MQ=255)

                                                  agaTAAAGTTGATGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTTCTAGCGTTGACTTTTGGTTATCTGGCGGCCGACggtgt                                       >  2:7840219/1‑89 (MQ=255)

                                                  agaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAgg                                       >  1:852126/1‑92 (MQ=255)

                                                  agaTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGt                                       >  2:3334348/1‑93 (MQ=255)

                                                    aTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTgg                                     <  1:5483446/93‑1 (MQ=255)

                                                    aTAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTgg                                     <  1:1139751/93‑1 (MQ=255)

                                                     tAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTggt                                    >  1:4137527/1‑92 (MQ=255)

                                                     tAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGa                                    >  1:2374902/1‑93 (MQ=255)

                                                     tAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGa                                    <  1:5309977/93‑1 (MQ=255)

                                                      aaaGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGAt                                   >  1:3117763/1‑93 (MQ=255)

                                                        aGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGAttt                                 <  1:2225524/93‑1 (MQ=255)

                                                        aGTTGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGAttt                                 >  2:719670/1‑93 (MQ=255)

                                                          ttGACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTgc                               <  1:7742870/93‑1 (MQ=255)

                                                            gACGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTgcgc                             <  2:8534946/93‑1 (MQ=255)

                                                              cGATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCt                           >  2:5082132/1‑93 (MQ=255)

                                                              cGATCGCAGCGTCAACGACCTTGATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCt                           <  1:3228855/93‑1 (MQ=255)

                                                               gATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTc                          <  1:963639/93‑1 (MQ=255)

                                                               gATCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTc                          <  1:3169575/93‑1 (MQ=255)

                                                                aTCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCg                         >  1:4755950/1‑93 (MQ=255)

                                                                aTCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCg                         >  1:3655990/1‑93 (MQ=255)

                                                                aTCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCg                         >  2:2498317/1‑93 (MQ=255)

                                                                aTCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCg                         >  2:2578114/1‑93 (MQ=255)

                                                                 tCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGa                        <  2:6775941/93‑1 (MQ=255)

                                                                 tCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGa                        >  1:1999257/1‑93 (MQ=255)

                                                                 tCGCAGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTCTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGa                        >  1:39702/1‑93 (MQ=255)

                                                                   gCAGAGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAg                      >  2:3707816/1‑93 (MQ=255)

                                                                    cAGCGTCACCGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAgg                     <  2:726784/93‑1 (MQ=255)

                                                                     aGCGTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGc                    <  1:3397876/93‑1 (MQ=255)

                                                                        gTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCAcc                 <  1:7248690/93‑1 (MQ=255)

                                                                        gTCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCAcc                 <  2:520748/93‑1 (MQ=255)

                                                                         tCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCt                >  1:7976328/1‑93 (MQ=255)

                                                                         tCAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCt                >  2:8620973/1‑93 (MQ=255)

                                                                          cAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTg               <  2:5865808/93‑1 (MQ=255)

                                                                          cAACGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTg               <  1:1759927/93‑1 (MQ=255)

                                                                             cGACCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGt            >  2:2750497/1‑93 (MQ=255)

                                                                               aCCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTa          >  2:794724/1‑93 (MQ=255)

                                                                               aCCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTa          >  2:5776204/1‑93 (MQ=255)

                                                                               aCCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTa          >  1:7179019/1‑93 (MQ=255)

                                                                               aCCTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTa          >  1:712361/1‑93 (MQ=255)

                                                                                ccTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTaa         >  2:939517/1‑93 (MQ=255)

                                                                                ccTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCCGGGTTaa         >  2:1687528/1‑93 (MQ=255)

                                                                                 cTTTATGACAGCACCGCCTGGCTGCAAGCGTTGACATTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAg        >  1:6743515/1‑93 (MQ=255)

                                                                                  tttATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAgg       >  1:3512927/1‑93 (MQ=255)

                                                                                   ttATGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAATGc      >  1:4767432/1‑93 (MQ=255)

                                                                                     aTGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTg    >  1:6542837/1‑93 (MQ=255)

                                                                                      tGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGa   <  1:6258116/93‑1 (MQ=255)

                                                                                      tGACAGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGa   <  2:2853268/93‑1 (MQ=255)

GATATCGCCGGTAGCCCGTTAACCACCAGCTACCAGTTCTACGGTACCCGAGATAAAGTTGACGATCGCAGCGTCAACGACCTTTATGACGGCACCGCCTGGCTGCAAGCGTTGACCTTTGGTTACCGGGCGGCTGACGTAGTGGATTTGCGCCTCGAAGGCACCTGGGTTAAGGCTGAC  >  NC_000913/708331‑708510

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