BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,419,074	A→G	N145S (AAC→AGC)	yfcG →	putative S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,419,074	0	A	G	99.4%	531.1 / ‑6.6	165	N145S (AAC→AGC)	yfcG	putative S‑transferase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (81/83); total (82/83)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 171 total aligned reads displayed.

TCATCATTTTAATCACGCAG‑CCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAACAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTGGACTC  >  NC_000913/2418985‑2419163

catcattttaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCa                                                                                         >  2:1210557‑M2/22‑93 (MQ=255)

 catcatTTTAATCACGCAGCCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCaa                                                                                        <  1:1176288/93‑1 (MQ=255)

 atccttttaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCaa                                                                                        >  2:2475188‑M2/21‑93 (MQ=255)

 atcattttaatcacgcagccCCCACCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCaa                                                                                        <  2:2629919‑M2/73‑1 (MQ=255)

  tcattttaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAg                                                                                       <  2:3714371‑M2/74‑1 (MQ=255)

  atcatTTTAATCACGCAG‑CCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAACAAGc                                                                                      >  2:1973578/1‑93 (MQ=255)

    attttaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCg                                                                                     <  2:3811462‑M2/76‑1 (MQ=255)

    attttaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCg                                                                                     <  2:1207095‑M2/76‑1 (MQ=255)

      ttgaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGc                                                                                   <  2:3486681‑M2/78‑1 (MQ=255)

       ttaatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCt                                                                                  <  2:2053911‑M2/79‑1 (MQ=255)

        taatcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTg                                                                                 <  2:4129442‑M2/80‑1 (MQ=255)

          atcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGa                                                                               >  1:1037141‑M2/12‑93 (MQ=255)

           tcacgcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGaa                                                                              >  2:1900696‑M2/11‑93 (MQ=255)

           cacgcagcccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAATTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGaa                                                                              >  2:4483074‑M2/11‑93 (MQ=255)

               gcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACt                                                                          <  2:2010896‑M2/87‑1 (MQ=255)

               gcagccCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACt                                                                          >  1:2762988‑M2/7‑93 (MQ=255)

                  accCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGc                                                                       <  2:3852685‑M2/90‑1 (MQ=255)

                    cCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGccc                                                                     >  1:3341707‑M2/2‑93 (MQ=255)

                    cCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGccc                                                                     <  1:3650541‑M2/92‑1 (MQ=255)

                    cCCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGccc                                                                     <  1:3072078‑M2/92‑1 (MQ=255)

                     ccccccAAACCATTCCTTACGCTATTGATCGTTATCAGGTTGAACCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCct                                                                    <  1:4205472/93‑1 (MQ=255)

                     ccccccAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCct                                                                    >  1:2058071/1‑93 (MQ=255)

                      cccccAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCctg                                                                   >  2:2906505/1‑93 (MQ=255)

                       ccccAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCctgg                                                                  >  1:2754463/1‑93 (MQ=255)

                        cccAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCctggc                                                                 <  1:4506217/93‑1 (MQ=255)

                        cccAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCctggc                                                                 <  1:1247749/93‑1 (MQ=255)

                        cccAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACGCGCCctggc                                                                 >  1:2901131/1‑93 (MQ=255)

                          cAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCctggctg                                                               >  2:372879/1‑93 (MQ=255)

                          cAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCctggctg                                                               <  2:3674391/93‑1 (MQ=255)

                            aaCCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTggg                                                             >  1:4045570/1‑93 (MQ=255)

                            aaCCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAACGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTggg                                                             >  2:503884/1‑93 (MQ=255)

                            aaCCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAACGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTggg                                                             >  2:503885/1‑93 (MQ=255)

                             aCCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTTAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGa                                                            <  2:3185789/93‑1 (MQ=255)

                             aCCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGa                                                            <  2:784107/93‑1 (MQ=255)

                              ccATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAg                                                           <  2:3824992/93‑1 (MQ=255)

                                aTTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGc                                                         >  2:4274218/1‑93 (MQ=255)

                                aTTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGc                                                         >  1:301469/1‑93 (MQ=255)

                                aTTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGc                                                         <  1:296347/93‑1 (MQ=255)

                                 ttccttccGCTATTGAACGTTATCAGGTTGCAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCg                                                        <  1:652428/93‑1 (MQ=255)

                                 ttCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCg                                                        <  1:652429/93‑1 (MQ=255)

                                  tCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCga                                                       <  1:1668984/93‑1 (MQ=255)

                                   ccTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGAGGCTGGAAAACTCGCCCTGGCTGGGAGGCgag                                                      >  1:4640815/1‑93 (MQ=255)

                                    cTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCgaga                                                     <  2:4085661/93‑1 (MQ=255)

                                     ttACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGaa                                                    <  1:3641315/93‑1 (MQ=255)

                                      tACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAAc                                                   >  1:1233729/1‑93 (MQ=255)

                                      tACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAAc                                                   >  1:3134976/1‑93 (MQ=255)

                                         gCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTAc                                                <  2:1538711/93‑1 (MQ=255)

                                          cTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACa                                               <  1:3462071/93‑1 (MQ=255)

                                           tATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAg                                              >  1:1711766/1‑93 (MQ=255)

                                            aTTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTAAAGc                                             >  2:69776/1‑93 (MQ=255)

                                                aaCGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTg                                         <  2:1027433/93‑1 (MQ=255)

                                                aaCGATATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTg                                         >  1:164849/1‑93 (MQ=255)

                                                     tATCAGGTTGACACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGat                                    <  2:1981397/93‑1 (MQ=255)

                                                      aTCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGata                                   >  2:2272300/1‑93 (MQ=255)

                                                       tCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGctat                                  >  2:111433/1‑93 (MQ=255)

                                                        cAGGTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATAtt                                 <  2:3150415/93‑1 (MQ=255)

                                                         aGGTTGAACCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTg                                <  1:1655098/93‑1 (MQ=255)

                                                         aGGTTGAAACTCTGCGTCTTTACAATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGataata                                <  2:4216635/93‑4 (MQ=255)

                                                          ggTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGGGGATATTGc                               >  1:1439620/1‑93 (MQ=255)

                                                          ggTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGc                               >  1:2619530/1‑93 (MQ=255)

                                                           gTTGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGcc                              >  1:27220/1‑93 (MQ=255)

                                                            ttGAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCct                             <  1:2800283/93‑1 (MQ=255)

                                                              gAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgc                           >  2:3396068/1‑93 (MQ=255)

                                                              gAAACTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgc                           >  2:3326181/1‑93 (MQ=255)

                                                               aaaCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgct                          <  2:1102923/93‑1 (MQ=255)

                                                                aaCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgctg                         <  1:1592799/93‑1 (MQ=255)

                                                                aaCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgctg                         >  2:1734839/1‑93 (MQ=255)

                                                                aaCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgctg                         <  1:4376583/93‑1 (MQ=255)

                                                                aaCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCctgctg                         >  1:559/1‑93 (MQ=255)

                                                                 aCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTgg                        >  2:3114943/1‑93 (MQ=255)

                                                                 aCTCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTgg                        >  2:901790/1‑93 (MQ=255)

                                                                   tCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGcc                      <  2:151574/93‑1 (MQ=255)

                                                                   tCAGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGcc                      <  2:40288/93‑1 (MQ=255)

                                                                     aGCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGt                    <  1:3992420/93‑1 (MQ=255)

                                                                     aGCGTCTTTACCAGGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGt                    <  2:528127/93‑1 (MQ=255)

                                                                      gCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTg                   <  2:1502091/93‑1 (MQ=255)

                                                                      gCGTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTg                   <  2:1253967/93‑1 (MQ=255)

                                                                        gTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTggg                 >  2:804827/1‑93 (MQ=255)

                                                                        gTCTTTACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTggg                 >  1:3397478/1‑93 (MQ=255)

                                                                           tttACCATGTCCTGAGCAAGCGGCTGGAAAACCCGCCCTGGCTGGGAGGCGCGACCCACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTa              >  2:4130215/1‑93 (MQ=255)

                                                                           tttACCATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTa              >  1:1009009/1‑93 (MQ=255)

                                                                           tttACCATCTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTa              <  1:1746188/93‑1 (MQ=255)

                                                                               ccATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGc          >  1:1329069/1‑93 (MQ=255)

                                                                               acATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGc          <  1:745080/92‑1 (MQ=255)

                                                                                cATGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGcc         >  1:3475786/1‑93 (MQ=255)

                                                                                 aTGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCt        >  2:3267738/1‑93 (MQ=255)

                                                                                 aTGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCt        >  1:1961076/1‑93 (MQ=255)

                                                                                 aTGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCt        >  1:2557049/1‑93 (MQ=255)

                                                                                 aTGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCt        <  1:232902/93‑1 (MQ=255)

                                                                                  tGTACTGAGCAAGCGGCTGGAACACTCGCCCTAGCTGGGAGGCGAGAACTACAGCATTGCGTATATTGCCTGCTGGCCGTGGGTTAATGCCTg       >  2:1597845/1‑93 (MQ=255)

                                                                                  tGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTg       <  2:2683958/93‑1 (MQ=255)

                                                                                  tGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTg       >  2:1597844/1‑93 (MQ=255)

                                                                                  tGTACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTg       <  2:431218/93‑1 (MQ=255)

                                                                                   gTACTGAGCAAGCGGCTGGAAAACTCGCCCTGACTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTgg      <  2:1149304/93‑1 (MQ=255)

                                                                                   gTACTGAACAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTgg      <  1:3279705/93‑1 (MQ=255)

                                                                                    tACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTggc     >  1:2832541/1‑92 (MQ=255)

                                                                                    tACTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTGGa     >  2:4252615/1‑93 (MQ=255)

                                                                                     aCTGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTGGAc    <  1:383459/93‑1 (MQ=255)

                                                                                       tGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTGGACTc  >  1:2245258/1‑93 (MQ=255)

                                                                                       tGAGCAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTGGACTc  >  2:1464170/1‑93 (MQ=255)

TCATCATTTTAATCACGCAG‑CCCCCCAAACCATTCCTTACGCTATTGAACGTTATCAGGTTGAAACTCAGCGTCTTTACCATGTACTGAACAAGCGGCTGGAAAACTCGCCCTGGCTGGGAGGCGAGAACTACAGCATTGCGGATATTGCCTGCTGGCCGTGGGTTAATGCCTGGACTC  >  NC_000913/2418985‑2419163

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