BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,372,457	A→G	V260A (GTG→GCG)	yjeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,372,457	0	A	G	99.8%	1676.1 / ‑5.2	516	V260A (GTG→GCG)	yjeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (223/291); total (223/293)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 547 total aligned reads displayed.

ACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCA  >  NC_000913/4372367‑4372546

aCTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTcgca                                                                                         >  1:538794/1‑93 (MQ=255)

aCTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTcgca                                                                                         <  1:8076184/93‑1 (MQ=255)

aCTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTcgca                                                                                         >  1:887986/1‑93 (MQ=255)

  tttATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATc                                                                                       >  1:7276745/1‑93 (MQ=255)

    tATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCAc                                                                                     <  1:3321707/93‑1 (MQ=255)

    tATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCAc                                                                                     <  2:2987942/93‑1 (MQ=255)

      tCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGt                                                                                   <  2:6950943/93‑1 (MQ=255)

       cGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTa                                                                                  <  1:3463189/93‑1 (MQ=255)

         aGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAAc                                                                                <  1:7005171/93‑1 (MQ=255)

            acaTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGa                                                                             <  2:3543742/93‑1 (MQ=255)

             caTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGaa                                                                            <  1:4122718/93‑1 (MQ=255)

             caTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGaa                                                                            <  1:7747686/93‑1 (MQ=255)

             caTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGaa                                                                            <  1:5011938/93‑1 (MQ=255)

             caTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGaa                                                                            <  2:9744025/93‑1 (MQ=255)

             caTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGaa                                                                            <  2:2108286/93‑1 (MQ=255)

              aTGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAAc                                                                           <  1:7501559/93‑1 (MQ=255)

                gCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGc                                                                         <  2:5853745/93‑1 (MQ=255)

                 cAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTTCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGct                                                                        >  2:8951337/1‑93 (MQ=255)

                  aGGTGATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGctc                                                                       <  2:3977363/93‑1 (MQ=255)

                   ggTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGctct                                                                      >  1:2012369/1‑93 (MQ=255)

                   ggTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGctct                                                                      >  1:1587746/1‑93 (MQ=255)

                   ggTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGctcg                                                                      >  2:6769268/1‑92 (MQ=255)

                    gTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTg                                                                     >  2:8281164/1‑93 (MQ=255)

                    gTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTg                                                                     <  2:9310712/93‑1 (MQ=255)

                      aaTACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGt                                                                   >  1:8745062/1‑93 (MQ=255)

                      aaTACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGt                                                                   >  1:1178205/1‑93 (MQ=255)

                       aTACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGtt                                                                  >  2:4162915/1‑93 (MQ=255)

                        tACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTc                                                                 <  1:8547957/93‑1 (MQ=255)

                         aCGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCa                                                                <  1:7944964/93‑1 (MQ=255)

                          cGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAg                                                               <  2:5086277/93‑1 (MQ=255)

                            gCAATACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCa                                                             <  2:2693929/93‑1 (MQ=255)

                              aTTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCaaa                                                           <  2:5614324/93‑1 (MQ=255)

                               ttACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAg                                                          >  1:5241789/1‑93 (MQ=255)

                                tccgccgcCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGt                                                         <  1:1348829/91‑1 (MQ=255)

                                 aCGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTa                                                        <  1:5711872/93‑1 (MQ=255)

                                  cGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGATCGCTCTGGTTCAGCAAAGTaa                                                       <  2:6118707/93‑1 (MQ=255)

                                  cGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTaa                                                       <  1:9803170/93‑1 (MQ=255)

                                  cGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTaa                                                       >  2:7397286/1‑93 (MQ=255)

                                   gccgccATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAAc                                                      <  2:8108905/93‑1 (MQ=255)

                                    ccGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACg                                                     >  2:638973/1‑93 (MQ=255)

                                    ccGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACg                                                     >  2:720140/1‑93 (MQ=255)

                                     cGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGc                                                    >  2:5751951/1‑93 (MQ=255)

                                     cGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGc                                                    >  1:425552/1‑93 (MQ=255)

                                      ggCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGcc                                                   >  2:761597/1‑93 (MQ=255)

                                       gCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGCACGCTCTGGTTCAGCAAAGTAACGCCt                                                  >  2:4850266/1‑93 (MQ=255)

                                       gCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCt                                                  >  1:2806292/1‑93 (MQ=255)

                                       gCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCt                                                  <  2:163137/93‑1 (MQ=255)

                                       gCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCt                                                  <  1:7639489/93‑1 (MQ=255)

                                       gCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCt                                                  <  1:6437985/93‑1 (MQ=255)

                                        cATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTa                                                 <  1:5927385/93‑1 (MQ=255)

                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               <  1:8045079/93‑1 (MQ=255)

                                          tCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAcc                                               <  2:1911276/93‑1 (MQ=255)

                                           cGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTAccc                                              <  2:3139047/93‑1 (MQ=255)

                                             aacaacGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCgg                                            >  1:4663046/1‑93 (MQ=255)

                                             aacaacGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAAGGCCTAccccc                                            >  1:9266669/1‑91 (MQ=255)

                                              acaacGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGgc                                           >  1:3167255/1‑93 (MQ=255)

                                               caacGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGgcg                                          >  1:1770389/1‑93 (MQ=255)

                                                 acGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCa                                        >  1:5654047/1‑93 (MQ=255)

                                                   gCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAAc                                      >  1:9021130/1‑93 (MQ=255)

                                                     aTTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACtt                                    <  2:4285586/93‑1 (MQ=255)

                                                      ttACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTa                                   >  2:5507840/1‑93 (MQ=255)

                                                      ttACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTa                                   <  1:6710511/93‑1 (MQ=255)

                                                      ttACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTa                                   >  2:10032646/1‑93 (MQ=255)

                                                      ttACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTa                                   >  2:7794168/1‑93 (MQ=255)

                                                       tACTCAGGTGTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACATAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAg                                  <  2:9867237/93‑1 (MQ=255)

                                                         cTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGcc                                <  2:9183586/93‑1 (MQ=255)

                                                         cTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGcc                                >  1:2399835/1‑93 (MQ=255)

                                                          tCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCa                               >  1:9511158/1‑93 (MQ=255)

                                                          tCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCa                               >  1:7091848/1‑93 (MQ=255)

                                                          tCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCa                               <  2:7693637/93‑1 (MQ=255)

                                                           cAGGTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCAt                              <  1:2498782/93‑1 (MQ=255)

                                                           cAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCAt                              <  1:3779648/93‑1 (MQ=255)

                                                              gTTTGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGc                           >  2:7903059/1‑93 (MQ=255)

                                                               tttGCCAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGcc                          <  1:4046028/93‑1 (MQ=255)

                                                                   ccAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGAc                      >  1:9396834/1‑93 (MQ=255)

                                                                    cAGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACg                     <  1:4554961/93‑1 (MQ=255)

                                                                     aGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACgg                    >  2:654330/1‑93 (MQ=255)

                                                                     aGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACgg                    >  2:1304123/1‑93 (MQ=255)

                                                                     aGCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACgg                    >  1:2133989/1‑93 (MQ=255)

                                                                      gCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGa                   >  1:10004394/1‑93 (MQ=255)

                                                                      gCGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGa                   >  1:3626472/1‑93 (MQ=255)

                                                                       cGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGaa                  <  2:8405526/93‑1 (MQ=255)

                                                                       cGTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGaa                  <  2:9407425/93‑1 (MQ=255)

                                                                        gTTTGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAAt                 >  2:5987255/1‑93 (MQ=255)

                                                                         tttGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATg                >  2:2305547/1‑93 (MQ=255)

                                                                          ttGTGCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGt               <  2:3400750/93‑1 (MQ=255)

                                                                           tgggCGGTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGtt              >  2:1523491/4‑93 (MQ=255)

                                                                            gtgCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGttt             >  1:1807879/1‑93 (MQ=255)

                                                                              gCGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCa           <  2:7548662/93‑1 (MQ=255)

                                                                               cGTTATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCAt          <  1:2654904/93‑1 (MQ=255)

                                                                                 ttATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCt        <  1:7348444/93‑1 (MQ=255)

                                                                                 ttATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCt        >  2:3771787/1‑93 (MQ=255)

                                                                                 ttATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCt        >  2:8628658/1‑93 (MQ=255)

                                                                                  tATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTa       <  1:4903679/93‑1 (MQ=255)

                                                                                  tATCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTa       <  1:2984434/93‑1 (MQ=255)

                                                                                    tCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAcc     <  1:6873482/93‑1 (MQ=255)

                                                                                    tCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAcc     <  2:4493118/93‑1 (MQ=255)

                                                                                    tCGTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTAcc     <  2:1125694/93‑1 (MQ=255)

                                                                                      gTTCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACATc   >  1:613359/1‑93 (MQ=255)

                                                                                       ttCGCATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCa  <  1:1129413/93‑1 (MQ=255)

ACTTTATCGAGCACATGCAGGTAATACGGCATTACGCCGGCATCGAACAACGCATTACTCAGGTTTGCCAGCGTTTGTGCGTTATCGTTCACATCACGTAACAGAACGCTCTGGTTCAGCAAAGTAACGCCTACCCGGCGCAACTTAGCCATCGCCTGACGGAATGTTTCATCTACCTCA  >  NC_000913/4372367‑4372546

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