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	97.3%	583.3 / 0.4	187	R97H (CGC→CAC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base G (3/2); new base A (81/97); total (85/102)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

Only 100 of 192 total aligned reads displayed.

TCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735326‑3735504

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

   gATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCgg                                                                                     <  1:489907/93‑1 (MQ=255)

     tGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGtc                                                                                   >  2:4025176/1‑92 (MQ=255)

        tttCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGc                                                                                >  1:3334437/1‑93 (MQ=255)

        tttCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGc                                                                                >  1:3334438/1‑93 (MQ=255)

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

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

            ccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACAGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCacc                                                                            >  2:1173856/1‑91 (MQ=255)

             cTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACCCCGGCACGCCAACa                                                                           >  2:3649933/1‑93 (MQ=255)

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

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

               gCAACCTTGTAGCAACACGACACCTCAATGGCGATTGTTCAGGCACGGCGAATATACGCTGCAAAAAGACATCTACTCCGGTACGCCAACAtt                                                                         >  2:2529834/1‑93 (MQ=255)

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

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

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

                  aCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgat                                                                      <  2:3039024/93‑1 (MQ=255)

                  aCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgat                                                                      <  1:4490363/93‑1 (MQ=255)

                   ccTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACCCCCACTCCGGTACGCCAACATtgatg                                                                     <  1:106260/93‑1 (MQ=255)

                    cTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgatga                                                                    >  1:2362121/1‑93 (MQ=255)

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

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

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

                       gTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGAttt                                                                 <  2:1131013/93‑1 (MQ=255)

                        tAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtc                                                                >  2:903/1‑93 (MQ=255)

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

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

                          gCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGCACGCCAACATTGATGATTtcca                                                              >  1:1577000/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                   cgccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcg                                                     <  1:1452936/93‑1 (MQ=255)

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

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

                                    gccTCAATGGCGATTGTTCAGGGACGGAGAATATACGCTGCAAATAGACACCCACTCCGGGACGCCAACATTGATGATTTACACCCAGAAcgc                                                    <  2:3452889/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                        cAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACgccgcc                                                >  2:54369/1‑93 (MQ=255)

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

                                             gCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAcc                                           >  1:4511103/1‑93 (MQ=255)

                                               gATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCgg                                         <  2:659196/93‑1 (MQ=255)

                                               gATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCgg                                         >  1:2136345/1‑93 (MQ=255)

                                                 ttGTTCAGGGATGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTa                                       <  2:2556185/93‑1 (MQ=255)

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

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

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

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

                                                  tGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAg                                      <  2:2487006/93‑1 (MQ=255)

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

                                                      cAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAg                                  <  2:4261728/93‑1 (MQ=255)

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

                                                        gggACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGcc                                <  2:2911496/93‑1 (MQ=255)

                                                        gggACGGCGAATATACGCTGCAAATAGACACCCCCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGcc                                >  1:4570667/1‑93 (MQ=255)

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

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

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

                                                           aCGGCGAATATACGCTGCAAATAGACGCCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTgg                             <  1:70401/93‑1 (MQ=255)

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

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

                                                              gCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTcc                          <  1:3700978/93‑1 (MQ=255)

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

                                                               cGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCg                         <  2:1447589/93‑1 (MQ=255)

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

                                                                 aaTATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTg                       <  2:3334510/93‑1 (MQ=255)

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

                                                                  atatACGCTGCAAATAGCCACCCACAACGGAACGCCAACCTTGATGATTTCCATCCAGAAAGCCGCCGAACCGGTAGCAAGACTAGACCgaga                      >  2:2876796/1‑93 (MQ=255)

                                                                  atatACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGa                      <  1:3137906/93‑1 (MQ=255)

                                                                  atatACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGa                      <  2:4305210/93‑1 (MQ=255)

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

                                                                   tataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCAGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     >  1:678802/1‑93 (MQ=255)

                                                                   tataCGCTGAAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     <  2:370434/93‑1 (MQ=255)

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

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

                                                                        gCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGccc                >  1:4013705/1‑93 (MQ=255)

                                                                        gCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGccc                <  2:860958/93‑1 (MQ=255)

                                                                          tGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGGCCGa              >  2:3002639/1‑93 (MQ=255)

                                                                          tGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGa              >  1:1805915/1‑93 (MQ=255)

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

                                                                             aaaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGaaag           >  1:905159/1‑92 (MQ=255)

                                                                             aaaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAAt           >  2:3364542/1‑93 (MQ=255)

                                                                             aaaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAAt           <  2:2201028/93‑1 (MQ=255)

                                                                             aaaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAAt           <  2:1340493/93‑1 (MQ=255)

                                                                             aaaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAAt           >  1:2742967/1‑93 (MQ=255)

                                                                              aaTAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGATCCGGTACCAAGCCTGGTCCGTGAATGCCCGAAATg          <  2:732031/93‑1 (MQ=255)

                                                                               aTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATgg         >  2:297314/1‑93 (MQ=255)

                                                                                tAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCATGAATGCCCGAAATggg        >  1:1309173/1‑93 (MQ=255)

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

                                                                                    cacCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGCACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATTggatgg    >  2:182327/1‑93 (MQ=255)

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

                                                                                      cccACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatggat  >  1:1692745/1‑93 (MQ=255)

                                                                                      cccACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatggat  >  2:1078581/1‑93 (MQ=255)

                                                                                      cccACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatggat  >  2:1680839/1‑93 (MQ=255)

TCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735326‑3735504

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