BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,114,866	T→C	N83S (AAC→AGC)	b2973 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,114,866	0	T	C	75.2%	478.2 / 204.7	295	N83S (AAC→AGC)	b2973	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (40/33); new base C (105/117); total (145/150)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 297 total aligned reads displayed.

GGATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATC  >  NC_000913/3114776‑3114958

ggATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTc                                                                                            <  1:2619258/93‑1 (MQ=255)

 gATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCa                                                                                           >  2:1478657/1‑93 (MQ=255)

 gATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCa                                                                                           <  2:1012862/93‑1 (MQ=255)

 gATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCa                                                                                           <  2:1288499/93‑1 (MQ=255)

 gATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACCTGCTCa                                                                                           <  2:6988028/93‑1 (MQ=255)

  aTAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAg                                                                                          >  1:661947/1‑93 (MQ=255)

     aCGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGc                                                                                       >  1:5993186/1‑93 (MQ=255)

      cGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGcc                                                                                      <  1:7349302/93‑1 (MQ=255)

      cGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGcc                                                                                      <  1:6557754/93‑1 (MQ=255)

       gTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCa                                                                                     <  2:2225368/93‑1 (MQ=255)

       gTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCa                                                                                     <  2:1366662/93‑1 (MQ=255)

        ttCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCAt                                                                                    <  2:2744717/93‑1 (MQ=255)

        ttCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCAt                                                                                    >  2:2502698/1‑93 (MQ=255)

        ttCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCAt                                                                                    <  2:5593962/93‑1 (MQ=255)

            tAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGAc                                                                                >  1:4335032/1‑93 (MQ=255)

            tAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCgnn                                                                                >  1:2185880/1‑91 (MQ=255)

            tAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGAc                                                                                <  1:2765108/93‑1 (MQ=255)

            tAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGAc                                                                                >  1:3040345/1‑93 (MQ=255)

            tAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGAc                                                                                >  1:616966/1‑93 (MQ=255)

             aGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACa                                                                               <  1:1848193/93‑1 (MQ=255)

             aGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACa                                                                               >  1:7275920/1‑93 (MQ=255)

             aGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACa                                                                               >  1:1216000/1‑93 (MQ=255)

              gACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAg                                                                              >  2:951690/1‑93 (MQ=255)

                cccAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGAc                                                                            >  1:344524/1‑93 (MQ=255)

                 ccAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGAcc                                                                           <  1:135725/93‑1 (MQ=255)

                 ccAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGAcc                                                                           <  1:5278350/93‑1 (MQ=255)

                 ccAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGAcc                                                                           <  1:7169632/93‑1 (MQ=255)

                  cAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCt                                                                          <  2:3654677/93‑1 (MQ=255)

                  cAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCt                                                                          <  1:5837661/93‑1 (MQ=255)

                   aaaTGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTg                                                                         >  1:1176564/1‑93 (MQ=255)

                    aaTGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGc                                                                        >  2:2993237/1‑93 (MQ=255)

                     aTGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCg                                                                       <  2:6664895/93‑1 (MQ=255)

                     aTGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCg                                                                       >  2:4658620/1‑93 (MQ=255)

                      tGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCgg                                                                      >  2:6165722/1‑93 (MQ=255)

                      tGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCgg                                                                      >  2:2254636/1‑93 (MQ=255)

                      tGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCAATCGCCAAACCTCCgc                                                                      >  2:4960677/1‑92 (MQ=255)

                         ccGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCAt                                                                   >  1:4959785/1‑93 (MQ=255)

                         ccGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCAt                                                                   >  2:3961557/1‑93 (MQ=255)

                         ccGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCAt                                                                   >  1:943120/1‑93 (MQ=255)

                           gTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATcc                                                                 >  2:4233056/1‑93 (MQ=255)

                           gTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATcc                                                                 >  2:3159366/1‑93 (MQ=255)

                               cgcgcTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACa                                                             <  1:2481426/93‑1 (MQ=255)

                                gcgcTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAg                                                            <  2:387839/93‑1 (MQ=255)

                                 cgcTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGa                                                           <  2:4096516/93‑1 (MQ=255)

                                     gctgACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGca                                                       <  1:3642439/93‑1 (MQ=255)

                                      ctgACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGcac                                                      <  1:4043769/93‑1 (MQ=255)

                                       tgACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGcaca                                                     >  1:2939828/1‑93 (MQ=255)

                                        gACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACaa                                                    >  2:1818927/1‑93 (MQ=255)

                                          cGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACaaaa                                                  >  1:3811859/1‑93 (MQ=255)

                                           gAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAAc                                                 <  1:1896334/93‑1 (MQ=255)

                                           gAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAAc                                                 <  1:2304767/93‑1 (MQ=255)

                                           gAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTCCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAAc                                                 <  1:7511620/93‑1 (MQ=255)

                                            aaCGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAAcc                                                >  2:6942947/1‑93 (MQ=255)

                                            aaCGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAAcc                                                >  2:6552808/1‑93 (MQ=255)

                                             aCGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCa                                               >  2:6784110/1‑93 (MQ=255)

                                             aCGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCa                                               <  2:7500416/93‑1 (MQ=255)

                                               gcgGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGa                                             <  2:1736286/93‑1 (MQ=255)

                                               gcgGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGa                                             <  1:6370006/93‑1 (MQ=255)

                                               gcgGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGa                                             <  2:4749598/93‑1 (MQ=255)

                                               gcgGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGa                                             <  2:5925931/93‑1 (MQ=255)

                                               gcgGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGa                                             >  1:423025/1‑93 (MQ=255)

                                                 ggTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGAcg                                           <  1:5726190/93‑1 (MQ=255)

                                                 ggTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGAcg                                           >  2:1878522/1‑93 (MQ=255)

                                                  gTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGAcgc                                          >  1:6416899/1‑93 (MQ=255)

                                                  gTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGAcgc                                          <  1:7232629/93‑1 (MQ=255)

                                                   ttCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGAcgcg                                         >  2:2838615/1‑93 (MQ=255)

                                                    tCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACgcgc                                        <  1:1670095/93‑1 (MQ=255)

                                                     cGGATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGct                                       <  2:2603770/93‑1 (MQ=255)

                                                      ggATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGctc                                      >  1:5004087/1‑93 (MQ=255)

                                                       gATACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGctct                                     <  2:4514215/93‑1 (MQ=255)

                                                        aTACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGctctc                                    <  1:5147695/93‑1 (MQ=255)

                                                        aTACCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGctctc                                    <  1:5445593/93‑1 (MQ=255)

                                                         tACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGctctct                                   >  2:6395227/1‑93 (MQ=255)

                                                          aCCCTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCtt                                  <  1:2948415/93‑1 (MQ=255)

                                                           cccTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTc                                 <  1:5206846/93‑1 (MQ=255)

                                                             cTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCt                               >  1:4571585/1‑93 (MQ=255)

                                                             cTTGCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCt                               <  2:3594552/93‑1 (MQ=255)

                                                              ttGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCtt                              <  1:7124440/93‑1 (MQ=255)

                                                                gCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTaa                            <  1:871872/93‑1 (MQ=255)

                                                                gCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTaa                            <  2:6297645/93‑1 (MQ=255)

                                                                gCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTaa                            <  1:3520997/93‑1 (MQ=255)

                                                                gCGGATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTaa                            <  1:6912330/93‑1 (MQ=255)

                                                                  ggATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGa                          <  2:5227678/93‑1 (MQ=255)

                                                                   gATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGAt                         <  1:5499720/93‑1 (MQ=255)

                                                                   gATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGAt                         >  2:5808910/1‑93 (MQ=255)

                                                                   gATCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGAt                         <  2:6626925/93‑1 (MQ=255)

                                                                    aTCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGAtt                        >  2:332703/1‑93 (MQ=255)

                                                                    aTCGTTATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGAtt                        <  2:5660489/93‑1 (MQ=255)

                                                                        ttATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCtc                    >  2:5472611/1‑93 (MQ=255)

                                                                        ttATTTACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCtc                    <  2:1409679/93‑1 (MQ=255)

                                                                           tttACTACCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCtcatc                 <  2:3807734/93‑1 (MQ=255)

                                                                            ttACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCtcatca                <  1:5598909/93‑1 (MQ=255)

                                                                             tactacCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCAc               <  1:3307278/93‑1 (MQ=255)

                                                                               ctacCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGt             >  1:6798225/1‑93 (MQ=255)

                                                                                tacCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGtt            >  2:1990367/1‑93 (MQ=255)

                                                                                tacCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGtt            >  1:488589/1‑93 (MQ=255)

                                                                                 acCGACTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGttt           <  2:2672793/93‑1 (MQ=255)

                                                                                     aCTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCa       <  2:7314129/93‑1 (MQ=255)

                                                                                      cTTGCTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCAt      <  2:82028/93‑1 (MQ=255)

                                                                                          ttCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATc  >  1:6954222/1‑93 (MQ=255)

GGATAACGTTCATAGACCCAAATGCCCGTTGCGCGCTGCTGACGAACGCGGTTCGGATACCCTTGCGGATCGTTATTTACTACCGACTTGTTCAGTGCCATCGACAGACCTGCGGCATCCAACAGACGCACAAAACCAGACGCGCTCTCTTCCTTAAGATTGCTCATCACGTTTTCCATGATC  >  NC_000913/3114776‑3114958

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