BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,735,415	G→A	R97H (CGC→CAC)	malS →	alpha‑amylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,735,415	0	G	A	98.4%	1250.1 / NA	384	R97H (CGC→CAC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base G (0/1); major base A (193/184); minor base C (1/5); total (194/190)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 396 total aligned reads displayed.

ATCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735325‑3735504

aTCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACt                                                                                         <  1:7452331/93‑1 (MQ=255)

 tCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTc                                                                                        >  1:5614544/1‑93 (MQ=255)

 tCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTc                                                                                        >  1:3063903/1‑93 (MQ=255)

  cAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTcc                                                                                       >  2:758550/1‑93 (MQ=255)

    gATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCgg                                                                                     >  1:4542722/1‑93 (MQ=255)

      tGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTa                                                                                   <  1:6534706/93‑1 (MQ=255)

       gCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTAc                                                                                  >  1:5307520/1‑93 (MQ=255)

        cTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACg                                                                                 >  2:4432557/1‑93 (MQ=255)

         tttCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGc                                                                                <  2:6321229/93‑1 (MQ=255)

              cTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACa                                                                           >  1:1583475/1‑93 (MQ=255)

              cTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACa                                                                           <  1:3766234/93‑1 (MQ=255)

               tGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAt                                                                          <  2:4941220/93‑1 (MQ=255)

               tGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAt                                                                          >  2:4707018/1‑93 (MQ=255)

               tGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAt                                                                          >  2:6729256/1‑93 (MQ=255)

                gCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAtt                                                                         <  1:164579/93‑1 (MQ=255)

                gCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAtt                                                                         <  1:5434225/93‑1 (MQ=255)

                  aaCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtga                                                                       <  2:1554065/93‑1 (MQ=255)

                      ttGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTgatgat                                                                   >  2:5049609/1‑93 (MQ=255)

                       tGTAGCAACGCGCCACCGCAATGGCGAGTGTTCAGGGACGGCGAATATACGCGGCAAATAGACCCCCACTCCGGTACGCCAACATTGATGAtt                                                                  <  2:115123/93‑1 (MQ=255)

                        gTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGAttt                                                                 >  2:3128511/1‑93 (MQ=255)

                          aGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGCTTtcc                                                               >  1:2540765/1‑93 (MQ=255)

                          aGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtcc                                                               <  1:648051/93‑1 (MQ=255)

                           gCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtcca                                                              >  1:6856880/1‑93 (MQ=255)

                             aaCACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatc                                                            >  2:5478073/1‑93 (MQ=255)

                             aaCACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatc                                                            <  2:7233924/93‑1 (MQ=255)

                               cacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatcca                                                          >  2:3781959/1‑93 (MQ=255)

                                acGCCGCCTCAAAGGCGATTGTTCAGGGACGGCGAATATACGCTGCACATAGACACCCACTCCGGTACGCCAACATAGATGATTTCCATCCAg                                                         <  2:7021248/93‑1 (MQ=255)

                                  gccgccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGaa                                                       <  2:1410966/93‑1 (MQ=255)

                                  gccgccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGaa                                                       <  1:1935207/93‑1 (MQ=255)

                                   ccgccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAc                                                      <  2:5246452/93‑1 (MQ=255)

                                     gccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgc                                                    <  1:6241689/93‑1 (MQ=255)

                                      ccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   <  1:2027293/93‑1 (MQ=255)

                                      ccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   <  1:4465562/93‑1 (MQ=255)

                                      ccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   <  2:5609729/93‑1 (MQ=255)

                                      ccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   <  1:4037350/93‑1 (MQ=255)

                                       cTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgccg                                                  >  1:1104095/1‑93 (MQ=255)

                                       cTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgccg                                                  <  2:3042497/93‑1 (MQ=255)

                                        tCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgccgc                                                 <  1:3832542/93‑1 (MQ=255)

                                        tCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgccgc                                                 <  1:4926045/93‑1 (MQ=255)

                                         cAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACgccgcc                                                <  1:5143066/93‑1 (MQ=255)

                                          aaTGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGccgccg                                               >  1:1417481/1‑93 (MQ=255)

                                            tGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTGCATCCAGAACGCCGCCGaa                                             >  1:7693407/1‑93 (MQ=255)

                                            tGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGaa                                             >  1:4797493/1‑93 (MQ=255)

                                             ggCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCTGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAc                                            <  1:1816149/93‑1 (MQ=255)

                                             ggCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAc                                            <  2:909801/93‑1 (MQ=255)

                                              gCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAcc                                           <  1:5823212/93‑1 (MQ=255)

                                              gCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAcc                                           <  2:5944651/93‑1 (MQ=255)

                                               cGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCg                                          <  2:4083956/93‑1 (MQ=255)

                                                gATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCGACGTTGATGATTTCCATCCAGAACGCCGCAGAAACgg                                         >  1:266376/1‑93 (MQ=255)

                                                 aTTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGt                                        <  2:1986444/93‑1 (MQ=255)

                                                  ttGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTa                                       <  1:5753276/93‑1 (MQ=255)

                                                   tGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAg                                      >  1:590058/1‑93 (MQ=255)

                                                     ttCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCa                                    >  1:887910/1‑93 (MQ=255)

                                                      tCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCaa                                   >  1:6958904/1‑93 (MQ=255)

                                                      tCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCaa                                   <  1:2116915/93‑1 (MQ=255)

                                                         gggACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGcc                                <  1:7203615/93‑1 (MQ=255)

                                                         gggACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGcc                                >  2:3989815/1‑93 (MQ=255)

                                                         gggACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGcc                                >  2:3267171/1‑93 (MQ=255)

                                                          ggACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCt                               <  2:769777/93‑1 (MQ=255)

                                                           gACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTg                              <  2:2803868/93‑1 (MQ=255)

                                                           gACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTg                              <  1:1735282/93‑1 (MQ=255)

                                                            aCGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTgg                             <  1:4328572/93‑1 (MQ=255)

                                                             cGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGt                            >  2:2197371/1‑93 (MQ=255)

                                                               gCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTcc                          >  1:3563468/1‑93 (MQ=255)

                                                               gCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTcc                          >  2:3443482/1‑93 (MQ=255)

                                                                 gAATATACGCTGCAAATAGACGCACATTCCGGTGCGTCAACATTAATGATTCCCATCCAGAGCGCCGCCGAACGGGTAGCAAGCCTGGTCCGt                        <  2:2438994/93‑1 (MQ=255)

                                                                 gAATATACGCTGCAAATAGACACCCACTCCGGTACGCCGACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGt                        >  2:802779/1‑93 (MQ=255)

                                                                 gAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGt                        >  1:2866798/1‑93 (MQ=255)

                                                                 gAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGt                        <  2:2332134/93‑1 (MQ=255)

                                                                  aaTATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTg                       >  2:7873627/1‑93 (MQ=255)

                                                                  aaTATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTg                       >  1:1535263/1‑93 (MQ=255)

                                                                   atatACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGa                      >  2:3578067/1‑93 (MQ=255)

                                                                   atatACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGCCCGTGa                      <  1:67189/93‑1 (MQ=255)

                                                                    tatcCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCCCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     <  1:2207167/93‑1 (MQ=255)

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     >  1:6050962/1‑93 (MQ=255)

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     >  1:4234909/1‑93 (MQ=255)

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     >  2:2913194/1‑93 (MQ=255)

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAG‑CTGGGCCGGGAAt                    >  1:1251856/1‑93 (MQ=255)

                                                                     ataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAAt                    >  1:3336332/1‑93 (MQ=255)

                                                                      taCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATg                   <  2:43460/93‑1 (MQ=255)

                                                                      taCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATg                   <  2:6646731/93‑1 (MQ=255)

                                                                      taCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATg                   >  1:6126103/1‑93 (MQ=255)

                                                                       aCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGc                  <  1:7285414/93‑1 (MQ=255)

                                                                        cGCTTCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGcc                 <  2:5657225/93‑1 (MQ=255)

                                                                         gCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGccc                >  1:3476734/1‑93 (MQ=255)

                                                                         gCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGccc                >  1:100919/1‑93 (MQ=255)

                                                                          cTGCAAATAGACACCCACCCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGCACCGGTAGCAAGCCTGGCCCGTGAATGCCCg               <  2:717496/93‑1 (MQ=255)

                                                                           tGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGCAGCAAGCCTGGTCCGTGAATGCCCGa              >  2:7398331/1‑93 (MQ=255)

                                                                            gCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGaa             >  2:1935038/1‑93 (MQ=255)

                                                                             cAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGaaa            <  2:4577213/93‑1 (MQ=255)

                                                                                aTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATgg         <  1:5473425/93‑1 (MQ=255)

                                                                                  aGACGCCCTCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGgga       <  2:1717199/93‑1 (MQ=255)

                                                                                  aGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGgga       <  2:5680815/93‑1 (MQ=255)

                                                                                  aGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGgga       <  1:2437970/93‑1 (MQ=255)

                                                                                   gACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggat      <  2:5832169/93‑1 (MQ=255)

                                                                                     cacCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatgg    <  2:1025993/93‑1 (MQ=255)

                                                                                      ccccACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatgga   <  2:6912502/92‑1 (MQ=255)

                                                                                      ccccACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatgga   <  1:2598365/92‑1 (MQ=255)

                                                                                      acCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatggt   >  1:5672151/1‑92 (MQ=255)

                                                                                      acCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatgga   >  2:6421486/1‑93 (MQ=255)

ATCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735325‑3735504

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