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.6%	702.6 / 1.2	219	coding (186/207 nt)	b4554	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (3/0); new base C (139/77); total (142/77)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 242 total aligned reads displayed.

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

ggcgCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGt                                                                                          <  1:1850589/93‑1 (MQ=255)

 gcgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑c                                                                                        >  1:2754052/1‑93 (MQ=255)

  cgcgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGtac                                                                                        >  2:3522500/1‑91 (MQ=255)

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

   gcgcCACTGCTGGATTCCCAGCCAAATCAAGCAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccc                                                                                      >  1:5131793/1‑93 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                     <  2:2061292/93‑2 (MQ=255)

    cgcCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑cccc                                                                                     >  1:1776774/1‑92 (MQ=255)

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

       cACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑ccccaac                                                                                  <  1:3576801/93‑5 (MQ=255)

        aCTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAAacc                                                                                 >  1:2850802/1‑93 (MQ=255)

           gctgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCa                                                                               >  1:6634692/1‑93 (MQ=255)

             tgGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACg                                                                             <  2:6732183/93‑1 (MQ=255)

              ggATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGc                                                                            >  1:3309210/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGcc                                                                           >  1:526029/1‑93 (MQ=255)

               gATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGcc                                                                           >  2:4273983/1‑93 (MQ=255)

                 ttCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCaa                                                                         >  2:3948935/1‑93 (MQ=255)

                 tccccAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCaa                                                                         <  1:5702120/91‑1 (MQ=255)

                  tCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCCCCCCAAc                                                                        >  1:4253462/1‑93 (MQ=255)

                     cAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAg                                                                     >  1:6249349/1‑93 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAGACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGc                                                                    >  2:4697031/1‑93 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGc                                                                    <  2:3315468/93‑1 (MQ=255)

                      aGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGc                                                                    >  2:5013434/1‑93 (MQ=255)

                          aaaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACaa                                                                <  2:5721916/93‑1 (MQ=255)

                           aaTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAat                                                               >  2:5551234/1‑93 (MQ=255)

                            aTCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAata                                                              >  2:619034/1‑93 (MQ=255)

                                aGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACCCCAACCAGCACAATATGaa                                                          >  2:5599282/1‑93 (MQ=255)

                                   aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCg                                                       >  1:418032/1‑93 (MQ=255)

                                   aaGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCg                                                       <  2:12832/93‑1 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  2:1899561/1‑93 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  2:2768496/1‑93 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  1:5981761/1‑93 (MQ=255)

                                    aGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGc                                                      >  1:2479483/1‑93 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     <  2:6144895/93‑1 (MQ=255)

                                     gCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGcc                                                     >  1:6007127/1‑93 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    <  1:4723367/93‑1 (MQ=255)

                                      ccGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCa                                                    >  2:2693532/1‑93 (MQ=255)

                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  <  1:400742/93‑1 (MQ=255)

                                        gcgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGt                                                  >  1:359637/1‑93 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAACCGCAATCACTGAGAGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 <  1:3292592/93‑1 (MQ=255)

                                         cgcCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTc                                                 >  2:6279739/1‑93 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:1477322/1‑93 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  2:441828/1‑93 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              <  1:1332194/93‑1 (MQ=255)

                                            cTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCgg                                              >  1:1069137/1‑93 (MQ=255)

                                              gCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTc                                            >  2:2910813/1‑93 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           >  2:2537872/1‑93 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           <  1:3769929/93‑1 (MQ=255)

                                               cccACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCt                                           <  2:4630088/93‑1 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         <  2:119919/93‑1 (MQ=255)

                                                 cacttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGa                                         >  1:2147877/1‑93 (MQ=255)

                                                  ccttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  1:5576939/92‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        <  2:5177683/93‑1 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:858317/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  2:1700128/1‑93 (MQ=255)

                                                  acttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGaa                                        >  1:5799595/1‑93 (MQ=255)

                                                   cttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAg                                       <  1:1217084/93‑1 (MQ=255)

                                                    ttcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGc                                      <  2:4311133/93‑1 (MQ=255)

                                                     tcaacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGcc                                     >  2:1738176/1‑93 (MQ=255)

                                                       aacacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCa                                   <  2:588507/93‑1 (MQ=255)

                                                         cacttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGc                                 >  2:2261975/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:1437487/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCa                                >  2:317608/1‑93 (MQ=255)

                                                          acttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCCCCAGTCCGGACTGAAGCACAGCa                                >  1:2581556/1‑93 (MQ=255)

                                                           cttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCaa                               >  2:4860372/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  1:1989918/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  2:6444636/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  2:6620555/1‑93 (MQ=255)

                                                            ttcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAAt                              >  2:2377333/1‑93 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             >  2:2166790/1‑93 (MQ=255)

                                                             tcaaACGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATa                             <  2:505929/93‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           <  1:4125473/93‑1 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  2:5259546/1‑93 (MQ=255)

                                                               aaaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGa                           >  1:999266/1‑93 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:5317982/1‑93 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:2495565/1‑93 (MQ=255)

                                                                aaCGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGAc                          >  2:2100322/1‑93 (MQ=255)

                                                                  cGCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGc                        >  2:3393711/1‑93 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       <  2:6517787/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       <  2:1969968/93‑1 (MQ=255)

                                                                   gCAATCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGcc                       >  2:135482/1‑93 (MQ=255)

                                                                      aTCACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAcc                    <  1:1255373/93‑1 (MQ=255)

                                                                        cACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  >  1:449492/1‑93 (MQ=255)

                                                                        cACTGAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgc                  >  1:4565404/1‑93 (MQ=255)

                                                                            gAGCGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCAccgcgccg              >  2:3686791/1‑93 (MQ=255)

                                                                               cGGGAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGa           >  1:2419256/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  1:1685674/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  2:37419/1‑93 (MQ=255)

                                                                                gggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGAt          >  2:5916066/1‑93 (MQ=255)

                                                                                 ggAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATa         >  2:4953852/1‑93 (MQ=255)

                                                                                  gAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        <  1:4986804/93‑1 (MQ=255)

                                                                                  gAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        >  2:6521975/1‑93 (MQ=255)

                                                                                  gAAAATAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATaa        <  2:4652055/93‑1 (MQ=255)

                                                                                   aaaaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAAc       <  2:3884537/93‑1 (MQ=255)

                                                                                     aaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     >  2:1849095/1‑93 (MQ=255)

                                                                                     aaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     >  1:4472159/1‑93 (MQ=255)

                                                                                     aaTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCg     >  2:409691/1‑93 (MQ=255)

                                                                                      aTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    >  2:4389538/1‑93 (MQ=255)

                                                                                      aTAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    >  2:758902/1‑93 (MQ=255)

                                                                                      aTAGCGTCCCCAACCCCACGCAAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGt    <  2:909053/93‑1 (MQ=255)

                                                                                       tAGCGTCCCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTg   <  1:693046/93‑1 (MQ=255)

GGCGCGCCACTGCTGGATTCCCAGCCAAATCAAGTAAGCCGCGCCTGCCCACTTCAACACTTCAAACGCAATCACTGAGCGGGAAAATAGCGT‑CCCAACCCCACGCCAACCAGCACAATATGAATCGCCAGTCCGGTCTGAAGCCCAGCAATAGACGCCACCGCGCCGCGATAACCGTGG  >  NC_000913/4006388‑4006567

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