BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	308,935	G→A	R106C (CGC→TGC)	yagY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	308,935	0	G	A	98.7%	498.2 / ‑5.2	150	R106C (CGC→TGC)	yagY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (70/77); total (70/80)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 154 total aligned reads displayed.

GGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGC  >  NC_000913/308845‑309023

gggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCAct                                                                                        >  2:3680942/1‑93 (MQ=255)

gggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCAct                                                                                        >  1:969815/1‑93 (MQ=255)

gggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCAct                                                                                        >  2:1687305/1‑93 (MQ=255)

 ggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActc                                                                                       <  1:3876142/93‑1 (MQ=255)

 ggcggAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCActc                                                                                       >  1:2106494/1‑93 (MQ=255)

  gcggAAGCAGTGGCCACAGCGCTGCGTTTTGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCCctct                                                                                      <  2:1621768/93‑1 (MQ=255)

    ggAGGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGTATCACTGAGGGCCTGATCAAACCAGACAATGCTGTAGTAGCACTCttt                                                                                    <  2:3379838/93‑1 (MQ=255)

     gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   <  1:1127183/93‑1 (MQ=255)

     gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   >  2:75544/1‑93 (MQ=255)

     gAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCtttt                                                                                   <  1:3985202/93‑1 (MQ=255)

      aaGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCttttt                                                                                  <  1:2482027/93‑1 (MQ=255)

       aGCATTGGCCTTAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACTAGACAATGCGGTAGTAGCACTCTTTTTc                                                                                 >  2:4362813/1‑93 (MQ=255)

        gCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCa                                                                                >  2:1205161/1‑93 (MQ=255)

         cAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCAt                                                                               <  2:3415380/93‑1 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTCGCACTCTTTTTCAtt                                                                              >  1:300704/1‑92 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                              <  2:602758/93‑1 (MQ=255)

          aGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATc                                                                              <  2:2174872/93‑1 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             >  1:3528137/1‑93 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             <  1:3383003/93‑1 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             >  1:2696322/1‑93 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             >  1:1186789/1‑93 (MQ=255)

           gTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCt                                                                             <  2:669785/93‑1 (MQ=255)

             ggCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGc                                                                           <  1:2841466/93‑1 (MQ=255)

              gCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCg                                                                          <  1:829757/93‑1 (MQ=255)

               ccACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCgg                                                                         <  2:4439849/93‑1 (MQ=255)

                 acaGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGt                                                                       >  2:3226046/1‑93 (MQ=255)

                    gcgcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTccc                                                                    <  2:4749014/93‑1 (MQ=255)

                    gcgcTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTccc                                                                    >  2:1305049/1‑93 (MQ=255)

                        tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                <  1:2177999/93‑1 (MQ=255)

                        tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                <  1:92417/93‑1 (MQ=255)

                        tGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtat                                                                >  1:516615/1‑93 (MQ=255)

                         gCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGGCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTtata                                                               <  2:2009770/93‑1 (MQ=255)

                          cGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATag                                                              <  2:228662/93‑1 (MQ=255)

                           ggTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATaga                                                             <  1:2165829/93‑1 (MQ=255)

                            gTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTCTagag                                                            >  1:1370143/1‑92 (MQ=255)

                              tGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaaga                                                          <  1:1117642/93‑1 (MQ=255)

                              tGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATagaaga                                                          <  1:4129005/93‑1 (MQ=255)

                               ggCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAgaagaa                                                         <  1:1550208/93‑1 (MQ=255)

                                gCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGaagaag                                                        <  2:1438146/93‑1 (MQ=255)

                                 cATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                       >  2:719706/1‑93 (MQ=255)

                                 cATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                       <  2:952341/93‑1 (MQ=255)

                                 cATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGc                                                       >  1:1449236/1‑93 (MQ=255)

                                  aTTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCg                                                      >  1:2246401/1‑93 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGCCAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    <  1:4076890/93‑1 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    <  2:2550888/93‑1 (MQ=255)

                                    tATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGa                                                    <  2:1257329/93‑1 (MQ=255)

                                      tCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATc                                                  <  1:4462810/93‑1 (MQ=255)

                                        gcgcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCAc                                                <  2:1545880/93‑1 (MQ=255)

                                          gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACtt                                              <  2:330956/93‑1 (MQ=255)

                                          gcTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACtt                                              <  1:3483353/93‑1 (MQ=255)

                                            tGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCg                                            >  1:3305579/1‑93 (MQ=255)

                                            tGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCg                                            >  2:214209/1‑93 (MQ=255)

                                             gcgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGc                                           <  2:4079081/93‑1 (MQ=255)

                                              cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                          <  1:3715320/93‑1 (MQ=255)

                                              cgcATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCt                                          <  1:3587591/93‑1 (MQ=255)

                                                 aTCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGc                                       >  2:2298811/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGACGAAGCGGATCACTTCGCTGGCt                                      >  2:4098945/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                      <  1:1417229/93‑1 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCt                                      >  1:4658274/1‑93 (MQ=255)

                                                  tCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATAACTTCGCTGGCt                                      >  2:184167/1‑93 (MQ=255)

                                                   cACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCtt                                     >  2:4169412/1‑93 (MQ=255)

                                                    aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                    >  1:3036914/1‑93 (MQ=255)

                                                    aCTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTg                                    >  1:4724790/1‑93 (MQ=255)

                                                     cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                   <  1:2103937/93‑1 (MQ=255)

                                                     cTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTgg                                   <  2:1177070/93‑1 (MQ=255)

                                                      tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  >  2:1105868/1‑93 (MQ=255)

                                                      tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  <  1:3831715/93‑1 (MQ=255)

                                                      tGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTggg                                  <  1:4001994/93‑1 (MQ=255)

                                                         gggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                               >  1:3620727/1‑93 (MQ=255)

                                                         gggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCgg                               >  1:4238522/1‑93 (MQ=255)

                                                          ggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                              <  2:4619009/93‑1 (MQ=255)

                                                          ggCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCggg                              >  2:1146826/1‑93 (MQ=255)

                                                           gCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGt                             <  1:3484884/93‑1 (MQ=255)

                                                           gCCTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGt                             >  1:2343793/1‑93 (MQ=255)

                                                            ccTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTa                            <  1:436912/93‑1 (MQ=255)

                                                             cTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTag                           >  1:521096/1‑93 (MQ=255)

                                                             cTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTag                           >  2:491221/1‑93 (MQ=255)

                                                             cTGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATATGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGATag                           >  2:418068/1‑93 (MQ=255)

                                                              tGATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagc                          >  1:1544959/1‑93 (MQ=255)

                                                               gATCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagca                         <  1:1278735/93‑1 (MQ=255)

                                                                aTCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcag                        <  1:2455188/93‑1 (MQ=255)

                                                                aTCAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTagcag                        >  2:3119310/1‑93 (MQ=255)

                                                                  cAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAgcagca                      <  1:1873711/93‑1 (MQ=255)

                                                                  cAAACCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAgcagca                      <  2:136869/93‑1 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     <  1:2783233/93‑1 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     <  2:368382/93‑1 (MQ=255)

                                                                   aaaCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCaa                     <  2:2593927/93‑1 (MQ=255)

                                                                     aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                   <  1:583723/93‑1 (MQ=255)

                                                                     aCCAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                   >  2:2051073/1‑93 (MQ=255)

                                                                     aCCAGACAATGCGGTAGTAGCACTCTTTTCCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGc                   <  2:4569618/93‑1 (MQ=255)

                                                                      ccAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCt                  <  1:727410/93‑1 (MQ=255)

                                                                      ccAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCt                  <  1:1894530/93‑1 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 >  1:416003/1‑93 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 <  1:1174522/93‑1 (MQ=255)

                                                                       cAGACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTg                 >  2:4610759/1‑93 (MQ=255)

                                                                         gACAATGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGc               >  2:1379995/1‑93 (MQ=255)

                                                                             aTGCGGTAGTAGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGga           <  1:4463230/93‑1 (MQ=255)

                                                                                  gtagtaGCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGgagtgag      <  2:4238522/93‑1 (MQ=255)

                                                                                    agtagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTa    <  1:3455786/93‑1 (MQ=255)

                                                                                      tagCACTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGc  >  2:4237751/1‑93 (MQ=255)

GGGCGGAAGCAGTGGCCACAGCGCTGCGGTTGGCATTATCGCGCTGCGCATCACTGAGGGCCTGATCAAACCAGACAATGCGGTAGTAGCGCTCTTTTTCATCTGCGGGTCCCTTATAGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGC  >  NC_000913/308845‑309023

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