BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,715,348	A→G	D141G (GAC→GGC)	yiaE →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,715,348	0	A	G	86.2%	779.7 / 108.4	298	D141G (GAC→GGC)	yiaE	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (21/20); new base G (119/138); total (140/158)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 310 total aligned reads displayed.

CTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCT  >  NC_000913/3715257‑3715440

cTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGAc                                                                                             >  1:5857186/1‑93 (MQ=255)

 tGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCt                                                                                            >  1:4216660/1‑93 (MQ=255)

 tGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCt                                                                                            <  2:5368732/93‑1 (MQ=255)

  gATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTg                                                                                           >  1:6637554/1‑93 (MQ=255)

  gATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTg                                                                                           <  2:6159259/93‑1 (MQ=255)

   aTGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTgg                                                                                          >  2:1682085/1‑93 (MQ=255)

   aTGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTgg                                                                                          <  1:4345739/93‑1 (MQ=255)

      gcgcTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTAc                                                                                       <  2:2054598/93‑1 (MQ=255)

      gcgcTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTAc                                                                                       <  1:5313942/93‑1 (MQ=255)

      gcgcTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTAc                                                                                       <  1:5960294/93‑1 (MQ=255)

        gcTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACgg                                                                                     >  1:2892781/1‑93 (MQ=255)

         cTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGc                                                                                    <  2:374296/93‑1 (MQ=255)

           ggTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCAc                                                                                  <  1:429001/93‑1 (MQ=255)

            gtgtTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAAGGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACt                                                                                 <  1:1858224/93‑1 (MQ=255)

              gttgtCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGa                                                                               <  2:624796/93‑1 (MQ=255)

               ttgtCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGAc                                                                              <  1:5383056/93‑1 (MQ=255)

                 gtCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGt                                                                            <  1:1606748/93‑1 (MQ=255)

                      ccGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGTCCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCAcc                                                                       <  2:3566154/93‑1 (MQ=255)

                        gCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCAccca                                                                     >  1:4732135/1‑91 (MQ=255)

                        gCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCAt                                                                     >  2:4748331/1‑93 (MQ=255)

                        gCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCAt                                                                     >  1:3123958/1‑93 (MQ=255)

                         cgcgtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATa                                                                    <  1:1407024/91‑1 (MQ=255)

                          tcgtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATaa                                                                   >  1:4485339/1‑93 (MQ=255)

                           cgtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATaaa                                                                  >  2:2869459/1‑93 (MQ=255)

                             tcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAac                                                                <  1:4030025/93‑1 (MQ=255)

                                ggTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACt                                                             >  2:5605505/1‑93 (MQ=255)

                                 gTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTg                                                            <  1:239266/93‑1 (MQ=255)

                                 gTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTg                                                            >  2:4417459/1‑93 (MQ=255)

                                 gTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTg                                                            >  2:1467457/1‑93 (MQ=255)

                                     tgGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCa                                                        <  2:780362/93‑1 (MQ=255)

                                     tgGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCa                                                        <  2:4163384/93‑1 (MQ=255)

                                       gAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCAtt                                                      <  2:4547623/93‑1 (MQ=255)

                                       gAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCAtt                                                      >  2:2013485/1‑93 (MQ=255)

                                       gAGGGAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCAtt                                                      <  1:1341332/93‑1 (MQ=255)

                                        aGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTg                                                     >  2:2397907/1‑93 (MQ=255)

                                          gTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTc                                                   >  1:213456/1‑93 (MQ=255)

                                            aGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCgg                                                 <  2:4251908/93‑1 (MQ=255)

                                            aGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCgg                                                 <  1:2338941/93‑1 (MQ=255)

                                             gCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCggg                                                >  2:6450481/1‑93 (MQ=255)

                                              cAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCCCTGACGTTCACCATAAAACACTGGGCATTGTCGGGa                                               <  2:120429/93‑1 (MQ=255)

                                               aGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGAt                                              <  1:3213430/93‑1 (MQ=255)

                                               aGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTAGGCATTGTCGTGag                                              >  1:2586400/1‑92 (MQ=255)

                                                gAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGCTg                                             >  2:1013044/1‑93 (MQ=255)

                                                gAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATg                                             <  1:2049433/93‑1 (MQ=255)

                                                gAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATg                                             >  1:3428438/1‑93 (MQ=255)

                                                gAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATg                                             <  1:3578324/93‑1 (MQ=255)

                                                 aaCGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATgg                                            >  1:3902734/1‑93 (MQ=255)

                                                 aaCGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATgg                                            <  2:1315098/93‑1 (MQ=255)

                                                   cGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGa                                          >  2:2236659/1‑93 (MQ=255)

                                                    gggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGAc                                         <  2:72917/93‑1 (MQ=255)

                                                    gggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGAc                                         <  2:3021342/93‑1 (MQ=255)

                                                    gggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGAc                                         >  2:5135943/1‑93 (MQ=255)

                                                    gggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGAc                                         <  1:1132791/93‑1 (MQ=255)

                                                     ggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACg                                        >  2:4957302/1‑93 (MQ=255)

                                                      gTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACgg                                       <  2:3503299/93‑1 (MQ=255)

                                                      gTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACgg                                       >  1:2517069/1‑93 (MQ=255)

                                                       tAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGa                                      <  1:2717662/93‑1 (MQ=255)

                                                        aaaaGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCTTTGTCGGGATGGGACGGAt                                     <  1:6497823/93‑1 (MQ=255)

                                                        aaaaGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGAt                                     <  1:1778733/93‑1 (MQ=255)

                                                        aaaaGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGAt                                     <  1:6342434/93‑1 (MQ=255)

                                                         ataGCAGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATc                                    <  1:6272129/91‑1 (MQ=255)

                                                         aaaGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATc                                    <  2:4010593/93‑1 (MQ=255)

                                                              aGGCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCAt                               >  2:2903633/1‑93 (MQ=255)

                                                               ggCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATg                              >  2:4098691/1‑93 (MQ=255)

                                                               ggCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATg                              >  2:1113291/1‑93 (MQ=255)

                                                                gCGAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATgg                             <  1:3364327/93‑1 (MQ=255)

                                                                  gAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGgcg                           >  1:4862921/1‑93 (MQ=255)

                                                                  gAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGgcg                           >  1:373557/1‑93 (MQ=255)

                                                                  gAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGgcg                           >  1:771810/1‑93 (MQ=255)

                                                                  gAATGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGgcg                           >  1:6084985/1‑93 (MQ=255)

                                                                   aaTGGACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGgcgc                          <  1:5092457/93‑1 (MQ=255)

                                                                    aTGGACCGCGGGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGAGGGGACGGATCGGCATGGctct                         <  1:1565155/93‑1 (MQ=255)

                                                                      ggACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCCGCATGGCGCTgg                       >  2:5998264/1‑93 (MQ=255)

                                                                       gACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCGTGGCGCTGGc                      >  1:3292036/1‑93 (MQ=255)

                                                                       gACCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGc                      <  1:704128/93‑1 (MQ=255)

                                                                       gACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGc                      >  1:5496651/1‑93 (MQ=255)

                                                                        aCCGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGAATGGCGCTGGca                     >  2:5304835/1‑93 (MQ=255)

                                                                         ccGCGAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGcac                    <  1:5548592/93‑1 (MQ=255)

                                                                           cggAGCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACaa                  <  1:4395137/91‑1 (MQ=255)

                                                                               gCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTg              >  2:280124/1‑93 (MQ=255)

                                                                               gCATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTg              >  1:6378527/1‑93 (MQ=255)

                                                                                cATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgc             <  1:4817783/93‑1 (MQ=255)

                                                                                cATAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgc             >  1:2981676/1‑93 (MQ=255)

                                                                                 aTAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgcg            <  2:6219403/93‑1 (MQ=255)

                                                                                  tAGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgcgc           >  2:5794217/1‑93 (MQ=255)

                                                                                  tAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgcgc           >  2:1918039/1‑93 (MQ=255)

                                                                                   aGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCa          >  2:3677659/1‑93 (MQ=255)

                                                                                   aGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCa          <  2:3155801/93‑1 (MQ=255)

                                                                                   aGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCa          <  2:1595060/93‑1 (MQ=255)

                                                                                   aGGCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCa          <  2:1556544/93‑1 (MQ=255)

                                                                                    ggCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCAc         >  2:3262576/1‑93 (MQ=255)

                                                                                    ggCCCGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCAc         >  2:2626594/1‑93 (MQ=255)

                                                                                      gccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACtt       <  2:6528690/92‑1 (MQ=255)

                                                                                      cccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACtt       >  1:4604853/1‑93 (MQ=255)

                                                                                      cccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACtt       >  2:849903/1‑93 (MQ=255)

                                                                                       ccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACttt      >  2:4835145/1‑93 (MQ=255)

                                                                                       ccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACttt      <  2:380303/93‑1 (MQ=255)

                                                                                       ccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACttt      <  2:256837/93‑1 (MQ=255)

                                                                                       ccGGGCTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACttt      <  2:1621871/93‑1 (MQ=255)

                                                                                          gACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGc   <  1:4175430/93‑1 (MQ=255)

CTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCT  >  NC_000913/3715257‑3715440

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