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.8%	1412.1 / NA	432	R97H (CGC→CAC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base G (0/1); major base A (210/216); minor base C (0/5); total (210/222)
Fisher's exact test for biased strand distribution p-value = 6.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 456 total aligned reads displayed.

ATCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735325‑3735504

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

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

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

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

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

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

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

       gCTTGCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTAc                                                                                  <  2:251639/93‑1 (MQ=255)

           tCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCa                                                                              >  1:5953861/1‑93 (MQ=255)

            cccTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCaa                                                                             >  2:5851688/1‑93 (MQ=255)

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

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

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

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

                  aaCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtga                                                                       >  2:7298749/1‑93 (MQ=255)

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

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

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

                    ccTTGTAGCAACACGCCGCCTCAATCGCGTTTAGTCAGGGACCGCGAATATACGCTGCAACTTGACACCCACTCCGGTGCGCGAACACtgatg                                                                     >  2:5656227/1‑93 (MQ=255)

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

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

                       tGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGAtt                                                                  >  1:3682898/1‑93 (MQ=255)

                       tGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGAtt                                                                  <  1:493222/93‑1 (MQ=255)

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

                        gTAGCAACACGCCGCCTCAATGGAGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGAttt                                                                 <  2:5557508/93‑1 (MQ=255)

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

                            cAACACGCCGCCTCAATGGCGATTGTTCAGTGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccat                                                             <  2:1772028/93‑1 (MQ=255)

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

                              acacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatcc                                                           <  1:8239933/93‑1 (MQ=255)

                              acacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatcc                                                           <  1:9783537/93‑1 (MQ=255)

                               cacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtccatcca                                                          <  2:7712999/93‑1 (MQ=255)

                                acGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAg                                                         >  1:3258444/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                           aTGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGa                                              >  2:6516009/1‑93 (MQ=255)

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

                                             ggCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAc                                            >  1:872442/1‑93 (MQ=255)

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

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

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

                                                gATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCAGAACCgg                                         >  2:2301827/1‑93 (MQ=255)

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

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

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

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

                                                    gTTCAGGGACGGCGAATATACGCTGCAAATAGAAACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGc                                     >  2:3314192/1‑93 (MQ=255)

                                                         gggACGGCGAATATACGCTGCAAATAGACGCCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGcc                                <  2:9459372/93‑1 (MQ=255)

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

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

                                                          ggACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGACt                               >  2:1045210/1‑93 (MQ=255)

                                                           gACCGCGAATATACGCTGCAAATAGACCCCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTg                              <  2:3722020/93‑1 (MQ=255)

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

                                                              ggCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTc                           <  2:3113612/93‑1 (MQ=255)

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

                                                              ggCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTc                           <  1:2142180/93‑1 (MQ=255)

                                                              ggCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTCCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTc                           <  1:5766083/93‑1 (MQ=255)

                                                              ggCGAATATACGCTGCAAATAAACACCCACTCCGGGACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGGAGCAAGCCTGGTc                           >  2:212073/1‑93 (MQ=255)

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

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

                                                                cGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCg                         >  2:383073/1‑93 (MQ=255)

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

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

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

                                                                    tataCGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGaa                     <  1:3452854/93‑1 (MQ=255)

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

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

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

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

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

                                                                      taCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAAtc                   >  1:4146861/1‑92 (MQ=255)

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

                                                                          cTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCg               <  1:6090651/93‑1 (MQ=255)

                                                                           tGCAAATAGACGCCCCCTCCGGTACGCCAGCATGGATGATTTCCATCCAGAACGCCGCCGAACAGGTAGCAAGCCTGGTCCGTGAATGCCCGa              <  2:1220747/93‑1 (MQ=255)

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

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

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

                                                                             cAAATAAGCACCCCCTCCGGTACCGCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGaaa            <  1:1878137/93‑1 (MQ=255)

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

                                                                               aaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATg          >  2:2005129/1‑93 (MQ=255)

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

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

                                                                                  aGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAATCCTGGTCCGTGAATGCCCGAAATGgga       <  2:6360914/93‑1 (MQ=255)

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

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

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

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

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

                                                                                     cacCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatgg    >  2:7865553/1‑93 (MQ=255)

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

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

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

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

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

                                                                                       cccACTCCGGTACGACAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatggat  <  1:2396461/93‑1 (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: `‑`