BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	138,888	T→C	T780A (ACG→GCG)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	138,888	0	T	C	99.7%	1039.6 / NA	302	T780A (ACG→GCG)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (171/130); minor base A (0/1); total (171/131)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

TGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTG  >  NC_000913/138798‑138977

tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCAc                                                                                         <  1:2037596/93‑1 (MQ=255)

tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCAc                                                                                         >  1:6063437/1‑93 (MQ=255)

tGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCAc                                                                                         >  1:2390662/1‑93 (MQ=255)

   acaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCaa                                                                                      >  2:6898915/1‑93 (MQ=255)

   acaacGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCaa                                                                                      >  1:547416/1‑93 (MQ=255)

      acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATg                                                                                   >  1:3233204/1‑93 (MQ=255)

      acGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATg                                                                                   <  2:5839329/93‑1 (MQ=255)

        gcgcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGaa                                                                                 <  1:7867078/93‑1 (MQ=255)

         cgcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAAc                                                                                <  2:1176445/93‑1 (MQ=255)

          gcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAAcc                                                                               >  2:3534016/1‑93 (MQ=255)

          gcACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCAACAAATGAAcc                                                                               >  2:7150064/1‑93 (MQ=255)

           cACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCg                                                                              <  2:2162315/93‑1 (MQ=255)

            aCGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGt                                                                             >  2:3919605/1‑93 (MQ=255)

                   cAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGcc                                                                      >  2:5634229/1‑93 (MQ=255)

                   cAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGcc                                                                      <  1:4145958/93‑1 (MQ=255)

                    aGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCt                                                                     >  2:2978000/1‑93 (MQ=255)

                    aGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCt                                                                     >  2:374745/1‑93 (MQ=255)

                     gCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTg                                                                    <  1:6606143/93‑1 (MQ=255)

                     gCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTg                                                                    >  2:4441167/1‑93 (MQ=255)

                     gCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTg                                                                    <  1:5055180/93‑1 (MQ=255)

                      cAAACTGAGCAAGTCTTACTTCACATCATCCTGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACAAAATGAACCGTGACCGCCTGc                                                                   <  1:1523785/93‑1 (MQ=255)

                       aaTCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCATCAAATGAACCGTGACCGCCTGCg                                                                  >  2:5993868/1‑93 (MQ=255)

                       aaTCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCg                                                                  >  1:6900894/1‑93 (MQ=255)

                       aaTCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCg                                                                  >  2:776249/1‑93 (MQ=255)

                       aaTCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCg                                                                  <  1:8171097/93‑1 (MQ=255)

                         tCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGga                                                                >  2:2462585/1‑93 (MQ=255)

                         tCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGga                                                                >  1:7083955/1‑93 (MQ=255)

                         tCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGga                                                                >  1:4797464/1‑93 (MQ=255)

                          cTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGgag                                                               >  1:5262716/1‑93 (MQ=255)

                           tGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGgaga                                                              >  1:5177506/1‑93 (MQ=255)

                           tGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGgaga                                                              <  2:1653444/93‑1 (MQ=255)

                           tGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGgaga                                                              <  1:4350308/93‑1 (MQ=255)

                            taGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGAt                                                             <  1:6432747/92‑1 (MQ=255)

                              gCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATca                                                           <  1:61383/93‑1 (MQ=255)

                                aaGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacc                                                         >  2:7601199/1‑93 (MQ=255)

                                 aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                        <  2:8358738/93‑1 (MQ=255)

                                 aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                        <  1:8530250/93‑1 (MQ=255)

                                 aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                        >  1:4075156/1‑93 (MQ=255)

                                 aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                        >  1:2948980/1‑93 (MQ=255)

                                 aGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcacca                                                        >  1:4914615/1‑93 (MQ=255)

                                  gTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATcaccac                                                       >  2:1964271/1‑93 (MQ=255)

                                    cTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACAt                                                     <  1:7780509/93‑1 (MQ=255)

                                     ttACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATa                                                    <  1:1935845/93‑1 (MQ=255)

                                      tACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATAc                                                   <  1:8746465/93‑1 (MQ=255)

                                      tACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATAc                                                   >  1:7645145/1‑93 (MQ=255)

                                      tACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATAc                                                   >  1:753623/1‑93 (MQ=255)

                                        cTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTg                                                 <  2:2862166/93‑1 (MQ=255)

                                        cTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTg                                                 <  2:8160469/93‑1 (MQ=255)

                                          tCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCt                                               <  2:7705411/93‑1 (MQ=255)

                                          tCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCt                                               <  2:4164047/93‑1 (MQ=255)

                                           cacaTCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCtt                                              >  2:5446318/1‑93 (MQ=255)

                                             catcatCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTAc                                            >  2:1818433/1‑93 (MQ=255)

                                             catcatCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTAc                                            >  1:3613671/1‑93 (MQ=255)

                                              atcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTAcc                                           <  1:4123901/93‑1 (MQ=255)

                                               tcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCa                                          >  2:5326901/1‑93 (MQ=255)

                                               tcatcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCa                                          <  1:3397022/93‑1 (MQ=255)

                                                catcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTT‑CCAtt                                        >  2:3200305/1‑93 (MQ=255)

                                                 atcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCAtt                                        >  2:4996011/1‑93 (MQ=255)

                                                 atcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCAtt                                        >  1:6959216/1‑93 (MQ=255)

                                                 atcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCAtt                                        <  1:4249790/93‑1 (MQ=255)

                                                 atcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCAtt                                        >  1:2197812/1‑93 (MQ=255)

                                                  tcCGGCAGCGCATAAGCCACAATATAGTCTCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTc                                       <  2:6495530/93‑1 (MQ=255)

                                                  tcCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTc                                       >  1:1035209/1‑93 (MQ=255)

                                                    cGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCAc                                     >  1:3957987/1‑93 (MQ=255)

                                                     ggCAGCGTATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACTGTTACCGCCTGCGGAGATCACCATATATTGCTTACCTTTCACt                                    >  2:7838123/1‑93 (MQ=255)

                                                        aGCGCATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCa                                 <  1:5865918/93‑1 (MQ=255)

                                                         gcgcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCAt                                >  1:1217552/1‑93 (MQ=255)

                                                         gcgcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCAt                                >  1:7280048/1‑93 (MQ=255)

                                                          cgcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATa                               <  1:4654162/93‑1 (MQ=255)

                                                          cgcAAAAGCCACAAAATAGTCGCCCATCTTCCCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATa                               <  2:7464356/93‑1 (MQ=255)

                                                           gcATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAg                              >  1:6645120/1‑93 (MQ=255)

                                                            cATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAgg                             <  2:4597803/93‑1 (MQ=255)

                                                            cATAAGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAgg                             <  2:7362469/93‑1 (MQ=255)

                                                               aaGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCa                          >  2:708675/1‑93 (MQ=255)

                                                                aGCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCAt                         <  1:2245221/93‑1 (MQ=255)

                                                                 gCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCAtt                        >  1:946900/1‑93 (MQ=255)

                                                                 gCCACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCAtt                        <  1:2753536/93‑1 (MQ=255)

                                                                  ccACAATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTg                       <  1:5119190/93‑1 (MQ=255)

                                                                    acaATATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGc                     <  2:3630144/93‑1 (MQ=255)

                                                                      aaTATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGt                   <  1:5210674/93‑1 (MQ=255)

                                                                      aaTATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGt                   >  1:2401879/1‑93 (MQ=255)

                                                                      aaTATAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGt                   >  2:6998153/1‑93 (MQ=255)

                                                                       atatAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGTGTa                  <  2:8629508/93‑1 (MQ=255)

                                                                       atatAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTa                  >  2:4771181/1‑93 (MQ=255)

                                                                       atatAGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTa                  >  1:7361078/1‑93 (MQ=255)

                                                                        tataGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAg                 >  1:8384498/1‑93 (MQ=255)

                                                                        tataGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAg                 <  2:2605034/93‑1 (MQ=255)

                                                                         ataGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGc                <  2:5378211/93‑1 (MQ=255)

                                                                         ataGTCGCCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGc                <  1:5782577/93‑1 (MQ=255)

                                                                           aGTCACCCATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCt              >  2:142447/1‑93 (MQ=255)

                                                                                  cATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGaccacc       >  1:4720549/1‑93 (MQ=255)

                                                                                  cATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGaccacc       >  2:524227/1‑93 (MQ=255)

                                                                                  cATCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGACTGcacacc       >  1:1668550/1‑93 (MQ=255)

                                                                                   aTCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCAccc      >  2:4533525/1‑93 (MQ=255)

                                                                                   aTCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCAccc      >  1:1844319/1‑93 (MQ=255)

                                                                                    tCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCAcccc     >  2:6373347/1‑92 (MQ=255)

                                                                                    tCTTCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCg     <  1:6447003/93‑1 (MQ=255)

                                                                                      ttCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCt   >  2:6742165/1‑93 (MQ=255)

                                                                                      ttCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCt   >  1:2456843/1‑93 (MQ=255)

                                                                                      ttCGCACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCt   >  1:8457733/1‑93 (MQ=255)

TGAACAACGCGCACGTTGTCAGCAATCTGAGCAAGTCTTACTTCACATCATCCGGCAGCGCATAAGCCACAATATAGTCGCCCATCTTCGTACCAAATGAACCGTGACCGCCTGCGGAGATCACCACATACTGCTTACCATTCACTTCATAGGTCATTGGCGTAGCCTGACCACCCGCTG  >  NC_000913/138798‑138977

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