BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,469,173	T→C	D122G (GAC→GGC)	b2354 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,469,173	0	T	C	99.6%	746.4 / ‑5.1	228	D122G (GAC→GGC)	b2354	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (103/124); total (103/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 237 total aligned reads displayed.

CACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTC  >  NC_000913/2469083‑2469262

cACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGccc                                                                                         <  1:3342049/93‑1 (MQ=255)

 aCCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCa                                                                                        >  2:4555816/1‑93 (MQ=255)

   cgcgCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGa                                                                                      >  2:3399459/1‑93 (MQ=255)

    gcgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGAt                                                                                     <  2:898236/93‑1 (MQ=255)

    gcgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGAt                                                                                     >  2:4150271/1‑93 (MQ=255)

    gcgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGAt                                                                                     <  1:169367/93‑1 (MQ=255)

     cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATa                                                                                    <  1:6423698/93‑1 (MQ=255)

      gcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATaa                                                                                   <  2:3888526/93‑1 (MQ=255)

       ccGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAAt                                                                                  <  1:2586840/93‑1 (MQ=255)

        cGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATc                                                                                 >  2:6690015/1‑93 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                <  2:1021964/93‑1 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                >  2:169969/1‑93 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                >  1:4022984/1‑93 (MQ=255)

          gATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCa                                                                               <  2:6086478/93‑1 (MQ=255)

           aTGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAg                                                                              <  2:2812085/93‑1 (MQ=255)

             gTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGcc                                                                            <  2:5615964/93‑1 (MQ=255)

             gTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGcc                                                                            <  1:1493887/93‑1 (MQ=255)

             gTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGcc                                                                            <  2:2412575/93‑1 (MQ=255)

               cATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCcac                                                                          >  1:3500537/1‑93 (MQ=255)

               cATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCcac                                                                          <  1:6537809/93‑1 (MQ=255)

                 tCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACaa                                                                        >  1:4790918/1‑93 (MQ=255)

                  cTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACaaa                                                                       <  1:2297383/93‑1 (MQ=255)

                    taccTGAACTGCCGGAGGAGTGGGCCATTCAATATCTGTCGCACTGGAAGTATCGGCCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAAtt                                                                     <  1:1118170/91‑1 (MQ=17)

                    ggCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAAtt                                                                     <  2:6570375/93‑1 (MQ=255)

                     gCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATta                                                                    >  1:1434288/1‑93 (MQ=255)

                       cTGAACTGCCGGATGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATtata                                                                  <  1:2861684/93‑1 (MQ=255)

                       cTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATtata                                                                  >  1:4060747/1‑93 (MQ=255)

                         gAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATAtt                                                                <  2:3116676/93‑1 (MQ=255)

                         gAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATAtt                                                                >  1:5044223/1‑93 (MQ=255)

                           accgCCGGGGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTgc                                                              <  1:2673432/90‑1 (MQ=255)

                            cTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTgcg                                                             <  1:3726668/93‑1 (MQ=255)

                             tGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTgcgc                                                            >  2:3785172/1‑93 (MQ=255)

                             tGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTgcgc                                                            <  1:516304/93‑1 (MQ=255)

                                 ggaggaGTAGCCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGt                                                        <  1:290079/93‑1 (MQ=255)

                                  gaggagTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGtt                                                       >  2:5795773/1‑93 (MQ=255)

                                   aggagTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTc                                                      <  1:3609551/93‑1 (MQ=255)

                                   aggagTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTc                                                      >  2:5094778/1‑93 (MQ=255)

                                    ggagTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTcc                                                     <  1:5216585/93‑1 (MQ=255)

                                       gTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCa                                                  <  1:241113/93‑1 (MQ=255)

                                       gTAGGCCATTCAATATCTGGCACACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCa                                                  <  2:2326400/93‑1 (MQ=255)

                                         aGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCAcc                                                >  2:1719797/1‑93 (MQ=255)

                                         aGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCAcc                                                <  1:4287185/93‑1 (MQ=255)

                                         aGGCCATTAAATATCTGGCAGATTGAAAGTATCGACAAGTTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCAcc                                                <  2:1245743/93‑1 (MQ=18)

                                          ggCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCt                                               >  1:4645367/1‑93 (MQ=255)

                                            ccATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCttt                                             <  2:849343/93‑1 (MQ=255)

                                             cATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTc                                            <  2:5384484/93‑1 (MQ=255)

                                             cATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTc                                            >  1:1299371/1‑93 (MQ=255)

                                              aTTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCa                                           >  2:4731053/1‑93 (MQ=255)

                                              aTTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCa                                           <  1:2136052/93‑1 (MQ=255)

                                               ttCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAg                                          >  2:2060293/1‑93 (MQ=255)

                                                 cAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgac                                        <  2:4877256/93‑1 (MQ=255)

                                                  aaTATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacg                                       >  2:1077280/1‑93 (MQ=255)

                                                   atatCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                      >  1:2958411/1‑93 (MQ=255)

                                                    tatCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacgac                                     <  1:4261679/93‑1 (MQ=255)

                                                     atCTGGCGCACTGGAAGTATCGACCAGCTCAAGTGCGCCCAGATAACCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGAcc                                    <  1:4138084/93‑1 (MQ=255)

                                                      tCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCa                                   >  1:4447728/1‑93 (MQ=255)

                                                        tGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAAt                                 >  1:5904781/1‑93 (MQ=255)

                                                        tGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAAt                                 >  1:6706650/1‑93 (MQ=255)

                                                         ggCGCACTGGAAGTATCGCCCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATa                                <  1:5950952/93‑1 (MQ=255)

                                                         ggCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATa                                <  2:6642751/93‑1 (MQ=255)

                                                          gcgcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAg                               >  1:6727598/1‑93 (MQ=255)

                                                          gcgcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAg                               <  2:34608/93‑1 (MQ=255)

                                                           cgcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGc                              <  2:1459059/93‑1 (MQ=255)

                                                           cgcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGc                              >  2:2694288/1‑93 (MQ=255)

                                                            gcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGcc                             <  2:6403404/93‑1 (MQ=255)

                                                             cACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCg                            >  1:5650975/1‑93 (MQ=255)

                                                               cTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCt                          <  2:301870/93‑1 (MQ=255)

                                                                tGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCtt                         >  1:3231331/1‑93 (MQ=255)

                                                                tGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCtt                         <  2:5900232/93‑1 (MQ=255)

                                                                tGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCtt                         <  2:3957583/93‑1 (MQ=255)

                                                                 ggAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCttt                        <  1:6309411/93‑1 (MQ=255)

                                                                 ggAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCttt                        <  2:5328358/93‑1 (MQ=255)

                                                                   aaGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAc                      <  2:5929061/93‑1 (MQ=255)

                                                                    aGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAcc                     <  1:6170113/93‑1 (MQ=255)

                                                                    aGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAcc                     <  2:4349416/93‑1 (MQ=255)

                                                                     gTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCa                    >  2:6317688/1‑93 (MQ=255)

                                                                     gTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCa                    >  1:3441133/1‑93 (MQ=255)

                                                                      tATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAg                   <  2:2334857/93‑1 (MQ=255)

                                                                        tCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGcc                 <  2:1149150/93‑1 (MQ=255)

                                                                          gACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCa               >  1:535053/1‑93 (MQ=255)

                                                                          gACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCa               <  2:6260347/93‑1 (MQ=255)

                                                                           acgaGCTCGAGAGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCAt              <  2:4200083/90‑1 (MQ=38)

                                                                            ccAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCAtt             <  2:5708308/93‑1 (MQ=255)

                                                                            ccAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCAtt             >  2:2900654/1‑93 (MQ=255)

                                                                             cAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTg            <  2:5038983/93‑1 (MQ=255)

                                                                             cAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTg            >  1:991922/1‑93 (MQ=255)

                                                                              aGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGt           >  1:5486522/1‑93 (MQ=255)

                                                                              aGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGt           <  1:4340213/93‑1 (MQ=255)

                                                                               gCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGtt          >  2:2496244/1‑93 (MQ=255)

                                                                                 tCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa        <  1:4586757/93‑1 (MQ=255)

                                                                                 tCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa        <  1:53668/93‑1 (MQ=255)

                                                                                  ccAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAc       >  2:4221983/1‑93 (MQ=255)

                                                                                  ccAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAc       <  1:1831731/93‑1 (MQ=255)

                                                                                  ccAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAc       <  2:3647102/93‑1 (MQ=255)

                                                                                   cAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACt      <  2:3053862/93‑1 (MQ=255)

                                                                                     gTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGt    <  1:1057845/93‑1 (MQ=255)

                                                                                     gTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGt    <  2:3330574/93‑1 (MQ=255)

                                                                                     gTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGt    >  1:574137/1‑93 (MQ=255)

                                                                                      tGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGtt   <  1:1849702/93‑1 (MQ=255)

                                                                                       gcgcCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCCCCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTc  <  1:3991426/93‑1 (MQ=255)

CACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTC  >  NC_000913/2469083‑2469262

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