BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,701,951	T→C	S128G (AGC→GGC)	b1624 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,701,951	0	T	C	98.6%	453.1 / ‑3.2	144	S128G (AGC→GGC)	b1624	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (69/73); total (70/74)
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.57e-01

Only 100 of 150 total aligned reads displayed.

GCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGG  >  NC_000913/1701861‑1702040

   gacggaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATccc                                                                                      <  2:4107355/93‑1 (MQ=255)

   gacggaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATccc                                                                                      <  2:3482705/93‑1 (MQ=255)

    acggaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCa                                                                                     <  2:2863544/93‑1 (MQ=255)

      ggaAGCGGTCAAAATGAGACTCAAAGTACGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGc                                                                                   <  1:4213555/93‑1 (MQ=255)

      ggaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGc                                                                                   <  1:2668427/93‑1 (MQ=255)

       gaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCg                                                                                  >  2:2149627/1‑93 (MQ=255)

       gaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCg                                                                                  <  2:4491843/93‑1 (MQ=255)

       gaAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCg                                                                                  >  2:3628068/1‑93 (MQ=255)

        aaGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGa                                                                                 <  1:2160386/93‑1 (MQ=255)

        aaGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGa                                                                                 >  2:2692010/1‑93 (MQ=255)

         aGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCTCAGCGAc                                                                                >  2:3830742/1‑93 (MQ=255)

         aGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGAc                                                                                >  1:1329620/1‑93 (MQ=255)

         aGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGAc                                                                                >  2:2731219/1‑93 (MQ=255)

         aGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGAc                                                                                >  2:3033229/1‑93 (MQ=255)

            ggTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGt                                                                             <  1:2000459/93‑1 (MQ=255)

            ggTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGt                                                                             >  1:639804/1‑93 (MQ=255)

            ggTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGt                                                                             <  2:35303/93‑1 (MQ=255)

             gTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGtt                                                                            <  1:3206210/93‑1 (MQ=255)

             gTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGtt                                                                            >  1:3649255/1‑93 (MQ=255)

              tCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTa                                                                           >  1:3794137/1‑93 (MQ=255)

               cAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTAt                                                                          <  1:818684/93‑1 (MQ=255)

               cAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTAt                                                                          <  2:4438892/93‑1 (MQ=255)

                 aaaTGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATgg                                                                        <  1:1649036/93‑1 (MQ=255)

                 aaaTGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATgg                                                                        <  1:3779197/93‑1 (MQ=255)

                   aTGAGACTCAAAGTAAGCAACTTCACCCAGCACGCATTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGaa                                                                      >  1:1318080/1‑93 (MQ=255)

                    tGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAAt                                                                     >  2:1869411/1‑93 (MQ=255)

                      agaCTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATac                                                                   <  1:1531753/93‑1 (MQ=255)

                       gaCTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATAGCCATCCCAGCGACGGTTATGGAATaca                                                                  >  1:3992255/1‑93 (MQ=255)

                       acaacAAAGAAACCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATaca                                                                  <  1:2174838/89‑1 (MQ=255)

                          tCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATacagac                                                               <  1:129856/93‑1 (MQ=255)

                          tCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATacagac                                                               <  1:129857/93‑1 (MQ=255)

                           cAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATacagaca                                                              <  2:4525938/93‑1 (MQ=255)

                           cAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATacagaca                                                              >  1:3140790/1‑93 (MQ=255)

                                tAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGcaca                                                         <  1:2740883/93‑1 (MQ=255)

                                tAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGcaca                                                         >  1:1923141/1‑93 (MQ=255)

                                 aaGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCacac                                                        <  2:2267054/93‑1 (MQ=255)

                                    cAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGcc                                                     <  2:3336518/93‑1 (MQ=255)

                                    cAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGcc                                                     <  2:386216/93‑1 (MQ=255)

                                     aaCTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGccc                                                    >  1:991773/1‑93 (MQ=255)

                                      aCTTCACCCAGCACGCCTTTCGCGAGTAAACCTTTTAGCGTCAAGAAATCTCCATCCCAGCGTCGGTTATGGAATACATACAGCAAACGcccc                                                   >  2:2872157/1‑93 (MQ=255)

                                      aCTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGcccc                                                   <  1:1016744/93‑1 (MQ=255)

                                      aCTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGcccc                                                   <  1:3219819/93‑1 (MQ=255)

                                      aCTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCCCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGcccc                                                   <  1:1016743/93‑1 (MQ=255)

                                       cTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCa                                                  >  1:2246243/1‑93 (MQ=255)

                                        ttCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAg                                                 >  2:2267904/1‑93 (MQ=255)

                                        ttCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGGTATGGAATACAGACAGCACACGCCCCAg                                                 >  1:3383198/1‑93 (MQ=255)

                                        ctcccccAGCCCGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAg                                                 <  2:1460661/89‑1 (MQ=255)

                                          cACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGc                                               <  1:3986636/93‑1 (MQ=255)

                                          cACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGc                                               >  1:2787496/1‑93 (MQ=255)

                                           aCCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCt                                              >  2:1722676/1‑93 (MQ=255)

                                            cccAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCtt                                             >  2:119511/1‑93 (MQ=255)

                                            cccAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCtt                                             <  2:2922898/93‑1 (MQ=255)

                                            cccAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCtt                                             <  1:118829/93‑1 (MQ=255)

                                            cccAGCACGCCTTCCGCGAGTAAACCTTTTAGCGGCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCtt                                             <  1:118827/93‑1 (MQ=255)

                                              cAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCtttt                                           <  1:4733796/93‑1 (MQ=255)

                                                  aCGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGcc                                       >  2:3463507/1‑93 (MQ=255)

                                                    gCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAg                                     >  1:2357324/1‑93 (MQ=255)

                                                        tCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCa                                 <  1:1393121/93‑1 (MQ=255)

                                                        tCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCa                                 <  1:2523220/93‑1 (MQ=255)

                                                        tCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCa                                 <  1:3474533/93‑1 (MQ=255)

                                                          cgcgAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCCCACGCCCCAGGCTTTTTGCCAGCGCATc                               <  2:3381328/93‑1 (MQ=255)

                                                          cgcgAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATc                               <  2:4320282/93‑1 (MQ=255)

                                                           gcgAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATcc                              >  2:2954370/1‑93 (MQ=255)

                                                           gcgAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATcc                              <  2:4037637/93‑1 (MQ=255)

                                                           gcgAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATcc                              >  1:957686/1‑93 (MQ=255)

                                                            cgAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCa                             <  2:2180431/93‑1 (MQ=255)

                                                             gAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGGATACAGACAGCACACGCCCCAGGCTTTTTGACGGCGCATTCAg                            >  2:4324091/1‑93 (MQ=255)

                                                             gAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAg                            >  2:3320793/1‑93 (MQ=255)

                                                             gAGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAg                            >  2:1676576/1‑93 (MQ=255)

                                                              aGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGTCGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGc                           <  1:3202974/93‑1 (MQ=255)

                                                              aGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGc                           >  1:2556970/1‑93 (MQ=255)

                                                              aGTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGc                           <  2:4152480/93‑1 (MQ=255)

                                                               gTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGct                          >  1:1133292/1‑93 (MQ=255)

                                                               gTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGct                          <  1:557941/93‑1 (MQ=255)

                                                               gTAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGct                          <  2:1500025/93‑1 (MQ=255)

                                                                tAAACCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGctc                         >  2:4544676/1‑93 (MQ=255)

                                                                 aagccTGTTAGCGTCAAGAATTCGCCATCCCAGCGACGGTGAGGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGctct                        <  2:2617432/90‑1 (MQ=255)

                                                                 aaaCCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGctct                        <  2:534955/93‑1 (MQ=255)

                                                                 aaaCCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGctct                        >  2:3262241/1‑93 (MQ=255)

                                                                 aaaCCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGctct                        >  1:730322/1‑93 (MQ=255)

                                                                  aaCCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCtctc                       >  2:202865/1‑93 (MQ=255)

                                                                  aaCCTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCtctc                       <  1:870068/93‑1 (MQ=255)

                                                                    ccTTTTAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTcgc                     >  1:2958622/1‑93 (MQ=255)

                                                                      ttttAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCgcgc                   >  1:2666336/1‑93 (MQ=255)

                                                                      ttttAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCgcgc                   >  1:500219/1‑93 (MQ=255)

                                                                      ttttAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCgcgc                   >  1:3717101/1‑93 (MQ=255)

                                                                        ttAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCtt                 <  2:3840843/93‑1 (MQ=255)

                                                                         tAGCGTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTtg                <  2:3841663/93‑1 (MQ=255)

                                                                             gTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGAc            <  2:3066042/93‑1 (MQ=255)

                                                                             gTCAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGAc            <  1:472206/93‑1 (MQ=255)

                                                                               cAAGAAATCGCCATGCCAGCGACGGTTATGGAATACATACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAg          <  1:3890551/93‑1 (MQ=255)

                                                                               cAAGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAg          <  1:3890552/93‑1 (MQ=255)

                                                                                aaGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAgt         <  1:1368370/93‑1 (MQ=255)

                                                                                aaGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAgt         >  2:3480969/1‑93 (MQ=255)

                                                                                 aGAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAgtg        >  2:512240/1‑93 (MQ=255)

                                                                                  gAAATCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAgtgt       >  1:3018675/1‑93 (MQ=255)

                                                                                     aTCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCAc    <  1:1094102/93‑1 (MQ=255)

                                                                                     aTCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCAc    <  2:2799029/93‑1 (MQ=255)

                                                                                      tCGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACg   >  1:4243821/1‑93 (MQ=255)

                                                                                       cGCCATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACgg  >  2:1142503/1‑93 (MQ=255)

GCGGACGGAAGCGGTCAAAATGAGACTCAAAGTAAGCAACTTCACCCAGCACGCCTTCCGCGAGTAAACCTTTTAGCGTCAAGAAATCGCTATCCCAGCGACGGTTATGGAATACAGACAGCACACGCCCCAGGCTTTTTGCCAGCGCATCCAGCTCTCGCGCTTGTGACAGTGTCACGG  >  NC_000913/1701861‑1702040

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