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	99.1%	786.0 / ‑5.8	235	R97H (CGC→CAC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (125/106); total (125/110)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 247 total aligned reads displayed.

ATCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735325‑3735504

 tCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTc                                                                                        <  1:3201985/93‑1 (MQ=255)

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

   aGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCg                                                                                      >  1:628836/1‑93 (MQ=255)

     aTGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGt                                                                                    >  1:1091099/1‑93 (MQ=255)

     aTGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGt                                                                                    >  2:2971863/1‑93 (MQ=255)

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

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

      tGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAAAATACGCTGCAAATAGACACCCACTCCGGTa                                                                                   >  2:1016495/1‑93 (MQ=255)

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

       gCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTAc                                                                                  <  2:1738426/93‑1 (MQ=255)

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

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

        cTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACg                                                                                 <  1:2654235/93‑1 (MQ=255)

          ttCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGcc                                                                               >  1:2309541/1‑93 (MQ=255)

            cccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCaa                                                                             <  2:401555/93‑1 (MQ=255)

            cccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCaa                                                                             <  2:2781237/93‑1 (MQ=255)

            cccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCaa                                                                             <  1:4167158/93‑1 (MQ=255)

             ccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACGCCCACTCCGGTACGCCAAc                                                                            >  1:155616/1‑93 (MQ=255)

             ccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAAc                                                                            >  1:2200089/1‑93 (MQ=255)

             ccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAAc                                                                            <  2:4294438/93‑1 (MQ=255)

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

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

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

                 cAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtg                                                                        >  1:242916/1‑93 (MQ=255)

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

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

                   aCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgat                                                                      >  1:2962570/1‑93 (MQ=255)

                   aCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgat                                                                      >  1:1697977/1‑93 (MQ=255)

                    ccTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgatg                                                                     <  1:3052942/93‑1 (MQ=255)

                    ccTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATAGACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgatg                                                                     >  1:3288495/1‑93 (MQ=255)

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

                      ttGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTgatgat                                                                   <  1:3801059/93‑1 (MQ=255)

                          aGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGTTTCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtcc                                                               >  1:4686131/1‑93 (MQ=255)

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

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

                           gCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtcca                                                              <  1:3406043/93‑1 (MQ=255)

                            cAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccat                                                             <  2:2912551/93‑1 (MQ=255)

                              acacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatcc                                                           >  2:1387974/1‑93 (MQ=255)

                              acacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatcc                                                           >  1:781851/1‑93 (MQ=255)

                                 cgccgcCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGa                                                        >  1:4024592/1‑93 (MQ=255)

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

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

                                    cgccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcg                                                     >  1:1388031/1‑93 (MQ=255)

                                     gccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgc                                                    >  1:2747099/1‑93 (MQ=255)

                                      ccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   >  1:3742051/1‑93 (MQ=255)

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

                                      ccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   >  2:437828/1‑93 (MQ=255)

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

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

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

                                        tCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgccgc                                                 >  2:2874288/1‑93 (MQ=255)

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

                                        tCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgccgc                                                 >  2:1235188/1‑93 (MQ=255)

                                           aTGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGa                                              <  1:3091630/93‑1 (MQ=255)

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

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

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

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

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

                                                  ttGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTa                                       >  1:2636260/1‑93 (MQ=255)

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

                                                    gTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGc                                     >  2:568286/1‑93 (MQ=255)

                                                    gTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGc                                     >  1:1931564/1‑93 (MQ=255)

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

                                                        aGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGc                                 >  2:3228143/1‑93 (MQ=255)

                                                        aGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGc                                 >  1:4378284/1‑93 (MQ=255)

                                                        aGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGc                                 <  1:3559106/93‑1 (MQ=255)

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

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

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

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

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

                                                          ggACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCt                               >  2:3900631/1‑93 (MQ=255)

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

                                                           gACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTg                              >  2:4623449/1‑93 (MQ=255)

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

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

                                                             cGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGt                            <  1:4251507/93‑1 (MQ=255)

                                                              ggCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGACCGCCGCCGAACCGGTAGCAAGCCTGGTc                           <  1:863832/93‑1 (MQ=255)

                                                              ggCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTc                           >  1:4643050/1‑93 (MQ=255)

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

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

                                                                   atatACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTAGTCCGTGa                      >  2:2340078/1‑93 (MQ=255)

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGGACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     >  1:3894029/1‑93 (MQ=255)

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     <  1:708468/93‑1 (MQ=255)

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

                                                                    tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTCCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     <  2:3948433/93‑1 (MQ=255)

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

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

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

                                                                            gCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGaa             <  1:1513265/93‑1 (MQ=255)

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

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

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

                                                                                  aGACACCCACTCCGGTACGCCAACATTAATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGgga       <  1:1624251/93‑1 (MQ=255)

                                                                                   gACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggat      >  2:3452921/1‑93 (MQ=255)

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

                                                                                   gACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggat      >  1:2920960/1‑93 (MQ=255)

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

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

ATCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735325‑3735504

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