BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,587	T→C	S116G (AGT→GGT)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,587	0	T	C	99.2%	397.5 / ‑4.5	126	S116G (AGT→GGT)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (60/65); total (60/66)
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.50e-01

Only 100 of 128 total aligned reads displayed.

GCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAT  >  NC_000913/1616497‑1616674

gCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCgg                                                                                       <  1:90258/93‑1 (MQ=255)

  gcTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCTCCGGGt                                                                                     <  2:814385/92‑1 (MQ=255)

  aCTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGt                                                                                     <  1:4261475/93‑1 (MQ=255)

  aCTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGt                                                                                     >  2:4327318/1‑93 (MQ=255)

    tgtgGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTc                                                                                   <  1:4672749/93‑1 (MQ=255)

      tgGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCAt                                                                                 >  1:2080135/1‑93 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                <  2:3703227/93‑1 (MQ=255)

       ggACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATc                                                                                >  2:1696696/1‑93 (MQ=255)

           gCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTcc                                                                            <  2:4301201/93‑1 (MQ=255)

           gCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTcc                                                                            <  2:54944/93‑1 (MQ=255)

            ccGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCa                                                                           >  2:330343/1‑93 (MQ=255)

             cGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAg                                                                          >  1:1052479/1‑93 (MQ=255)

              gAGCTAATAATCATCGCTATGGTTCGCGGACCGGCAGTGCTTGGCAGTGCCGGCGGCACAAAGGCAATATTGGCACGGGGTTCATCTTCCAGc                                                                         >  1:2263001/1‑93 (MQ=255)

              gAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGc                                                                         >  1:2487598/1‑93 (MQ=255)

              gAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGc                                                                         <  1:1194587/93‑1 (MQ=255)

               aGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGct                                                                        >  1:3035480/1‑93 (MQ=255)

               aGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTAATCTTCCAggt                                                                        >  1:1611547/1‑91 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                       >  1:736727/1‑93 (MQ=255)

                gCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctc                                                                       >  2:3477647/1‑93 (MQ=255)

                 cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctct                                                                      >  1:2471574/1‑93 (MQ=255)

                 cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctct                                                                      >  1:1101794/1‑93 (MQ=255)

                 cTAATAATCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGctct                                                                      <  2:4310183/93‑1 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     >  2:4696650/1‑93 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     <  2:2208273/93‑1 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     >  2:1558444/1‑93 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCtt                                                                     >  2:2123300/1‑93 (MQ=255)

                  taataaTCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGCTCATCTTCCAACTCtt                                                                     >  1:4596428/1‑93 (MQ=255)

                   aataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                    >  2:2570506/1‑93 (MQ=255)

                   aataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                    >  2:577983/1‑93 (MQ=255)

                   aataatCATCCCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTc                                                                    <  1:2300350/93‑1 (MQ=255)

                    ataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTcc                                                                   >  1:3894967/1‑93 (MQ=255)

                    ataatCATCGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTcc                                                                   <  1:3662806/93‑1 (MQ=255)

                       atcatcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGAc                                                                >  2:3843997/1‑93 (MQ=255)

                         catcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACtt                                                              <  2:1886456/93‑1 (MQ=255)

                         catcGCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACtt                                                              <  1:222241/93‑1 (MQ=255)

                             gCTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtt                                                          >  1:410965/1‑93 (MQ=255)

                             gCTATGATTCCCGGACCGGAAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCtt                                                          <  1:901150/93‑1 (MQ=255)

                              cTATGGTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCttt                                                         <  1:474475/93‑1 (MQ=255)

                                aTGGTTCCCGGACCGGCAGTGTTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGCATTGCTTGTc                                                       >  1:4696337/1‑93 (MQ=255)

                                  ggTTCCCGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGc                                                     >  1:1385605/1‑93 (MQ=255)

                                       ccGGACCGGCAGTGCTTGGCATTGCCAGCGTCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCg                                                <  2:3281877/93‑1 (MQ=255)

                                       ccGGACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCg                                                <  1:1783073/93‑1 (MQ=255)

                                        cGGACCGGCAGTGCTTGTCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCgg                                               <  2:1786532/93‑1 (MQ=255)

                                         ggACCGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGt                                              >  1:3338593/1‑93 (MQ=255)

                                            ccGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaa                                           <  2:3632791/93‑1 (MQ=255)

                                            ccGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaa                                           <  2:3632790/93‑1 (MQ=255)

                                             cGGCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaat                                          >  1:2123994/1‑93 (MQ=255)

                                              ggCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatc                                         >  1:4691020/1‑93 (MQ=255)

                                              ggCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatc                                         >  2:2511498/1‑93 (MQ=255)

                                              ggCAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaata                                         >  2:2511497/1‑92 (MQ=255)

                                                cAGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaataaa                                       >  2:3524835/1‑93 (MQ=255)

                                                 aGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaat                                      <  2:2780738/93‑1 (MQ=255)

                                                 aGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaat                                      >  2:4564271/1‑93 (MQ=255)

                                                 aGTGCTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaat                                      <  2:1716157/93‑1 (MQ=255)

                                                  gtgcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGaatcaatc                                     >  1:4544031/4‑93 (MQ=255)

                                                   tgcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCg                                    >  2:3583530/3‑93 (MQ=255)

                                                    gcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGc                                   <  2:3574338/92‑1 (MQ=255)

                                                    gcTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGc                                   <  1:4235361/92‑1 (MQ=255)

                                                      ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                 >  2:3974071/1‑93 (MQ=255)

                                                      ttGGCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCtt                                 >  1:2729212/1‑93 (MQ=255)

                                                        ggCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                               <  1:1154478/93‑1 (MQ=255)

                                                        ggCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTc                               <  1:661098/93‑1 (MQ=255)

                                                         gCATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCt                              <  1:2639656/93‑1 (MQ=255)

                                                          cATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTg                             <  1:2810425/93‑1 (MQ=255)

                                                          cATTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTg                             <  1:505813/93‑1 (MQ=255)

                                                           aTTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGt                            >  2:3068625/1‑93 (MQ=255)

                                                           aTTGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGt                            >  2:1663213/1‑93 (MQ=255)

                                                            ttGCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGtt                           >  1:415298/1‑93 (MQ=255)

                                                              tccAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGc                         <  2:2843110/92‑1 (MQ=255)

                                                              gCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGc                         <  2:4617616/93‑1 (MQ=255)

                                                              gCCAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGc                         <  1:807866/93‑1 (MQ=255)

                                                                cAGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCgg                       >  2:1117262/1‑93 (MQ=255)

                                                                 aGCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGa                      >  2:3224752/1‑93 (MQ=255)

                                                                  gCGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaa                     >  1:4687533/1‑93 (MQ=255)

                                                                   cGGCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGaaa                    <  1:1313971/93‑1 (MQ=255)

                                                                    ggCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAag                   >  1:1887501/1‑93 (MQ=255)

                                                                    ggCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAag                   >  2:3137062/1‑93 (MQ=255)

                                                                     gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                  <  2:1729357/93‑1 (MQ=255)

                                                                     gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                  >  1:3963909/1‑93 (MQ=255)

                                                                     gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                  <  1:4401599/93‑1 (MQ=255)

                                                                     gCACAAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAaga                  <  1:7647/93‑1 (MQ=255)

                                                                      cacaAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGATTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAagag                 <  1:4598041/93‑1 (MQ=255)

                                                                        caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCa               <  1:4461691/93‑1 (MQ=255)

                                                                        caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCa               >  2:4641186/1‑93 (MQ=255)

                                                                        caAAGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGCGAATCAATCGCTTTCTGTTGCGGAAAGAGCa               >  1:1083069/1‑93 (MQ=255)

                                                                         aaaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCAt              >  1:3225647/1‑93 (MQ=255)

                                                                         aaaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCAt              >  1:670038/1‑93 (MQ=255)

                                                                          aaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc             >  2:2291773/1‑93 (MQ=255)

                                                                          aaGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATc             >  2:1296232/1‑93 (MQ=255)

                                                                           aGGCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATcc            <  2:3995376/93‑1 (MQ=255)

                                                                            ggCAATATTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCg           <  2:1616148/93‑1 (MQ=255)

                                                                                atatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaaa       >  2:3851835/1‑93 (MQ=255)

                                                                                atatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGaaaa       <  1:1148001/93‑1 (MQ=255)

                                                                                 tatTGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAAc      <  1:910883/93‑1 (MQ=255)

                                                                                   ttGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCg    <  2:2908158/93‑1 (MQ=255)

                                                                                    tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa   <  2:3104381/93‑1 (MQ=255)

                                                                                    tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa   <  2:3038450/93‑1 (MQ=255)

                                                                                    tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa   <  1:386375/93‑1 (MQ=255)

                                                                                    tGGCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGa   <  1:4404135/93‑1 (MQ=255)

                                                                                     ggCACCGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAt  >  1:3850723/1‑93 (MQ=255)

GCACTGTGGACGCCGAGCTAATAATCATCGCTATGGTTCCCGGACCGGCAGTACTTGGCATTGCCAGCGGCACAAAGGCAATATTGGCACTGGGTTCATCTTCCAGCTCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGAT  >  NC_000913/1616497‑1616674

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