BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,431,142	T→C	H154R (CAC→CGC)	ytfG ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,431,142	0	T	C	99.5%	669.0 / ‑4.8	216	H154R (CAC→CGC)	ytfG	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (108/107); total (109/107)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 239 total aligned reads displayed.

TCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAG  >  NC_000913/4431052‑4431230

  gCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTc                                                                                      <  2:1247711/93‑1 (MQ=255)

   cTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTcc                                                                                     <  1:2816920/93‑1 (MQ=255)

    tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                    <  2:695942/93‑1 (MQ=255)

    tAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCa                                                                                    >  2:440191/1‑93 (MQ=255)

     aaTCCCGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAg                                                                                   <  2:660879/93‑1 (MQ=255)

      aTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGt                                                                                  >  2:3856674/1‑93 (MQ=255)

      aTCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGt                                                                                  >  1:3438905/1‑93 (MQ=255)

       tCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtg                                                                                 >  2:230301/1‑93 (MQ=255)

       tCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtg                                                                                 <  1:1997762/93‑1 (MQ=255)

         aCGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGtgct                                                                               >  1:2222560/1‑93 (MQ=255)

             gTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCg                                                                           <  2:430900/93‑1 (MQ=255)

             gTGCCGCAGCTGCCGCATAATCTGCCCGCGGTGCTGACGCGATTTTGCCATCGGCCGCCGCACCGATAAATACGCAGCGTTCCAGTGCTGCCg                                                                           >  2:3774511/1‑93 (MQ=255)

                ccGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCgggg                                                                        <  2:3361078/93‑1 (MQ=255)

                   caccTGCCGCATAATCGGCCCGCGTGGCTGAGGCGATTTTGCCATGGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcgc                                                                     <  2:2396276/90‑1 (MQ=255)

                   cAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGgcgc                                                                     <  1:1398908/93‑1 (MQ=255)

                       tGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGc                                                                 >  1:100927/1‑93 (MQ=255)

                          cGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAg                                                              <  2:3715283/93‑1 (MQ=255)

                          cGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAg                                                              >  1:3143208/1‑93 (MQ=255)

                           gCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAgg                                                             >  1:4200506/1‑93 (MQ=255)

                           gCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAgg                                                             >  2:2979853/1‑93 (MQ=255)

                            cATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGt                                                            >  1:2113225/1‑93 (MQ=255)

                            cATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGt                                                            >  1:4056626/1‑93 (MQ=255)

                             aTAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGGGGTa                                                           >  2:128745/1‑93 (MQ=255)

                             aTAATCGGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTa                                                           <  1:2883714/93‑1 (MQ=255)

                              tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAg                                                          >  2:1049287/1‑93 (MQ=255)

                              tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAg                                                          >  1:1864562/1‑93 (MQ=255)

                              tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGGAg                                                          >  1:1444726/1‑93 (MQ=255)

                              tAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCCGCCGGGGCGCTGGCGAGGTAg                                                          >  1:3532514/1‑93 (MQ=255)

                               aaTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGt                                                         >  1:2675819/1‑93 (MQ=255)

                               aaTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGt                                                         <  1:1213247/93‑1 (MQ=255)

                               aaTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGt                                                         >  2:1140576/1‑93 (MQ=255)

                                aTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGtt                                                        <  1:348106/93‑1 (MQ=255)

                                aTCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGtt                                                        >  2:2499824/1‑93 (MQ=255)

                                 tCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGttt                                                       <  1:4415435/93‑1 (MQ=255)

                                  cTGCCCGCGTTGCTGAGGCGATTTTGCTATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGtttt                                                      <  2:1604214/93‑1 (MQ=255)

                                  cTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGtttt                                                      <  2:1312718/93‑1 (MQ=255)

                                   tGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTc                                                     <  1:2435091/93‑1 (MQ=255)

                                        gcgTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGt                                                <  1:1349229/93‑1 (MQ=255)

                                          gTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAc                                              >  2:3016676/1‑93 (MQ=255)

                                           ttGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTTGCGAGGGAGTTTTCGCTGTAcc                                             >  2:2109777/1‑93 (MQ=255)

                                           tggCTGAGGCGATTTTGCCACCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTAcc                                             <  2:4158511/91‑1 (MQ=255)

                                             gCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAg                                           <  1:574341/93‑1 (MQ=255)

                                             gCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAg                                           >  1:3893038/1‑93 (MQ=255)

                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                        >  2:3305573/1‑93 (MQ=255)

                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                        >  2:2905222/1‑93 (MQ=255)

                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                        >  2:901787/1‑93 (MQ=255)

                                                gAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGGGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCg                                        >  2:176216/1‑93 (MQ=255)

                                                 aGGCTATTTTGCCCTCGCCCCCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCGGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGt                                       <  2:2792033/93‑1 (MQ=255)

                                                 aGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGt                                       >  1:2962083/1‑93 (MQ=255)

                                                 aGGCGATTTGGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGt                                       <  2:4075857/93‑1 (MQ=255)

                                                   gCGTTTTTCCCATCCCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTg                                     <  2:2645504/93‑1 (MQ=255)

                                                    gaaTTTTGCCACCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgc                                    <  2:1724250/91‑1 (MQ=255)

                                                    caaTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgc                                    <  2:3035993/91‑1 (MQ=255)

                                                    cGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgc                                    <  2:3792514/93‑1 (MQ=255)

                                                    cGATTTTGCCATCCCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgc                                    <  1:3953826/93‑1 (MQ=255)

                                                     gATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcg                                   >  1:3051909/1‑93 (MQ=255)

                                                      aTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgc                                  >  1:3758657/1‑93 (MQ=255)

                                                      aTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgc                                  >  1:2482774/1‑93 (MQ=255)

                                                      aTTTTGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgc                                  <  1:3819582/93‑1 (MQ=255)

                                                       ttttGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTgcgcg                                 >  1:778306/1‑92 (MQ=255)

                                                         ttGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGAGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagc                               >  1:4264181/1‑93 (MQ=255)

                                                          tcccATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagca                              <  2:454933/91‑1 (MQ=255)

                                                          tGCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCgcagca                              <  1:3155764/93‑1 (MQ=255)

                                                           gcaaTCCCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGcagcag                             <  2:1413861/90‑1 (MQ=255)

                                                           gCCATCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGcagcag                             <  1:861736/93‑1 (MQ=255)

                                                            ccAGCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAgg                            <  1:2053960/93‑1 (MQ=255)

                                                              aTCGCCCGCCGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGgt                          <  2:3193185/93‑1 (MQ=255)

                                                                 gcccgccGCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTa                       >  1:112101/1‑93 (MQ=255)

                                                                   ccgccgCACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTaaa                     <  2:3564994/93‑1 (MQ=255)

                                                                    cgccgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAAc                    >  2:3407225/1‑93 (MQ=255)

                                                                      ccgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGa                  <  1:4182542/93‑1 (MQ=255)

                                                                      ccgcACCGATAAATACGCCGCGTTCCAGGGTTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGa                  <  2:2008900/93‑1 (MQ=255)

                                                                       cgcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGAt                 <  2:48132/93‑1 (MQ=255)

                                                                        gcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATg                >  1:3386688/1‑93 (MQ=255)

                                                                        gcACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATg                >  1:1684309/1‑93 (MQ=255)

                                                                         cACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc               <  2:4240058/93‑1 (MQ=255)

                                                                         cACCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc               <  1:2242209/93‑1 (MQ=255)

                                                                         cACCGATAAATACCCCCCGTTCCATTCCTGCCGGTGCGCTGGCGAGGTAGTTTTCGCTCTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGc               <  2:3599431/93‑1 (MQ=255)

                                                                          aCCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGCTGcc              >  2:2286523/1‑93 (MQ=255)

                                                                          aCCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGcc              >  1:3377568/1‑93 (MQ=255)

                                                                          aCCGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGcc              >  1:547740/1‑93 (MQ=255)

                                                                           ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa             >  2:206530/1‑93 (MQ=255)

                                                                           ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa             <  1:2628428/93‑1 (MQ=255)

                                                                           ccGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCa             >  1:3228924/1‑93 (MQ=255)

                                                                            cGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAg            >  2:74312/1‑93 (MQ=255)

                                                                            cGATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAg            >  1:2857650/1‑93 (MQ=255)

                                                                             taTAAATACGCCGCGTTCCAGTGCTTCCGGGTCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGa           <  2:3623918/92‑1 (MQ=255)

                                                                             gATAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCAAGa           >  1:2542769/1‑93 (MQ=255)

                                                                              aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa          >  2:3798447/1‑93 (MQ=255)

                                                                              aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa          >  2:786768/1‑93 (MQ=255)

                                                                              aTAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGaa          <  1:4579374/93‑1 (MQ=255)

                                                                               tAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAt         <  2:821196/93‑1 (MQ=255)

                                                                               tAAATACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAAt         >  2:893954/1‑93 (MQ=255)

                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGTCGAGGTAGTTTTCGCTGTGCCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc        <  1:1608106/93‑1 (MQ=255)

                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc        <  2:4624117/93‑1 (MQ=255)

                                                                                aaaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATc        <  2:3239717/93‑1 (MQ=255)

                                                                                 aaTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATca       <  1:1172757/93‑1 (MQ=255)

                                                                                  aTACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcag      <  2:112239/93‑1 (MQ=255)

                                                                                   tACGCCGCGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATgagc     >  1:1067264/1‑93 (MQ=255)

                                                                                     cgccgcGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATcagcca   >  1:4008590/1‑93 (MQ=255)

TCGCTAATCACGCGTGCCGCAGCTGCCGCATAATCTGCCCGCGTTGCTGAGGCGATTTTGCCATCGCCCGCCGCACCGATAAATACGCCGTGTTCCAGTGCTGCCGGGGCGCTGGCGAGGTAGTTTTCGCTGTACCAGCCGTTGCGCAGCAGGGTGTAAACGATGCCAGAATCAGCCAG  >  NC_000913/4431052‑4431230

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