BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,006,483	+C	coding (183/207 nt)	b4554 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,006,480	1	.	C	98.5%	726.6 / 1.5	204	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (114/87); total (117/87)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 221 total aligned reads displayed.

GCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006389‑4006567

gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                         <  2:2302841/93‑1 (MQ=255)

gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                         <  2:5356414/93‑1 (MQ=255)

  gcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccc                                                                                       <  2:3304621/93‑1 (MQ=255)

     ccACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaa                                                                                    <  2:4022865/93‑4 (MQ=255)

      cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                   >  2:5689633/1‑89 (MQ=255)

      cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                   >  1:73696/1‑89 (MQ=255)

      cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACGTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                   >  1:4560178/1‑89 (MQ=255)

        ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AAccc                                                                                 <  1:4935043/93‑1 (MQ=255)

        ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AAccc                                                                                 <  1:839972/93‑1 (MQ=255)

        ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑CAaccc                                                                                >  2:2328007/1‑93 (MQ=255)

        ctgctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑CAaccc                                                                                >  2:1477029/1‑93 (MQ=255)

          gctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCa                                                                               >  2:5576641/1‑93 (MQ=255)

          gctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCa                                                                               <  1:4309714/93‑1 (MQ=255)

          gctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCa                                                                               >  1:3183534/1‑93 (MQ=255)

             ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGc                                                                            >  1:4428693/1‑93 (MQ=255)

             ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGc                                                                            >  2:254386/1‑93 (MQ=255)

               aTTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCa                                                                          <  2:111957/93‑1 (MQ=255)

                ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCaa                                                                         >  1:2760036/1‑93 (MQ=255)

                   ccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCa                                                                      <  2:537317/93‑1 (MQ=255)

                     aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGc                                                                    >  1:5118161/1‑93 (MQ=255)

                       ccAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AAACCCACGCCAACCAGcaa                                                                  >  1:3353793/1‑92 (MQ=255)

                        cAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGcaca                                                                 <  1:4961092/93‑1 (MQ=255)

                         aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACaa                                                                >  2:1059472/1‑93 (MQ=255)

                          aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAat                                                               >  2:4457515/1‑93 (MQ=255)

                            tCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAatat                                                             >  1:3540198/1‑93 (MQ=255)

                              aaGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGa                                                           <  1:2426974/93‑1 (MQ=255)

                               aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGaa                                                          >  1:2177621/1‑93 (MQ=255)

                                 tAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATc                                                        <  2:4881658/93‑1 (MQ=255)

                                   aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGTGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGc                                                      <  2:5656140/93‑1 (MQ=255)

                                   aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  2:305532/1‑93 (MQ=255)

                                    gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  1:5450399/93‑1 (MQ=255)

                                     ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  1:4057637/93‑1 (MQ=255)

                                       gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  2:4712601/1‑93 (MQ=255)

                                       gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  1:1394375/1‑93 (MQ=255)

                                        cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  2:1994415/93‑1 (MQ=255)

                                         gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                <  1:5068422/93‑1 (MQ=255)

                                         gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  1:990656/1‑93 (MQ=255)

                                         gcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTcc                                                >  2:3585883/1‑93 (MQ=255)

                                          ccTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCg                                               >  2:4219120/1‑93 (MQ=255)

                                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  2:3451004/1‑93 (MQ=255)

                                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:6056973/1‑93 (MQ=255)

                                           cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCTAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:2279640/1‑93 (MQ=255)

                                            tGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGt                                             >  1:5515440/1‑93 (MQ=255)

                                             gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            >  1:2641335/1‑93 (MQ=255)

                                             gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            <  1:2118767/93‑1 (MQ=255)

                                             gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            >  1:4635155/1‑93 (MQ=255)

                                              cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  2:2941972/1‑93 (MQ=255)

                                              cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           <  2:3894208/93‑1 (MQ=255)

                                              cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  2:3074240/1‑93 (MQ=255)

                                               ccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCta                                          >  1:1661383/1‑92 (MQ=255)

                                                cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  2:5226547/1‑93 (MQ=255)

                                                cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         <  2:762238/93‑1 (MQ=255)

                                                 acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  2:427191/93‑1 (MQ=255)

                                                 acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:1647996/1‑93 (MQ=255)

                                                 acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  1:1501612/1‑93 (MQ=255)

                                                 acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:3248584/1‑93 (MQ=255)

                                                  cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       >  2:1709788/1‑93 (MQ=255)

                                                   ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      >  2:2608016/1‑93 (MQ=255)

                                                   gtcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      <  2:2813882/92‑1 (MQ=255)

                                                    tcaacacatcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  1:563667/1‑93 (MQ=255)

                                                     caacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                                    <  2:4592599/93‑1 (MQ=255)

                                                     caacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGccc                                    <  2:5725716/93‑1 (MQ=255)

                                                      aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATTAATCGCCAGTCCGGTCTGAAGCCCa                                   <  1:4223205/93‑1 (MQ=255)

                                                      aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  2:5289873/93‑1 (MQ=255)

                                                      aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   >  2:3786557/1‑93 (MQ=255)

                                                       acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  >  2:115953/1‑93 (MQ=255)

                                                       acacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAg                                  >  2:5872681/1‑93 (MQ=255)

                                                         acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  1:1585826/1‑93 (MQ=255)

                                                           ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              <  1:1222751/93‑1 (MQ=255)

                                                            tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  1:1064190/93‑1 (MQ=255)

                                                            tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  2:3952202/93‑1 (MQ=255)

                                                              aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  2:5410082/93‑1 (MQ=255)

                                                               aaCGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          <  2:3676386/93‑1 (MQ=255)

                                                                 cGCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGACCAGCAATAGACGc                        >  1:4764575/1‑93 (MQ=255)

                                                                  gCAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  1:5309271/1‑93 (MQ=255)

                                                                   cAATCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCa                      >  1:3521385/1‑93 (MQ=255)

                                                                    aaTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:3869009/1‑93 (MQ=255)

                                                                    aaTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     <  1:5550612/93‑1 (MQ=255)

                                                                    aaTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAc                     >  1:1183451/1‑93 (MQ=255)

                                                                     aTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    <  2:3875728/93‑1 (MQ=255)

                                                                     aTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  1:2449149/1‑93 (MQ=255)

                                                                     aTCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    >  2:6104170/1‑93 (MQ=255)

                                                                      tCACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCAGGACTGACGCCCAGCAATAGACGacacac                   >  2:905926/1‑91 (MQ=255)

                                                                       cACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  <  2:4904932/93‑1 (MQ=255)

                                                                       cACTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  <  2:2134524/93‑1 (MQ=255)

                                                                        aCTGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcg                 >  2:2210082/1‑93 (MQ=255)

                                                                          tGAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgcc               <  2:1344917/93‑1 (MQ=255)

                                                                           gAGCGGGAAAATAGCGTCCCC‑AACCCCACGCCACCCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              <  2:4388995/93‑1 (MQ=255)

                                                                             gCGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccgcg            >  1:3767794/1‑93 (MQ=255)

                                                                              cGGGAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGa           <  1:2536080/93‑1 (MQ=255)

                                                                                ggAAAATAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCGGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  1:5149068/1‑93 (MQ=255)

                                                                                  aaaaTAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       >  1:2911242/1‑93 (MQ=255)

                                                                                  aaaaTAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       >  1:564373/1‑93 (MQ=255)

                                                                                   aaaTAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAcc      >  2:3989657/1‑93 (MQ=255)

                                                                                     aTAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  2:5559199/93‑1 (MQ=255)

                                                                                      tAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  2:50657/93‑1 (MQ=255)

                                                                                      tAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  1:1735864/93‑1 (MQ=255)

                                                                                      tAGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCACCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   >  1:1900482/1‑93 (MQ=255)

                                                                                      aaGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  1:2036193/92‑1 (MQ=255)

                                                                                       aGCGTCCCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTgg  <  2:1414993/93‑1 (MQ=255)

GCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCC‑AACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006389‑4006567

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