BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,343,823	C→T	G315E (GGG→GAG)	yciR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,343,823	0	C	T	98.5%	458.5 / ‑5.4	134	G315E (GGG→GAG)	yciR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (56/76); total (57/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

AAGGGGAGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGC  >  NC_000913/1343733‑1343912

 aGGGGAGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCAc                                                                                        <  2:4454996/93‑1 (MQ=255)

  ggggAGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCAcc                                                                                       <  2:4351837/93‑1 (MQ=255)

    ggAGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTg                                                                                     <  2:559187/93‑1 (MQ=255)

      aGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGa                                                                                   <  1:281344/93‑1 (MQ=255)

      aGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGa                                                                                   <  1:975080/93‑1 (MQ=255)

       gCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGAc                                                                                  >  1:4062967/1‑93 (MQ=255)

             aGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGAGCCa                                                                            >  1:3627396/1‑93 (MQ=255)

               cGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAac                                                                          <  1:1096617/93‑1 (MQ=255)

               cGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAac                                                                          >  1:3787638/1‑93 (MQ=255)

               cGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAac                                                                          <  2:1679141/93‑1 (MQ=255)

                 ggTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAacac                                                                        <  1:563011/93‑1 (MQ=255)

                  gTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACAcc                                                                       <  1:3027410/93‑1 (MQ=255)

                  gTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACAcc                                                                       <  1:2458580/93‑1 (MQ=255)

                   tCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACATCt                                                                      <  2:1408834/93‑1 (MQ=255)

                    cAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTg                                                                     <  1:1754299/93‑1 (MQ=255)

                    cAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTg                                                                     <  1:4564028/93‑1 (MQ=255)

                     aaaaTTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTgg                                                                    <  2:1779190/93‑1 (MQ=255)

                       aaTTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTc                                                                  <  2:2147801/93‑1 (MQ=255)

                        aTTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCa                                                                 <  1:4110821/93‑1 (MQ=255)

                        aTTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCa                                                                 <  1:3558616/93‑1 (MQ=255)

                         ttCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCAt                                                                <  2:2864972/93‑1 (MQ=255)

                          tCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATg                                                               >  1:2130259/1‑93 (MQ=255)

                           cGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGt                                                              <  1:1282513/93‑1 (MQ=255)

                           cGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGt                                                              <  1:1282512/93‑1 (MQ=255)

                           cGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGt                                                              <  2:747898/93‑1 (MQ=255)

                            gTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGtt                                                             >  2:934694/1‑93 (MQ=255)

                             tGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTc                                                            >  2:1360805/1‑93 (MQ=255)

                              gATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCg                                                           <  1:1233336/93‑1 (MQ=255)

                              gATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCg                                                           >  1:1669753/1‑93 (MQ=255)

                              gATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCg                                                           >  1:2883223/1‑93 (MQ=255)

                               aTGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCga                                                          <  1:3488086/93‑1 (MQ=255)

                               aTGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCga                                                          <  1:781673/93‑1 (MQ=255)

                               aTGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCga                                                          <  2:1699417/93‑1 (MQ=255)

                               aTGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCga                                                          <  2:2843312/93‑1 (MQ=255)

                               aTGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCga                                                          >  2:1832393/1‑93 (MQ=255)

                               aTGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCga                                                          >  1:4406015/1‑93 (MQ=255)

                                tGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCgag                                                         <  1:392921/93‑1 (MQ=255)

                                  ccATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGAc                                                       >  2:3403312/1‑93 (MQ=255)

                                  ccATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGAc                                                       <  1:541138/93‑1 (MQ=255)

                                  ccATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGAc                                                       >  2:2535553/1‑93 (MQ=255)

                                   cATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACa                                                      <  2:4443701/93‑1 (MQ=255)

                                    aTTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAg                                                     <  1:4351748/93‑1 (MQ=255)

                                    aTTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAg                                                     >  1:3273590/1‑93 (MQ=255)

                                     ttGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGc                                                    >  2:1948845/1‑93 (MQ=255)

                                     ttGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGc                                                    >  2:3775395/1‑93 (MQ=255)

                                     ttGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGc                                                    <  2:23013/93‑1 (MQ=255)

                                     ttGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGc                                                    >  2:103151/1‑93 (MQ=255)

                                     ttGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGc                                                    >  2:104331/1‑93 (MQ=255)

                                      tGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCt                                                   <  2:261113/93‑1 (MQ=255)

                                      tGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCt                                                   >  1:190380/1‑93 (MQ=255)

                                       tcTTCCAGCGCGCTGTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCCCCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCtt                                                  <  2:3837302/92‑1 (MQ=255)

                                          tCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTaaa                                               >  2:3096866/1‑93 (MQ=255)

                                           ccAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTaaaa                                              <  1:2750083/93‑1 (MQ=255)

                                            cAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAAt                                             <  2:897750/93‑1 (MQ=255)

                                              gcgcgcTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCACCAACACCTGGTCATGTTCGAGACAGCTTAAAATAg                                           <  2:4335472/93‑1 (MQ=255)

                                               cgcgcTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGc                                          >  2:4088319/1‑93 (MQ=255)

                                                 cgcTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCa                                        <  1:2269309/93‑1 (MQ=255)

                                                  gcTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCaa                                       >  1:4607905/1‑93 (MQ=255)

                                                  gcTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCaa                                       <  2:280609/93‑1 (MQ=255)

                                                  gcTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCCGCACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCaa                                       >  2:3260621/1‑93 (MQ=255)

                                                   cTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAAt                                      <  2:140638/93‑1 (MQ=255)

                                                    tttGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATg                                     <  1:1528545/93‑1 (MQ=255)

                                                     ttGTGAGGTGTTGGATGCCAGTACCAGAAAAACATCCTCACCAGGACGCACAAACACCTGGTAATGTTCGAGACAGCTTAAAATAGCCAATGa                                    <  2:4337399/93‑1 (MQ=255)

                                                      tgtgAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGac                                   >  1:3948921/1‑93 (MQ=255)

                                                      tgtgAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTCAAATAGCCCATGag                                   >  1:3948922/1‑92 (MQ=255)

                                                        tgAGGTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACCCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGacgc                                 <  2:1860082/93‑3 (MQ=255)

                                                         gAGGTGTTGGATGCCAGTACCAGAAACTCACCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcg                                <  1:1914263/93‑1 (MQ=255)

                                                           ggTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACCCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTc                              <  1:3926173/93‑1 (MQ=255)

                                                           ggTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTc                              >  2:2349446/1‑93 (MQ=255)

                                                            gtgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcg                             >  1:2211437/1‑93 (MQ=255)

                                                            gtgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcg                             >  2:125694/1‑93 (MQ=255)

                                                            gtgtTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACCCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTcg                             <  1:4408754/93‑1 (MQ=255)

                                                             tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                            >  2:3704123/1‑93 (MQ=255)

                                                             tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                            <  1:3290874/93‑1 (MQ=255)

                                                             tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                            >  1:1598231/1‑93 (MQ=255)

                                                             tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                            >  2:4063619/1‑93 (MQ=255)

                                                               ttGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcgtcgcgt                          >  1:211419/1‑93 (MQ=255)

                                                               gtGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcgtcgcgt                          <  1:1007347/92‑1 (MQ=255)

                                                                 ggATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaa                        <  2:4291614/93‑1 (MQ=255)

                                                                  gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAGCACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                       <  2:2677291/93‑1 (MQ=255)

                                                                  gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCTAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                       <  2:2445167/93‑1 (MQ=255)

                                                                  gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                       >  2:2082897/1‑93 (MQ=255)

                                                                  gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                       <  2:600974/93‑1 (MQ=255)

                                                                   aTGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaata                      >  1:1295272/1‑93 (MQ=255)

                                                                     gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                    <  2:2957054/93‑1 (MQ=255)

                                                                     gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                    <  1:2593472/93‑1 (MQ=255)

                                                                      ccAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACGCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACt                   <  2:4435158/93‑1 (MQ=255)

                                                                        aGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTgg                 >  2:3157887/1‑93 (MQ=255)

                                                                        aGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTgg                 <  1:205192/93‑1 (MQ=255)

                                                                         gTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGt                <  2:1136503/93‑1 (MQ=255)

                                                                         gTACCAGAAACTCATACTCACCTGGACGCGCCAGCACCTGGACATATTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTgat                <  1:4206106/93‑3 (MQ=255)

                                                                           aCCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCa              <  2:1059050/93‑1 (MQ=255)

                                                                             cAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCAcc            >  2:1906031/1‑93 (MQ=255)

                                                                              aGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACca           >  1:3078876/1‑93 (MQ=255)

                                                                              aGAAACTCATCCTCACCTGGACGCCCCAACACCTGGCCATGTTCGAGACAGCTTAAAATAGCCAATAACGCGTCGCGTAAAAACTGGTCACca           >  2:1244631/1‑93 (MQ=255)

                                                                                 aaCTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaac        <  1:1726878/93‑1 (MQ=255)

                                                                                 aaCTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaac        <  2:3345857/93‑1 (MQ=255)

                                                                                     cATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAAt    <  1:1385777/93‑1 (MQ=255)

                                                                                      aTCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATg   >  1:1703405/1‑93 (MQ=255)

                                                                                       tCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGc  >  2:3433024/1‑93 (MQ=255)

AAGGGGAGCCGTAAGCGGGTCAAAATTCGTGATGCCATTGCTTCCAGCGCGCTTTGTGAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGC  >  NC_000913/1343733‑1343912

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