BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,147,941	A→G	V356A (GTC→GCC)	b2070 ←	putative chaperonin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,941	0	A	G	89.0%	680.0 / 56.0	263	V356A (GTC→GCC)	b2070	putative chaperonin
Reads supporting (aligned to +/- strand): ref base A (10/19); new base G (109/125); total (119/144)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 263 total aligned reads displayed.

AAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACC  >  NC_000913/2147851‑2148031

   ggAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTc                                                                                       >  1:4083114/1‑93 (MQ=255)

   ggAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTc                                                                                       <  2:100448/93‑1 (MQ=255)

            cTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGcaca                                                                              >  2:1542853/1‑93 (MQ=255)

             tGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACaa                                                                             >  1:7058464/1‑93 (MQ=255)

             tGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACaa                                                                             <  2:4355254/93‑1 (MQ=255)

              gCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACaaa                                                                            >  1:354227/1‑93 (MQ=255)

                 aGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCGATCACCATCGCCCAGGCTTTCGCCAGCACAAAAtt                                                                         <  2:4526740/93‑1 (MQ=255)

                  gTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTc                                                                        <  2:7366979/93‑1 (MQ=255)

                   tCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTcc                                                                       >  1:1951008/1‑93 (MQ=255)

                   tCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTcc                                                                       >  1:1276393/1‑93 (MQ=255)

                   tCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTcc                                                                       <  2:4206011/93‑1 (MQ=255)

                   tCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTcc                                                                       >  1:5422511/1‑93 (MQ=255)

                        tGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCccgccg                                                                  >  2:555427/1‑93 (MQ=255)

                          gAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggt                                                                <  2:726377/93‑1 (MQ=255)

                          gAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggt                                                                <  2:3283570/93‑1 (MQ=255)

                           aaaTGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggtg                                                               <  2:1961158/93‑1 (MQ=255)

                            aaTGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCGAGCACAAAATTCCCGCCggtgg                                                              <  1:1025558/93‑1 (MQ=255)

                            aaTGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggtgg                                                              >  2:7034311/1‑93 (MQ=255)

                            aaTGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCggtgg                                                              <  2:6040946/93‑1 (MQ=255)

                             aTGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCggtggt                                                             >  2:1902838/1‑93 (MQ=255)

                                aTCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAAt                                                          <  1:5724200/93‑1 (MQ=255)

                                  cGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGa                                                        >  2:2401811/1‑93 (MQ=255)

                                   gggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAc                                                       <  1:2728002/93‑1 (MQ=255)

                                    ggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGagc                                                      >  1:2797519/1‑91 (MQ=255)

                                    ggTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGAcc                                                      <  2:6992361/93‑1 (MQ=255)

                                     gTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACct                                                     <  1:570461/93‑1 (MQ=255)

                                       aaTGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACctct                                                   <  1:2312086/93‑1 (MQ=255)

                                       aaTGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACctct                                                   >  2:6648980/1‑93 (MQ=255)

                                        aTGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCtt                                                  <  1:6911301/93‑1 (MQ=255)

                                        aTGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCtt                                                  <  2:5170349/93‑1 (MQ=255)

                                         tGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTc                                                 <  2:372468/93‑1 (MQ=255)

                                         tGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTc                                                 <  1:1848047/93‑1 (MQ=255)

                                           gCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGc                                               >  2:2881788/1‑93 (MQ=255)

                                             aGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCg                                             >  1:4536538/1‑93 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                            >  1:3215226/1‑93 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                            <  2:5452729/93‑1 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                            >  2:6136378/1‑93 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                            <  2:112730/93‑1 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                            <  2:2221751/93‑1 (MQ=255)

                                              gTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCga                                            >  1:6724703/1‑93 (MQ=255)

                                                tGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCgaga                                          <  1:5798344/93‑1 (MQ=255)

                                                tGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCgaga                                          >  1:4193342/1‑93 (MQ=255)

                                                 gCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGAt                                         <  1:5596821/93‑1 (MQ=255)

                                                 gCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGAt                                         <  2:305435/93‑1 (MQ=255)

                                                 gCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGAt                                         >  1:781992/1‑93 (MQ=255)

                                                   aaTGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCa                                       >  1:2119421/1‑93 (MQ=255)

                                                   aaTGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCa                                       <  1:1438333/93‑1 (MQ=255)

                                                   aaTGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCa                                       >  2:4672188/1‑93 (MQ=255)

                                                    aggCTGGCGGTGCTACTGCCGCCATCCCCATCGCCCAGGCTTTCGCCAGCACAAAATCCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCaa                                      <  1:6531284/91‑1 (MQ=255)

                                                    aTGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCaa                                      <  2:3573342/93‑1 (MQ=255)

                                                     tGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAAt                                     >  1:1633438/1‑93 (MQ=255)

                                                     tGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAAt                                     >  1:2503825/1‑93 (MQ=255)

                                                      gCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATa                                    <  2:4474382/93‑1 (MQ=255)

                                                       cTGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaa                                   <  2:3543156/93‑1 (MQ=255)

                                                        tGGCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaa                                  >  2:4698228/1‑93 (MQ=255)

                                                         ggCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATaaaa                                 >  1:6377356/1‑93 (MQ=255)

                                                          gCGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAAt                                <  1:5360079/93‑1 (MQ=255)

                                                           cGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATa                               <  1:436318/93‑1 (MQ=255)

                                                           cGGTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATa                               <  1:5622840/93‑1 (MQ=255)

                                                            ggTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAg                              >  2:4239536/1‑93 (MQ=255)

                                                            ggTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAg                              <  2:6790660/93‑1 (MQ=255)

                                                             gTGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGa                             <  2:4832658/93‑1 (MQ=255)

                                                              tGCTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGAt                            >  1:3597502/1‑93 (MQ=255)

                                                              tGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGAt                            >  1:4477031/1‑93 (MQ=255)

                                                                cTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATaa                          <  1:1113254/93‑1 (MQ=255)

                                                                cTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATaa                          <  1:5915726/93‑1 (MQ=255)

                                                                cTACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATaa                          <  1:7222837/93‑1 (MQ=255)

                                                                cTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATaa                          <  2:5071676/93‑1 (MQ=255)

                                                                cTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATaa                          <  1:2889846/93‑1 (MQ=255)

                                                                 tACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAac                         >  2:2330729/1‑93 (MQ=255)

                                                                 tACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAac                         <  1:3160253/93‑1 (MQ=255)

                                                                 tACTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAac                         >  1:3081916/1‑93 (MQ=255)

                                                                  aCTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAaca                        >  1:51989/1‑93 (MQ=255)

                                                                   cTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                       <  2:1104214/93‑1 (MQ=255)

                                                                   cTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                       >  1:4765529/1‑93 (MQ=255)

                                                                   cTGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAacac                       >  2:2576375/1‑93 (MQ=255)

                                                                    tGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACAcc                      >  2:6340724/1‑93 (MQ=255)

                                                                    tGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACAcc                      <  1:3995359/93‑1 (MQ=255)

                                                                    tGACGCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACAcc                      <  1:2344477/93‑1 (MQ=255)

                                                                        gCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGa                  >  1:4817638/1‑93 (MQ=255)

                                                                        gCCATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGa                  <  2:2778885/93‑1 (MQ=255)

                                                                        cccATCCCCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGa                  <  1:2559557/92‑1 (MQ=255)

                                                                         ccATCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGAt                 <  2:155892/93‑1 (MQ=255)

                                                                            tCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTAt              >  2:1407534/1‑93 (MQ=255)

                                                                            tCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTAt              >  2:3885053/1‑93 (MQ=255)

                                                                            tCACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCCCCGGGGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTAt              >  1:5089264/1‑93 (MQ=255)

                                                                             cACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATc             <  2:861116/93‑1 (MQ=255)

                                                                             cACCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATc             <  2:3169326/93‑1 (MQ=255)

                                                                              aCCATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATcc            >  2:581748/1‑93 (MQ=255)

                                                                               ccATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCa           <  2:4865861/93‑1 (MQ=255)

                                                                               ccATCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCa           <  1:2196245/93‑1 (MQ=255)

                                                                                 aTCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgc         <  1:4597697/93‑1 (MQ=255)

                                                                                 aTCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgc         <  1:6198758/93‑1 (MQ=255)

                                                                                  tCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcg        <  2:4585453/93‑1 (MQ=255)

                                                                                  tCGCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAgcg        >  2:3965899/1‑93 (MQ=255)

                                                                                    gCCCAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGca      >  2:2115823/1‑93 (MQ=255)

                                                                                      ccAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGcaca    >  2:2673341/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac   >  1:3763547/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac   >  1:3388938/1‑93 (MQ=255)

                                                                                       cAGGCTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCacac   <  2:2747177/93‑1 (MQ=255)

AAAGGAAGCGGCCTGCCAGTCGGTTGGAAATGATCGGGTAATGGCAGTTGCAATGCTGGCGGTGCTACTGACGCCATCACCATCGCCCAGACTTTCGCCAGCACAAAATTCCCGCCGGTGGTAATGACCTCTTCGCCGAGATCAATAAAATAGATAACACCTTGATTATCCAGCGCACACC  >  NC_000913/2147851‑2148031

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