BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,862,881	T→C	L208P (CTA→CCA)	ygbM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,862,881	0	T	C	99.0%	608.2 / ‑3.4	190	L208P (CTA→CCA)	ygbM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (84/104); total (85/105)
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.16e-01

Only 100 of 201 total aligned reads displayed.

ATACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862791‑2862970

 tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        >  1:614675/1‑93 (MQ=255)

 tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        <  2:1077520/93‑1 (MQ=255)

  aCTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCt                                                                                       <  1:1356835/93‑1 (MQ=255)

  aCTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCt                                                                                       <  2:1467781/93‑1 (MQ=255)

   cTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGCAATATGCGCACGTACAAATTGCCGGACCACCTg                                                                                      >  2:3408981/1‑93 (MQ=255)

    ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGa                                                                                     >  2:4026083/1‑93 (MQ=255)

    ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     >  2:1780867/1‑93 (MQ=255)

    ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACATGa                                                                                     <  2:814967/93‑1 (MQ=255)

      ttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATc                                                                                   <  2:3948402/93‑1 (MQ=255)

       tCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCg                                                                                  <  1:4154322/93‑1 (MQ=255)

        cATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGa                                                                                 <  2:3169742/93‑1 (MQ=255)

         aTGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGAc                                                                                >  1:2111254/1‑93 (MQ=255)

          tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               >  1:3798225/1‑93 (MQ=255)

            cgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATg                                                                             <  1:4231343/93‑1 (MQ=255)

            cgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATg                                                                             <  1:234274/93‑1 (MQ=255)

             ccAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGa                                                                            <  1:3369239/92‑1 (MQ=255)

              cAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGaa                                                                           <  2:3493341/93‑1 (MQ=255)

              cAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGaa                                                                           <  1:283201/93‑1 (MQ=255)

               aaaaaGTAGATGGCAACCTGACGGATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGagc                                                                          >  1:500431/1‑91 (MQ=255)

               aaaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAc                                                                          >  2:4572876/1‑93 (MQ=255)

                aaaaGTAGATGGCAGCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAcc                                                                         <  1:4535275/93‑1 (MQ=255)

                aaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAcc                                                                         >  2:335500/1‑93 (MQ=255)

                 aaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCg                                                                        <  1:4115631/93‑1 (MQ=255)

                 aaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCg                                                                        >  1:223604/1‑93 (MQ=255)

                  aaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCgg                                                                       >  1:2538205/1‑93 (MQ=255)

                   aGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGa                                                                      <  2:2133796/93‑1 (MQ=255)

                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  2:2280006/1‑93 (MQ=255)

                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  1:1451109/1‑93 (MQ=255)

                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  2:1082176/1‑93 (MQ=255)

                    gTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGAc                                                                     >  1:3836342/1‑93 (MQ=255)

                      aGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGa                                                                   <  1:404034/93‑1 (MQ=255)

                      aGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGa                                                                   <  2:4502048/93‑1 (MQ=255)

                       gATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCTGACCACCTGATCGACATGAACCGGACGAt                                                                  <  2:3879911/93‑1 (MQ=255)

                       gATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGAt                                                                  >  2:1908502/1‑93 (MQ=255)

                        aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 <  1:2035397/93‑1 (MQ=255)

                        aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 <  1:1584251/93‑1 (MQ=255)

                             aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            >  2:806166/1‑93 (MQ=255)

                             aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGAACACCTGATCGAAATGCACCGGACGAAGGAGaa                                                            >  1:2507809/1‑93 (MQ=255)

                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCCCCTGATCGACATGAACCGGACGATGGAGaaa                                                           <  1:1208713/93‑1 (MQ=255)

                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           <  2:1615275/93‑1 (MQ=255)

                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           <  1:1645448/93‑1 (MQ=255)

                              aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           >  1:1008105/1‑93 (MQ=255)

                                cTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAAta                                                         >  1:2713003/1‑92 (MQ=255)

                                cTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATc                                                         >  2:1704103/1‑93 (MQ=255)

                                cTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATc                                                         >  1:1020771/1‑93 (MQ=255)

                                 tGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCa                                                        >  1:3250853/1‑93 (MQ=255)

                                 tGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCa                                                        >  2:1135982/1‑93 (MQ=255)

                                  gACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATaaa                                                       >  1:1570438/1‑93 (MQ=255)

                                  gACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCaa                                                       >  2:671214/1‑93 (MQ=255)

                                     gCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTa                                                    >  2:21235/1‑93 (MQ=255)

                                       aTTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATc                                                  <  2:2502420/93‑1 (MQ=255)

                                        tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 >  1:3118647/1‑93 (MQ=255)

                                        tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 <  1:389751/93‑1 (MQ=255)

                                          tAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGt                                               <  1:2948546/93‑1 (MQ=255)

                                           aaTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTg                                              >  1:4425852/1‑93 (MQ=255)

                                           aaTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTg                                              <  1:3568068/93‑1 (MQ=255)

                                            aTCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTgg                                             >  2:4305601/1‑93 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            <  2:121177/93‑1 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            <  1:3480651/93‑1 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            >  1:1304119/1‑93 (MQ=255)

                                             tCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGc                                            >  2:201490/1‑93 (MQ=255)

                                                gcgACCACGCCCGGAAATATGCGCACGGACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGt                                         <  1:1409555/93‑1 (MQ=255)

                                                  gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGTTGGAGAAATCAACTATCCGTGGCTGTTc                                       <  2:4408989/93‑1 (MQ=255)

                                                   aCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGGAATCAACTATCCGTGGCTGTTcc                                      >  2:2589184/1‑93 (MQ=255)

                                                     tACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGc                                    >  2:2799027/1‑93 (MQ=255)

                                                          cGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGtt                               >  2:4078119/1‑93 (MQ=255)

                                                          cGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGtt                               >  1:2087764/1‑93 (MQ=255)

                                                           gggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              >  1:4296686/1‑93 (MQ=255)

                                                            ggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCg                             >  2:629989/1‑93 (MQ=255)

                                                            ggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCg                             >  2:1244693/1‑93 (MQ=255)

                                                             gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            <  1:4345276/93‑1 (MQ=255)

                                                              aaaTATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGAt                           <  1:150969/93‑1 (MQ=255)

                                                                atatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGa                         <  1:2311645/93‑1 (MQ=255)

                                                                 tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        <  1:3949473/93‑1 (MQ=255)

                                                                   tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  1:525182/93‑1 (MQ=255)

                                                                   gtcgcgcgTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  2:2424481/91‑1 (MQ=255)

                                                                    gcgcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTg                     <  2:3286167/93‑1 (MQ=255)

                                                                     cgcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTgg                    <  1:4633541/93‑1 (MQ=255)

                                                                     cgcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTgg                    >  2:2264818/1‑93 (MQ=255)

                                                                     cgcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTgg                    >  2:4107216/1‑93 (MQ=255)

                                                                      gcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTggg                   >  2:2563133/1‑93 (MQ=255)

                                                                      gcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTggg                   >  1:19108/1‑93 (MQ=255)

                                                                       cACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGa                  >  1:4464561/1‑93 (MQ=255)

                                                                        aCGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGat                 <  1:1934043/93‑1 (MQ=255)

                                                                         cGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGata                <  2:1670433/93‑1 (MQ=255)

                                                                           tACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATc              >  2:2731497/1‑93 (MQ=255)

                                                                           tACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATc              >  1:3636427/1‑93 (MQ=255)

                                                                           tACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATc              <  1:3738149/93‑1 (MQ=255)

                                                                            aCAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCa             <  2:484854/93‑1 (MQ=255)

                                                                               aaTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAggg          >  2:4627435/1‑93 (MQ=255)

                                                                                aTTGCCGGACCACTGGATCGACATGAACCGGAAGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGc         <  1:1160095/93‑1 (MQ=255)

                                                                                aTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGGGGTGGGATATCAgggt         >  1:2798325/1‑92 (MQ=255)

                                                                                 ttGCCGTACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGGGGGATATCCGGGCt        >  2:3586828/1‑93 (MQ=255)

                                                                                 ttGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCt        <  2:33132/93‑1 (MQ=255)

                                                                                   gCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTgg      <  1:2191124/93‑1 (MQ=255)

                                                                                   gCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTgg      <  2:4658478/93‑1 (MQ=255)

                                                                                   gCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTgg      <  1:4668469/93‑1 (MQ=255)

                                                                                      ggACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATc   <  2:2934297/93‑1 (MQ=255)

                                                                                      ggACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATc   >  2:4714838/1‑93 (MQ=255)

                                                                                       gACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCg  <  2:738625/93‑1 (MQ=255)

ATACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862791‑2862970

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