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.6%	825.9 / NA	242	R97H (CGC→CAC)	malS	alpha‑amylase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (113/128); minor base C (0/1); total (113/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 252 total aligned reads displayed.

ATCAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCGCTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGGGATGGAT  >  NC_000913/3735325‑3735504

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

  cAGATGCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTcc                                                                                       <  1:345336/93‑1 (MQ=255)

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

   aGATGCTTTCCCTCCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCg                                                                                      <  1:2763336/93‑1 (MQ=255)

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

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

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

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

       gCTTTCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCGAATAGACACCCACTCCGGTAc                                                                                  >  1:1387829/1‑93 (MQ=255)

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

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

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

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

          ttCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGcc                                                                               <  1:5171014/93‑1 (MQ=255)

           tCCCTGCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCa                                                                              <  1:1406105/93‑1 (MQ=255)

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

               tGCAGCCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAt                                                                          <  2:2663774/93‑1 (MQ=255)

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

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

                gCAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACAtt                                                                         >  1:2863488/1‑93 (MQ=255)

                 cAACCTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtg                                                                        <  1:5695501/93‑1 (MQ=255)

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

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

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

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

                    ccTTGTAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATtgatg                                                                     >  1:3502360/1‑93 (MQ=255)

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

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

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

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

                         tAGCAACACGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTtc                                                                <  1:2316955/93‑1 (MQ=255)

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

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

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

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

                              acacGCCGCCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACCCCCACTCCGGTACGCCAACATTGATGATTtccatcc                                                           <  1:1738408/93‑1 (MQ=255)

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

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

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

                                 cgccgcCTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGa                                                        <  2:2537889/93‑1 (MQ=255)

                                   ccgccTCAATGGCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAc                                                      >  2:5737160/1‑93 (MQ=255)

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

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

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

                                      ccTCAATGGCGATTGTGCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAAcgcc                                                   <  1:211776/93‑1 (MQ=255)

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

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

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

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

                                              gCGATTGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAAcc                                           >  2:5346840/1‑93 (MQ=255)

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

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

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

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

                                                   tGTTCAGGGACGGCGAATATACGCTGCAAATAGACACCCACTCCGGTAGGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAg                                      <  1:1885887/93‑1 (MQ=255)

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

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

                                                      tCAGGGACGGCGAATATACGCTGTAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCaa                                   <  1:1601011/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                     ataCGCTGCAAATAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAAt                    <  1:1457171/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                              aaaTAGACACCCACTCCGGTACGCCAACATTGATGATTCCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAAt           <  1:2682676/93‑1 (MQ=255)

                                                                               aaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATg          <  2:2178780/93‑1 (MQ=255)

                                                                               aaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATg          <  2:4155493/93‑1 (MQ=255)

                                                                               aaTAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATg          <  1:1560488/93‑1 (MQ=255)

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

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

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

                                                                                 tAGACACCCACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATggg        <  2:6292081/93‑1 (MQ=255)

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

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

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

                                                                                       cccACTCCGGTACGCCAACATTGATGATTTCCATCCAGAACGCCGCCGAACCGGTAGCAAGCCTGGTCCGTGAATGCCCGAAATGggatggat  >  2:3894060/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` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`