BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,934,800	T→C	L42P (CTG→CCG)	yebK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,934,800	0	T	C	85.3%	578.9 / 93.5	252	L42P (CTG→CCG)	yebK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (21/16); new base C (107/108); total (128/124)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-02

Only 100 of 253 total aligned reads displayed.

GGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACA  >  NC_000913/1934708‑1934890

    cATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGaa                                                                                        <  1:3386405/93‑1 (MQ=255)

    cATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGaa                                                                                        <  1:6137498/93‑1 (MQ=255)

    cATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGaa                                                                                        >  1:1078259/1‑93 (MQ=255)

      tttGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGc                                                                                      >  2:5618763/1‑93 (MQ=255)

       ttGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGcc                                                                                     >  1:3853425/1‑93 (MQ=255)

           gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATg                                                                                 >  1:1285350/1‑93 (MQ=255)

           gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATg                                                                                 >  2:6651380/1‑93 (MQ=255)

           gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATg                                                                                 >  1:1896282/1‑93 (MQ=255)

           gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATg                                                                                 >  2:3189358/1‑93 (MQ=255)

           gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATg                                                                                 <  1:2158724/93‑1 (MQ=255)

             aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGtt                                                                               >  1:4942854/1‑93 (MQ=255)

             aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGtt                                                                               <  2:5557967/93‑1 (MQ=255)

             aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGtt                                                                               >  2:4978114/1‑93 (MQ=255)

              aaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTa                                                                              <  2:4387166/93‑1 (MQ=255)

              aaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTa                                                                              <  2:2996793/93‑1 (MQ=255)

               aTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAg                                                                             <  1:681687/93‑1 (MQ=255)

                tCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGc                                                                            >  2:6046595/1‑93 (MQ=255)

                  agagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGa                                                                          >  1:4282053/1‑93 (MQ=255)

                  agagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGa                                                                          <  2:7153196/93‑1 (MQ=255)

                  agagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGa                                                                          <  2:7496163/93‑1 (MQ=255)

                   gagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGaa                                                                         >  2:5972561/1‑93 (MQ=255)

                     gcgcAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAAcc                                                                       <  2:6368883/93‑1 (MQ=255)

                     gcgcAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAAcc                                                                       <  1:4632040/93‑1 (MQ=255)

                       gcAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGa                                                                     >  1:6241947/1‑93 (MQ=255)

                        cAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGAc                                                                    <  1:3467429/93‑1 (MQ=255)

                        cAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGAc                                                                    <  2:5029744/93‑1 (MQ=255)

                        cAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGAc                                                                    >  2:2198886/1‑93 (MQ=255)

                         aaaGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACg                                                                   <  1:3743196/93‑1 (MQ=255)

                           aGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGt                                                                 <  1:1322148/93‑1 (MQ=255)

                            gTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTg                                                                >  2:360370/1‑93 (MQ=255)

                            gTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTg                                                                >  1:414825/1‑93 (MQ=255)

                              tGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGaa                                                              <  1:2575774/93‑1 (MQ=255)

                              tGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGaa                                                              >  2:3925421/1‑93 (MQ=255)

                                ccGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATc                                                            <  1:1626079/93‑1 (MQ=255)

                                ccGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATc                                                            >  2:4055639/1‑93 (MQ=255)

                                 cGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCg                                                           >  2:899870/1‑93 (MQ=255)

                                  gAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGt                                                          >  1:2680791/1‑93 (MQ=255)

                                    ggTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGttt                                                        >  2:7044891/1‑93 (MQ=255)

                                    ggTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGttt                                                        >  1:531510/1‑93 (MQ=255)

                                     gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTc                                                       <  1:1027376/93‑1 (MQ=255)

                                     gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTc                                                       >  2:2503432/1‑93 (MQ=255)

                                      tCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCt                                                      <  2:1842773/93‑1 (MQ=255)

                                      tCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCt                                                      <  2:3768635/93‑1 (MQ=255)

                                        aTTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGt                                                    >  1:557307/1‑93 (MQ=255)

                                         ttCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTc                                                   >  1:1695283/1‑93 (MQ=255)

                                         ttCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTc                                                   >  1:5064566/1‑93 (MQ=255)

                                         ttCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTc                                                   <  2:1950109/93‑1 (MQ=255)

                                          tCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCg                                                  <  1:2941959/93‑1 (MQ=255)

                                            tGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgca                                                <  1:5262871/93‑1 (MQ=255)

                                                   gCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGa                                         <  1:4817798/93‑1 (MQ=255)

                                                   gCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGa                                         >  2:54577/1‑93 (MQ=255)

                                                    cccGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGac                                        <  2:1295783/93‑1 (MQ=255)

                                                     ccGATAACGCGGTCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGaca                                       >  2:4656173/1‑93 (MQ=255)

                                                     ccGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGaca                                       <  1:7206/93‑1 (MQ=255)

                                                      cGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGacac                                      >  1:3264726/1‑93 (MQ=255)

                                                        aTAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgc                                    <  1:398137/93‑1 (MQ=255)

                                                        aTAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgc                                    <  2:6354874/93‑1 (MQ=255)

                                                         tAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcg                                   >  2:6894524/1‑93 (MQ=255)

                                                         tAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcg                                   <  1:653158/93‑1 (MQ=255)

                                                          aaCGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcgc                                  <  2:2829922/93‑1 (MQ=255)

                                                           aCGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcgcg                                 >  2:1757631/1‑93 (MQ=255)

                                                           aCGCGATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcgcg                                 <  1:3005248/93‑1 (MQ=255)

                                                            cgcgATCCATTCGAGTATTGGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGGGAATCGTTTCTGTCGCAGCATGGACAcgcgcg                                 >  1:1147404/1‑93 (MQ=255)

                                                             gcgTTCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGt                               >  1:6973165/1‑93 (MQ=255)

                                                             gcgATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGt                               <  2:7109311/93‑1 (MQ=255)

                                                              cgATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGtt                              <  2:2643890/93‑1 (MQ=255)

                                                               gATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGttt                             >  2:1665607/1‑93 (MQ=255)

                                                               gATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGttt                             >  1:321291/1‑93 (MQ=255)

                                                               gATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGttt                             >  2:2626229/1‑93 (MQ=255)

                                                               gATCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGttt                             <  1:6100600/93‑1 (MQ=255)

                                                                 tCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTc                           >  1:882249/1‑93 (MQ=255)

                                                                 tCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTc                           >  1:4513553/1‑93 (MQ=255)

                                                                 tCCATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTc                           <  2:5040628/93‑1 (MQ=255)

                                                                  ccATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                          <  2:5163462/93‑1 (MQ=255)

                                                                  ccATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                          <  2:2712424/93‑1 (MQ=255)

                                                                  ccATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                          <  1:502744/93‑1 (MQ=255)

                                                                  ccATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                          >  1:2249973/1‑93 (MQ=255)

                                                                  ccATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                          <  2:6932721/93‑1 (MQ=255)

                                                                   cATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCt                         <  1:887918/93‑1 (MQ=255)

                                                                   cATTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCt                         >  1:7501751/1‑93 (MQ=255)

                                                                    aTTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                        >  1:5913184/1‑93 (MQ=255)

                                                                    aTTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                        <  1:6654382/93‑1 (MQ=255)

                                                                    aTTCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                        >  2:2392552/1‑93 (MQ=255)

                                                                     ttCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                       <  2:397759/93‑1 (MQ=255)

                                                                     ttCGAGTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                       <  2:350486/93‑1 (MQ=255)

                                                                          gTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTa                  >  1:2068304/1‑93 (MQ=255)

                                                                          gTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTa                  >  1:7122522/1‑93 (MQ=255)

                                                                          gTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTa                  <  2:303414/93‑1 (MQ=255)

                                                                          gTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTa                  >  2:3374678/1‑93 (MQ=255)

                                                                          gTATT‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTa                  >  2:1120210/1‑93 (MQ=255)

                                                                           tATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaa                 >  2:2227957/1‑93 (MQ=255)

                                                                             tt‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAAc               >  1:994441/1‑93 (MQ=255)

                                                                             tt‑GCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAAc               <  2:4246026/93‑1 (MQ=255)

                                                                                 ctgctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTc            <  1:677514/93‑1 (MQ=255)

                                                                                  tgctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCa           >  1:1205013/1‑93 (MQ=255)

                                                                                   gctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCAt          <  1:78397/93‑1 (MQ=255)

                                                                                    ctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATc         <  2:5076539/93‑1 (MQ=255)

                                                                                       tGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTgg      <  2:4820878/93‑1 (MQ=255)

                                                                                       tGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTgg      >  1:7320702/1‑93 (MQ=255)

                                                                                        ggCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGc     >  2:7526179/1‑93 (MQ=255)

GGAACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATT‑GCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCACA  >  NC_000913/1934708‑1934890

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