BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,735,666	T→C	S181P (TCC→CCC)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,735,666	0	T	C	84.1%	747.6 / 132.5	333	S181P (TCC→CCC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base T (25/28); new base C (149/131); total (174/159)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 341 total aligned reads displayed.

ATTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTACGGA  >  NC_000913/3735576‑3735758

aTTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACAccc                                                                                            >  2:3848975/1‑93 (MQ=255)

 ttGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCg                                                                                           <  2:769340/93‑1 (MQ=255)

 ttGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCg                                                                                           >  2:7217131/1‑93 (MQ=255)

  tGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGc                                                                                          <  1:6305247/93‑1 (MQ=255)

    cGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGccc                                                                                        <  1:4142909/93‑1 (MQ=255)

     gATAGTGAAGAACGGTCAAATAATGTTACACCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGcccc                                                                                       <  1:1618761/93‑1 (MQ=255)

     gATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGcccc                                                                                       >  1:1967385/1‑93 (MQ=255)

     gATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGcccc                                                                                       >  1:3146725/1‑93 (MQ=255)

      aTAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCt                                                                                      >  2:7387952/1‑93 (MQ=255)

      aTAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCt                                                                                      <  1:3825694/93‑1 (MQ=255)

       tAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCtt                                                                                     <  2:4011926/93‑1 (MQ=255)

        aGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCttt                                                                                    <  1:1539858/93‑1 (MQ=255)

         gTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTc                                                                                   <  1:5429463/93‑1 (MQ=255)

           gaagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGa                                                                                 >  1:5222079/1‑93 (MQ=255)

            aagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGAc                                                                                >  2:5030674/1‑93 (MQ=255)

            aagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGAc                                                                                <  2:7124458/93‑1 (MQ=255)

            aagaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGAc                                                                                >  2:581161/1‑93 (MQ=255)

              gaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTg                                                                              <  1:78466/93‑1 (MQ=255)

              gaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTg                                                                              >  1:3286077/1‑93 (MQ=255)

              gaaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTg                                                                              >  2:3932215/1‑93 (MQ=255)

               aaCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTgg                                                                             >  1:4270424/1‑93 (MQ=255)

                aCGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGc                                                                            >  2:2434410/1‑93 (MQ=255)

                  ggTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCAt                                                                          >  2:5655152/1‑93 (MQ=255)

                   gTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATa                                                                         >  1:3958546/1‑93 (MQ=255)

                    tCAAATAATGTTACACCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                        <  1:6889264/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATaa                                                                        <  2:2830344/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                        >  2:3675708/1‑93 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                        <  1:2842542/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                        <  2:673152/93‑1 (MQ=255)

                    tCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATaa                                                                        >  1:5315731/1‑93 (MQ=255)

                     cAAATAATGTTACACCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                       <  1:2637056/93‑1 (MQ=255)

                     cAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                       >  2:6914128/1‑93 (MQ=255)

                     cAAATAATGTTACAACCCGCTGCCACCAACAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAAc                                                                       <  1:1046332/93‑1 (MQ=255)

                      aaaTAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACg                                                                      >  1:2329514/1‑93 (MQ=255)

                      aaaTAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACg                                                                      >  1:5625125/1‑93 (MQ=255)

                       aataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                     <  1:1977498/93‑1 (MQ=255)

                       aataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                     >  1:5906501/1‑93 (MQ=255)

                       aataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGc                                                                     >  2:3323829/1‑93 (MQ=255)

                        ataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGcc                                                                    >  2:3290993/1‑93 (MQ=255)

                        ataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGcc                                                                    >  1:7092057/1‑93 (MQ=255)

                        ataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGcc                                                                    >  1:7014177/1‑93 (MQ=255)

                        ataatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGcc                                                                    >  1:3326467/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCa                                                                   <  2:488021/93‑1 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCa                                                                   <  2:3116977/93‑1 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                   >  2:1740565/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                   >  2:2286312/1‑93 (MQ=255)

                         taatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCa                                                                   <  1:6027261/93‑1 (MQ=255)

                          aatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCAc                                                                  <  2:4262471/93‑1 (MQ=255)

                          aatGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCAc                                                                  <  1:5095830/93‑1 (MQ=255)

                           atGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACg                                                                 >  2:2543931/1‑93 (MQ=255)

                           atGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACg                                                                 >  2:6205957/1‑93 (MQ=255)

                            tGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACgg                                                                >  1:4881200/1‑93 (MQ=255)

                              ttACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGtt                                                              >  2:597836/1‑93 (MQ=255)

                                aCAACCCGCTGCCACCAGCAACGGCTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTa                                                            >  1:1664561/1‑93 (MQ=255)

                                   aCCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACtt                                                         <  1:2061446/93‑1 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                      >  1:1016942/1‑93 (MQ=255)

                                      cGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTtgt                                                      >  2:6241971/1‑93 (MQ=255)

                                        cTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGc                                                    >  2:6616314/1‑93 (MQ=255)

                                         tGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCt                                                   <  1:4788263/93‑1 (MQ=255)

                                         tGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCt                                                   >  2:5952174/1‑93 (MQ=255)

                                         tGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCt                                                   >  1:4478937/1‑93 (MQ=255)

                                          gCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTg                                                  <  1:1389579/93‑1 (MQ=255)

                                          gCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTg                                                  <  2:4775646/93‑1 (MQ=255)

                                           ccaccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGa                                                 <  1:2762094/93‑1 (MQ=255)

                                              ccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAg                                              >  2:5262701/1‑93 (MQ=255)

                                              ccaGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAg                                              <  2:1843710/93‑1 (MQ=255)

                                               caGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGa                                             >  1:2984611/1‑93 (MQ=255)

                                                aGCAACGGTTTCCTCCTGCTGGAACGGGCAGAAACTGACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGAt                                            <  1:2498561/93‑1 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATc                                           >  1:3594019/1‑93 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATc                                           >  2:1081860/1‑93 (MQ=255)

                                                 gCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATc                                           <  2:1377958/93‑1 (MQ=255)

                                                      ggTTTACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTc                                      <  1:2782191/93‑1 (MQ=255)

                                                       gTTTACTCCTGCTGGAACGGGCAGAAACTGACCCCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCg                                     <  1:3746731/93‑1 (MQ=255)

                                                        tttACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGa                                    >  1:4618571/1‑93 (MQ=255)

                                                        tttACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGa                                    <  1:6424614/93‑1 (MQ=255)

                                                          tACTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGaaa                                  <  1:3469425/93‑1 (MQ=255)

                                                            cTCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAAc                                >  2:5464450/1‑93 (MQ=255)

                                                             tCCTGCTGGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACg                               >  1:5278863/1‑93 (MQ=255)

                                                               ctgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGc                             >  1:1822114/1‑93 (MQ=255)

                                                               ctgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGc                             <  1:4838141/93‑1 (MQ=255)

                                                                tgctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCg                            >  2:2719004/1‑93 (MQ=255)

                                                                  ctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGAt                          <  2:7008822/93‑1 (MQ=255)

                                                                  ctgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGAt                          >  2:603319/1‑93 (MQ=255)

                                                                   tgGAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATc                         >  1:6835161/1‑93 (MQ=255)

                                                                    ggAACGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATcc                        >  2:4117089/1‑93 (MQ=255)

                                                                      aaCGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCa                      <  1:7201291/93‑1 (MQ=255)

                                                                       aCGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAg                     <  1:823215/93‑1 (MQ=255)

                                                                        cGGGCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGt                    <  2:4672181/93‑1 (MQ=255)

                                                                          ggCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTaa                  >  2:1184946/1‑93 (MQ=255)

                                                                           gCAGAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCCCGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAAt                 <  2:7253769/93‑1 (MQ=255)

                                                                             aGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGa               <  2:2654883/93‑1 (MQ=255)

                                                                              gAAACTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGAc              <  2:6373996/93‑1 (MQ=255)

                                                                                  cTGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCaga          <  2:432123/93‑1 (MQ=255)

                                                                                   tGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCagag         >  2:4547543/1‑93 (MQ=255)

                                                                                   tGACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCagag         >  1:2171469/1‑93 (MQ=255)

                                                                                    gACACCCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGt        <  1:6610413/93‑1 (MQ=255)

                                                                                    gACACACCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGt        <  1:4397240/93‑1 (MQ=255)

                                                                                      cacaCCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTa      <  1:1757451/93‑1 (MQ=255)

                                                                                       acacCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAc     <  1:6136796/93‑1 (MQ=255)

                                                                                       acacCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTAc     <  2:1005685/93‑1 (MQ=255)

ATTGCGATAGTGAAGAACGGTCAAATAATGTTACAACCCGCTGCCACCAGCAACGGTTTACTCCTGCTGGAACGGGCAGAAACTGACACATCCGCCCCTTTCGACTGGCATAACGCCACGGTTTACTTTGTGCTGACAGATCGTTTCGAAAACGGCGATCCCAGTAATGACCAGAGTTACGGA  >  NC_000913/3735576‑3735758

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