BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,575,878	A→G	Y344H (TAT→CAT)	mcrB ←	McrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,575,878	0	A	G	99.4%	1565.3 / ‑4.1	478	Y344H (TAT→CAT)	mcrB	McrB
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (200/275); total (200/278)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 494 total aligned reads displayed.

CTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTT  >  NC_000913/4575788‑4575967

cTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGt                                                                                          <  2:7252938/93‑1 (MQ=255)

      ttttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTcaacaa                                                                                    >  1:4352787/1‑93 (MQ=255)

        ttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACg                                                                                  <  1:9602858/93‑1 (MQ=255)

        ttATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACg                                                                                  <  1:5752737/93‑1 (MQ=255)

          aTTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGc                                                                                <  2:2877387/93‑1 (MQ=255)

           ttCAGTAAAAAATGCCGGAACTGTGGTGTATCAAACGCTGGCTCAGTATCTATGAAAGAAAATCGTCTGCGTAGAACATGGTCAACAACGGcc                                                                               <  2:1215952/93‑1 (MQ=255)

            tCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCa                                                                              <  1:7154293/93‑1 (MQ=255)

            tCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCa                                                                              <  1:5251563/93‑1 (MQ=255)

             cAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCag                                                                             >  1:1036743/1‑93 (MQ=255)

              aGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCaga                                                                            <  2:7407781/93‑1 (MQ=255)

              aGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCaga                                                                            <  2:3600935/93‑1 (MQ=255)

              aGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCaga                                                                            >  2:1688589/1‑93 (MQ=255)

              aGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCaga                                                                            <  1:7271083/93‑1 (MQ=255)

               gTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagag                                                                           >  2:9687385/1‑93 (MQ=255)

               gTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagag                                                                           <  1:5282384/93‑1 (MQ=255)

                tAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagaga                                                                          >  1:4062208/1‑93 (MQ=255)

                 aaaaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCagagag                                                                         <  1:8857373/93‑1 (MQ=255)

                  aaaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCCACAACGGCCAGAGAGc                                                                        >  2:2032471/1‑93 (MQ=255)

                   aaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCg                                                                       >  2:9500246/1‑93 (MQ=255)

                   aaaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCg                                                                       >  2:2587362/1‑93 (MQ=255)

                    aaaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGa                                                                      >  1:8140437/1‑93 (MQ=255)

                     aaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGAt                                                                     >  1:8023587/1‑93 (MQ=255)

                     aaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGAt                                                                     <  1:4128796/93‑1 (MQ=255)

                     aaTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGAt                                                                     <  2:6672399/93‑1 (MQ=255)

                      aTTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATc                                                                    <  2:1921375/93‑1 (MQ=255)

                       ttCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCg                                                                   <  1:6888672/93‑1 (MQ=255)

                       ttCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCg                                                                   >  1:5807102/1‑93 (MQ=255)

                        tCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGTGATCgg                                                                  <  1:6618943/93‑1 (MQ=255)

                         ccGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGc                                                                 >  1:9169974/1‑93 (MQ=255)

                          cGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCa                                                                >  1:10006512/1‑93 (MQ=255)

                          cGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCa                                                                <  1:5354431/93‑1 (MQ=255)

                          cGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCa                                                                <  1:238002/93‑1 (MQ=255)

                           ggAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAg                                                               >  1:2091902/1‑93 (MQ=255)

                            gAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGt                                                              >  2:6438184/1‑93 (MQ=255)

                            gAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGt                                                              >  1:8173565/1‑93 (MQ=255)

                               cTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTAtt                                                           <  1:120627/93‑1 (MQ=255)

                                tgtgGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTc                                                          >  2:4269399/1‑93 (MQ=255)

                                 gtggtgTATCAAACCCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCa                                                         <  1:6340631/93‑1 (MQ=255)

                                 gtggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCa                                                         >  2:6168859/1‑93 (MQ=255)

                                  tggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑t                                                       <  2:2386972/93‑1 (MQ=255)

                                  tggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑t                                                       <  1:7902718/93‑1 (MQ=255)

                                   ggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGTAGTATTCA‑tt                                                      <  2:9333655/93‑1 (MQ=255)

                                   ggtgTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑tt                                                      <  1:8477766/93‑1 (MQ=255)

                                    gtgtATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTa                                                     <  2:5127734/93‑1 (MQ=255)

                                    gtgtATCAAAACCTGGCTCAATATCTATGAAAAAAAAACGTCTGCGTAGGGGATGGGCAACAACGGCCAGAGAGAGATCGGCAGTATTCA‑TTa                                                     >  2:6769939/1‑93 (MQ=255)

                                     tgtATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTaa                                                    >  2:1069786/1‑93 (MQ=255)

                                      gtATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTaaa                                                   <  2:1434028/93‑1 (MQ=255)

                                       tATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAAc                                                  <  1:4717249/93‑1 (MQ=255)

                                       tATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAAc                                                  <  1:4577271/93‑1 (MQ=255)

                                       tATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAAc                                                  >  2:135514/1‑93 (MQ=255)

                                        aTCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAAcc                                                 <  2:417353/93‑1 (MQ=255)

                                         tCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCg                                                >  1:3964229/1‑93 (MQ=255)

                                         tCAAAACCTGGCTCAATATCTATGAAAAAAAAAAGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGGATTCATTTAAAcc                                                 >  2:3081284/1‑93 (MQ=255)

                                           aaaaCCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCgat                                              >  1:8656981/1‑93 (MQ=255)

                                            aaaCCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCgatg                                             <  2:5725471/93‑1 (MQ=255)

                                            aaaCCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCAATCGGCAGTATTCA‑TTAAACCgatg                                             <  1:5660561/93‑1 (MQ=255)

                                             aaCCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCgatga                                            <  1:187457/93‑1 (MQ=255)

                                                 tGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGAtata                                        <  1:9511509/93‑1 (MQ=255)

                                                 tGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGAtata                                        <  1:6233909/93‑1 (MQ=255)

                                                 tGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGAtata                                        <  2:8644221/93‑1 (MQ=255)

                                                 tGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGAtata                                        >  2:8881909/1‑93 (MQ=255)

                                                 tGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGAtata                                        <  1:4786621/93‑1 (MQ=255)

                                                  ggCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATaa                                       >  1:1696289/1‑93 (MQ=255)

                                                   gCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATaaa                                      <  1:763658/93‑1 (MQ=255)

                                                        atatCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATtc                                 <  1:2460711/93‑1 (MQ=255)

                                                        atatCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATtc                                 <  2:802635/93‑1 (MQ=255)

                                                          atCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATtctc                               >  1:1120379/1‑93 (MQ=255)

                                                          atCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATtctc                               <  2:8090840/93‑1 (MQ=255)

                                                          atCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATtctc                               <  1:2011218/93‑1 (MQ=255)

                                                           tCTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTcc                              >  2:8274142/1‑93 (MQ=255)

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCg                             <  2:7558269/93‑1 (MQ=255)

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCg                             >  1:1479334/1‑93 (MQ=255)

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCg                             <  1:934502/93‑1 (MQ=255)

                                                            cTATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCg                             <  1:5953445/93‑1 (MQ=255)

                                                             tATGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCgg                            >  2:3826353/1‑93 (MQ=255)

                                                               tGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGa                          <  2:4467278/93‑1 (MQ=255)

                                                               tGAAAGAAAATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGa                          >  2:4920872/1‑93 (MQ=255)

                                                                gaaagaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGAc                         <  2:9033255/93‑1 (MQ=255)

                                                                  aagaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACAt                       >  1:2695093/1‑93 (MQ=255)

                                                                    gaaaATCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAg                     >  1:3467455/1‑93 (MQ=255)

                                                                     aaaaTCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGa                    <  1:3825362/93‑1 (MQ=255)

                                                                      aaaTCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGaa                   <  1:3700872/93‑1 (MQ=255)

                                                                         tCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCg                <  2:5456287/93‑1 (MQ=255)

                                                                         tCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCg                <  1:4866337/93‑1 (MQ=255)

                                                                         tCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCg                <  1:2859839/93‑1 (MQ=255)

                                                                         tCGTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCg                <  1:5516349/93‑1 (MQ=255)

                                                                           gTCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGtt              <  2:335478/93‑1 (MQ=255)

                                                                            tCTGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttc             <  1:6712257/93‑1 (MQ=255)

                                                                              tGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttctt           >  1:3226759/1‑93 (MQ=255)

                                                                              tGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttctt           <  1:7756548/93‑1 (MQ=255)

                                                                              tGCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttctt           >  2:4706389/1‑93 (MQ=255)

                                                                               gCGTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGttcttc          <  2:7066400/93‑1 (MQ=255)

                                                                                ggTAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCa         <  2:7500782/92‑1 (MQ=255)

                                                                                  tAGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATc       >  1:7006492/1‑93 (MQ=255)

                                                                                   aGGGCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCg      >  1:1322378/1‑93 (MQ=255)

                                                                                     ggCATGGTCAACAACGGCCAGAGAGCGATCGGGAGTATTCA‑TTAAAACGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtt    >  2:595035/1‑93 (MQ=255)

                                                                                     ggCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtt    >  2:6540218/1‑93 (MQ=255)

                                                                                     ggCATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtt    >  2:5324357/1‑93 (MQ=255)

                                                                                       cATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtttt  <  1:2220423/93‑1 (MQ=255)

                                                                                       cATGGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGtttt  <  1:160819/93‑1 (MQ=255)

CTGCTTTTTTATTCAGTAAAAAATTCCGGAACTGTGGTGTATCAAAACCTGGCTCAATATCTATGAAAGAAAATCGTCTGCGTAGGGCATAGTCAACAACGGCCAGAGAGCGATCGGCAGTATTCA‑TTAAACCGATGATATAAACATTCTCCGGGACATAGAATCGTTCTTCATCGTTTT  >  NC_000913/4575788‑4575967

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