BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,903,375	A→G	V274A (GTG→GCG)	bglG ←	BglG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,903,375	0	A	G	96.1%	2087.1 / 42.9	699	V274A (GTG→GCG)	bglG	BglG
Reads supporting (aligned to +/- strand): ref base A (12/15); new base G (338/332); total (352/347)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 728 total aligned reads displayed.

GACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCC  >  NC_000913/3903285‑3903464

gACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGcgcgcg                                                                                         >  2:1682071/1‑93 (MQ=255)

   gTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCgcgcgcgc                                                                                      <  2:5337991/93‑1 (MQ=255)

      aaTCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGctct                                                                                   <  1:5384380/93‑1 (MQ=255)

      aaTCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGctct                                                                                   <  1:4486435/93‑1 (MQ=255)

       aTCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCta                                                                                  >  1:501758/1‑93 (MQ=255)

       aTCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCta                                                                                  >  2:404027/1‑93 (MQ=255)

         cAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCtata                                                                                <  1:3463987/93‑1 (MQ=255)

         cAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGAGCGCGCTCtata                                                                                >  1:8176573/1‑93 (MQ=255)

          aaGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTatat                                                                               <  2:2104330/93‑1 (MQ=255)

          aaGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTatat                                                                               >  1:3140058/1‑93 (MQ=255)

          aaGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTatat                                                                               >  1:970787/1‑93 (MQ=255)

             cAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTa                                                                            <  2:600400/93‑1 (MQ=255)

             cAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTa                                                                            <  2:7350610/93‑1 (MQ=255)

             cAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTa                                                                            >  1:8703489/1‑93 (MQ=255)

             cAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTa                                                                            <  2:9429640/93‑1 (MQ=255)

              aaTTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTAt                                                                           >  1:221340/1‑93 (MQ=255)

              aaTTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTAt                                                                           >  2:3787725/1‑93 (MQ=255)

               aTTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATg                                                                          >  1:6788887/1‑93 (MQ=255)

               aTTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATg                                                                          <  2:4517132/93‑1 (MQ=255)

                 ttttttCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGc                                                                        >  2:9919306/1‑93 (MQ=255)

                     ttCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTaaa                                                                    <  2:9278490/93‑1 (MQ=255)

                     ttCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTaaa                                                                    <  1:2003429/93‑1 (MQ=255)

                        aGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAAc                                                                 >  2:4970999/1‑93 (MQ=255)

                         ggTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACa                                                                >  2:6362599/1‑93 (MQ=255)

                         ggTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACa                                                                <  2:6238294/93‑1 (MQ=255)

                          gTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACAt                                                               <  2:5148288/93‑1 (MQ=255)

                           ttttGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATa                                                              <  1:9114147/93‑1 (MQ=255)

                           ttttGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATa                                                              <  1:6057232/93‑1 (MQ=255)

                            tttGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATaa                                                             <  2:745089/93‑1 (MQ=255)

                            tttGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATaa                                                             <  2:7139026/93‑1 (MQ=255)

                            tttGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATaa                                                             <  1:8788871/93‑1 (MQ=255)

                                cccGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAAtctc                                                         <  1:1544499/93‑1 (MQ=255)

                                cccGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAAtctc                                                         >  1:7905115/1‑93 (MQ=255)

                                cccGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAAtctc                                                         <  2:8258895/93‑1 (MQ=255)

                                 acGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATctct                                                        >  1:4838376/2‑93 (MQ=255)

                                  cGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTg                                                       >  2:2257585/1‑93 (MQ=255)

                                   gCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGc                                                      <  2:2842159/93‑1 (MQ=255)

                                     ttAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCgg                                                    <  2:7657955/93‑1 (MQ=255)

                                     ttAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCgg                                                    >  1:9280986/1‑93 (MQ=255)

                                       aGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCggga                                                  >  1:8441362/1‑92 (MQ=255)

                                        gTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTg                                                 >  2:9151954/1‑93 (MQ=255)

                                         tGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGa                                                >  2:109283/1‑93 (MQ=255)

                                          gCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGaa                                               >  2:2389415/1‑93 (MQ=255)

                                            ggTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAAt                                             >  2:8030639/1‑93 (MQ=255)

                                            ggTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAAt                                             >  1:6354102/1‑93 (MQ=255)

                                             gTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAAtt                                            >  2:9041525/1‑93 (MQ=255)

                                             gTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAAtt                                            <  1:3224458/93‑1 (MQ=255)

                                             gTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAAtt                                            >  1:7615709/1‑93 (MQ=255)

                                                aCAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTa                                         >  2:6634657/1‑93 (MQ=255)

                                                aCAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTa                                         <  1:4124699/93‑1 (MQ=255)

                                                 cAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTAc                                        <  2:9782102/93‑1 (MQ=255)

                                                 cAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTAc                                        >  2:3392532/1‑93 (MQ=255)

                                                   aTCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGt                                      >  2:8566788/1‑93 (MQ=255)

                                                   aTCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGt                                      <  2:5510270/93‑1 (MQ=255)

                                                   aTCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGt                                      <  1:7007186/93‑1 (MQ=255)

                                                   aTCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGt                                      >  2:1955512/1‑93 (MQ=255)

                                                   aTCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGt                                      <  2:8332601/93‑1 (MQ=255)

                                                    tCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGtt                                     >  2:3252951/1‑93 (MQ=255)

                                                     ccTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTg                                    >  1:8540065/1‑93 (MQ=255)

                                                      cTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGa                                   >  1:7150651/1‑93 (MQ=255)

                                                      cTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGa                                   >  2:3923537/1‑93 (MQ=255)

                                                       tttACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTNCGGGTGAAATTTTACGGTGAt                                  >  1:6239887/1‑93 (MQ=255)

                                                       tttACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGAt                                  <  2:4901723/93‑1 (MQ=255)

                                                         tACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATAc                                >  2:2657032/1‑93 (MQ=255)

                                                         tACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATAc                                >  2:9581351/1‑93 (MQ=255)

                                                         tACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATAc                                <  2:4617568/93‑1 (MQ=255)

                                                         tACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATAc                                >  1:5881336/1‑93 (MQ=255)

                                                         tACTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATAc                                <  2:4253488/93‑1 (MQ=255)

                                                          aCTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACt                               >  1:1211694/1‑93 (MQ=255)

                                                          aCTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACt                               <  2:1855077/93‑1 (MQ=255)

                                                          aCTCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACt                               >  1:6616736/1‑93 (MQ=255)

                                                            tCAGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGc                             >  1:3607026/1‑93 (MQ=255)

                                                             cAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCa                            <  1:2351258/93‑1 (MQ=255)

                                                              aGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCaa                           >  2:303045/1‑93 (MQ=255)

                                                              aGTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCaa                           <  2:7780853/93‑1 (MQ=255)

                                                              aGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGCTACTGCaa                           >  1:5210800/1‑93 (MQ=255)

                                                               gTAATAATATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCaaa                          >  2:162782/1‑93 (MQ=255)

                                                                 aataatATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAAcc                        >  2:6816082/1‑93 (MQ=255)

                                                                 aataatATTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAAcc                        >  2:6989841/1‑93 (MQ=255)

                                                                   taataTTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCaa                      >  1:9897861/1‑93 (MQ=255)

                                                                   taataTTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCaa                      <  1:2945104/93‑1 (MQ=255)

                                                                   taataTTTCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCaa                      >  1:8384928/1‑93 (MQ=255)

                                                                        tttCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATaaaa                 >  1:1704404/1‑93 (MQ=255)

                                                                         ttCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATaaaaa                <  2:5653379/93‑1 (MQ=255)

                                                                         ttCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATaaaaa                <  2:3491990/93‑1 (MQ=255)

                                                                          tCAGTGTTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAAt               <  1:3878851/93‑1 (MQ=255)

                                                                          tCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAAt               <  1:474556/93‑1 (MQ=255)

                                                                             gtgtTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGc            <  1:5698573/93‑1 (MQ=255)

                                                                               gtgcTTGGCGCGCGCGTTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTGTACGTTGATACTGCAAACCATTAAAAAGGGGGa          <  2:1499380/90‑1 (MQ=255)

                                                                               gtTCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGa          >  1:5580223/1‑93 (MQ=255)

                                                                                ttCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGAt         <  2:7872645/93‑1 (MQ=255)

                                                                                 tCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATc        >  1:4234099/1‑93 (MQ=255)

                                                                                 tCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATc        >  1:7373779/1‑93 (MQ=255)

                                                                                 tCTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATc        <  1:3320146/93‑1 (MQ=255)

                                                                                  cTTTGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATcc       <  2:8903340/93‑1 (MQ=255)

                                                                                   tttGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCg      <  2:1293142/93‑1 (MQ=255)

                                                                                    ttGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGt     <  1:2654040/93‑1 (MQ=255)

                                                                                    tggcgcgcgcGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGCGTAATTTTACGGTGATACTGCAAACCAATAATAATGGCGATCCGt     >  1:8935736/3‑93 (MQ=255)

                                                                                     tGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGtt    >  1:9586401/1‑93 (MQ=255)

                                                                                     tGCGCGCGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGtt    <  1:1008961/93‑1 (MQ=255)

GACGTGAATCAAGCAATTTTTTTCAGGTTTTGCCCGCTTAGTGCGGTAACAATCCTTTACTCAGTAATAATATTTCAGTGTTCTTTGCGCACGCGCTCTATATTTATGGCTAAAAACATAATCTCTGCGGGTGAAATTTTACGTTGATACTGCAAACCAATAAAAATGGCGATCCGTTCC  >  NC_000913/3903285‑3903464

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: `‑`