BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,339,805	C→T	18.9%	Q239Q (CAG→CAA)	adiA ←	arginine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,339,805	0	C	T	18.9%	448.0 / 91.3	207	Q239Q (CAG→CAA)	adiA	arginine decarboxylase
Reads supporting (aligned to +/- strand): ref base C (91/77); new base T (23/16); total (114/93)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

Only 100 of 210 total aligned reads displayed.

CATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAA  >  NC_000913/4339706‑4339903

cATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGcc                                                                                                    >  1:1915544/1‑100 (MQ=255)

cATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGcc                                                                                                    <  2:310080/100‑1 (MQ=255)

 atatAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCt                                                                                                   >  2:1859916/1‑100 (MQ=255)

  tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTg                                                                                                  <  1:2892464/100‑1 (MQ=255)

  tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTg                                                                                                  <  1:2034938/100‑1 (MQ=255)

  tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGACTg                                                                                                  >  2:2850370/1‑100 (MQ=255)

   ataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGc                                                                                                 <  2:775902/100‑1 (MQ=255)

    taGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCa                                                                                                <  1:2187336/100‑1 (MQ=255)

     aGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCat                                                                                               <  1:2655289/100‑1 (MQ=255)

     aGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCat                                                                                               <  2:825815/100‑1 (MQ=255)

     aGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCa                                                                                                <  2:2171986/99‑1 (MQ=255)

      gACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatg                                                                                              >  1:2649243/1‑100 (MQ=255)

      gACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatg                                                                                              <  1:1116529/100‑1 (MQ=255)

      gACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatg                                                                                              >  2:2198562/1‑100 (MQ=255)

        ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCGGTTACGGTGAACGACCAGGACATCGTTATCGGTCATGCAGGCTTGCatgac                                                                                            <  2:1718521/100‑2 (MQ=255)

        ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGCCCACGGCATCGTTATCGGTCATGCAAGCCTGCatgat                                                                                            >  2:2573056/1‑100 (MQ=255)

        ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgac                                                                                            <  1:1026732/100‑2 (MQ=255)

        ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatgat                                                                                            >  1:2605477/1‑100 (MQ=255)

         cGGTGTCGCGCCTGTCAGCATCAAACCTGGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgacc                                                                                           <  2:555369/100‑3 (MQ=255)

          ggTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                          <  1:335178/100‑4 (MQ=255)

          ggTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                          <  1:2495781/100‑4 (MQ=255)

           gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccga                                                                                         <  2:769304/100‑5 (MQ=255)

           gTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGt                                                                                         >  2:645677/1‑100 (MQ=255)

           gTTGCGCGCCTGTCAGCATCAAACCTTGTTCGATGGCTGTATGGCAGTTACGGTCAACGTCCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                          <  2:236119/99‑4 (MQ=255)

             ttCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGc                                                                                       >  1:300802/1‑100 (MQ=255)

             ttCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGc                                                                                       <  2:2465962/100‑1 (MQ=255)

               cgcgCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                     <  2:940201/100‑1 (MQ=255)

               cgcgCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCgg                                                                                     >  1:1843835/1‑100 (MQ=255)

                gcgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGt                                                                                    <  1:2954121/100‑1 (MQ=255)

                 cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                   <  1:2803290/100‑1 (MQ=255)

                 cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                   >  1:425392/1‑100 (MQ=255)

                  gcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTa                                                                                  <  2:1528050/100‑1 (MQ=255)

                   ccTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTag                                                                                 <  2:2803976/100‑1 (MQ=255)

                    cTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTaga                                                                                <  1:142074/100‑1 (MQ=255)

                    cTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTaga                                                                                <  1:2210387/100‑1 (MQ=255)

                     tGTCAGCATCAAACCTTGTTCGATGGATTTATCGCAGTTACGGTCAACGACCACGACATCGTTACGGGTCATGCAAGCCTGCATGATCGTGAGGTTagag                                                                               >  2:1559190/1‑100 (MQ=255)

                        cAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCg                                                                            >  2:608345/1‑100 (MQ=255)

                            aTCAAACCTTGTTCGATGGATTTATGTCAGTTACGGTTAACGCCCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGaa                                                                         <  2:2305731/99‑1 (MQ=255)

                                aaCCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTAcc                                                                    <  1:3117651/100‑1 (MQ=255)

                                  ccTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcga                                                                  >  1:118068/1‑100 (MQ=255)

                                    ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                                <  2:1915533/100‑1 (MQ=255)

                                    ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                                <  2:2288797/100‑1 (MQ=255)

                                    ttGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgac                                                                 >  2:864546/1‑99 (MQ=255)

                                      gTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCgacgac                                                              <  2:349314/100‑1 (MQ=255)

                                         cGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAc                                                           >  2:2909186/1‑100 (MQ=255)

                                          gATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAcc                                                          >  2:900701/1‑100 (MQ=255)

                                          gATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAcc                                                          <  1:1799814/100‑1 (MQ=255)

                                          gATGGATATATGGCAGGTACGGTCAACGACCCCGACCTCGTTAGAGGTCCTGCAAGCCTGCATGATTGTGCGGCTAGAGCCGGAAGTACctacgcccac                                                           >  2:696986/1‑94 (MQ=255)

                                              gATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                      <  1:3097806/100‑1 (MQ=255)

                                              gATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                      >  2:1258007/1‑100 (MQ=255)

                                                   aTGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                                 >  1:87719/1‑100 (MQ=255)

                                                   aTGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAg                                                 <  1:141555/100‑1 (MQ=255)

                                                     ggCAGTTACGGTCAACGACCGCGCCATCGTTATCGGTCATGCAAGCCTGCATGACGGTGCGGTTAGAGCCGGAAGTCGCGACGACTCCGGACCAGGAg                                                 >  2:1325228/1‑98 (MQ=255)

                                                     ggCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCg                                               >  1:50244/1‑100 (MQ=255)

                                                        aGTTACGGTCAACGACCACGACATCGTTATCGGTCAGGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                            <  2:199979/100‑1 (MQ=255)

                                                         gTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCg                                           >  2:1001469/1‑100 (MQ=255)

                                                           tACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCgtc                                         <  1:2580251/100‑3 (MQ=255)

                                                           tACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGc                                         >  2:3130053/1‑100 (MQ=255)

                                                            aCGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGTCTACCGACCAGGAGCGATCGGc                                         <  2:1380049/99‑1 (MQ=255)

                                                             cGGTCAACGACTACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCgg                                          <  2:2604519/97‑1 (MQ=255)

                                                             cGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCAc                                       <  1:3568590/100‑1 (MQ=255)

                                                              ggTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCAcc                                      >  1:3789811/1‑100 (MQ=255)

                                                               gTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGCCCAGGAGCGATCGGCACCa                                     >  2:3717713/1‑100 (MQ=255)

                                                               gTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCa                                     >  1:2757656/1‑100 (MQ=255)

                                                                tCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCaa                                    >  1:3414125/1‑100 (MQ=255)

                                                                tCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCa                                     <  2:3208969/99‑1 (MQ=255)

                                                                 cAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCTGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCaaa                                   <  1:3640413/100‑1 (MQ=255)

                                                                     gACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcg                               <  2:1028544/100‑1 (MQ=255)

                                                                      aCCACGAGATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgc                              >  1:2004655/1‑100 (MQ=255)

                                                                       ccACGACATCGTTATCGGTCATGCAGGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgc                              >  2:1970142/1‑99 (MQ=255)

                                                                       ccACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATAcgcg                             <  1:3446846/100‑1 (MQ=255)

                                                                          cGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGCGCGATCGGCACCAAATACGCGTGc                          >  1:2222759/1‑100 (MQ=255)

                                                                          cGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGc                          <  2:1112618/100‑1 (MQ=255)

                                                                            aCATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGCGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCgg                        <  2:2446968/100‑1 (MQ=255)

                                                                             cATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGc                       >  1:1203217/1‑100 (MQ=255)

                                                                               aaccgcaccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCat                     >  1:1868535/11‑100 (MQ=38)

                                                                               aaccgcaccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCat                     >  1:242872/11‑100 (MQ=38)

                                                                                 gTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCatat                   <  2:3367428/100‑1 (MQ=255)

                                                                                 gTTATCGGTCATGCAAGCCTGAATGATGGTGCGGTTAGAGCCGGAAGTACCGTCGACTACCGACCAGGAGCGATCGGCACCAAATGCGCGTGCGGCatat                   >  2:2240745/1‑100 (MQ=255)

                                                                                      ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc              >  1:136777/4‑100 (MQ=255)

                                                                                      ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc              >  1:1487798/4‑100 (MQ=255)

                                                                                      cGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCTACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttc              >  2:2434264/1‑100 (MQ=255)

                                                                                       ggTCATTCAAGCCTGCATGATGGGGCGGTTAGAGTCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg             <  2:1488021/100‑1 (MQ=255)

                                                                                       ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg             <  1:556789/100‑1 (MQ=255)

                                                                                       ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg             <  1:209388/100‑1 (MQ=255)

                                                                                       ggTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg             <  1:672510/100‑1 (MQ=255)

                                                                                       catCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttca             >  1:93610/3‑99 (MQ=255)

                                                                                        atCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgc            >  1:2291188/2‑100 (MQ=255)

                                                                                        atCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgc            >  1:83830/2‑100 (MQ=255)

                                                                                           aTGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgcttt         >  1:1493876/1‑100 (MQ=255)

                                                                                           aTGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgcttt         >  1:1626212/1‑100 (MQ=255)

                                                                                            tGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttc        >  1:3476637/1‑100 (MQ=255)

                                                                                            tGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttc        >  2:2536113/1‑100 (MQ=255)

                                                                                             gCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcg       >  1:759572/1‑100 (MQ=255)

                                                                                             gCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcg       <  2:2024566/100‑1 (MQ=255)

                                                                                              cAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcgc      >  2:355651/1‑100 (MQ=255)

                                                                                              cAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcgc      <  2:1594280/100‑1 (MQ=255)

                                                                                                aGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCa    <  1:3191368/100‑1 (MQ=255)

                                                                                                 gCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCaa   >  2:302698/1‑100 (MQ=255)

                                                                                                  ccTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCaaa  >  2:1182933/1‑100 (MQ=255)

CATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAA  >  NC_000913/4339706‑4339903

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