BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,856,333	T→G	17.6%	E199A (GAG→GCG) ‡	ydjH ←	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,856,333	0	T	G	17.6%	263.3 / 58.6	121	E199A (GAG→GCG)	ydjH	putative kinase
Reads supporting (aligned to +/- strand): ref base T (48/51); new base G (13/8); total (61/60)
Fisher's exact test for biased strand distribution p-value = 3.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

Only 100 of 122 total aligned reads displayed.

CCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATA  >  NC_000913/1856234‑1856431

ccGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCt                                                                                                    >  2:3599792/1‑100 (MQ=255)

 cGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTc                                                                                                   >  2:404215/1‑100 (MQ=255)

  gTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCg                                                                                                  >  2:3448694/1‑100 (MQ=255)

   tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCgg                                                                                                 <  1:416662/100‑1 (MQ=255)

    cTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGc                                                                                                <  2:3283855/100‑1 (MQ=255)

       tACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCTGCaaa                                                                                             >  1:1075320/1‑100 (MQ=255)

         cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCTGCaaaa                                                                                            >  1:1307047/1‑99 (MQ=255)

         cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                           >  1:1357779/1‑100 (MQ=255)

         cccGTTTCAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                           <  2:896135/100‑1 (MQ=255)

          ccGTTTTAATCGCCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTGCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAtt                                                                                          <  2:1442725/100‑1 (MQ=255)

          ccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAtt                                                                                          <  1:1402946/100‑1 (MQ=255)

          ccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAtt                                                                                          <  1:3265883/100‑1 (MQ=255)

           cGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTa                                                                                         <  2:207248/100‑1 (MQ=255)

           cGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTa                                                                                         <  2:1822432/100‑1 (MQ=255)

            gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAg                                                                                        >  2:1886709/1‑100 (MQ=255)

            gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                        <  1:3531786/100‑2 (MQ=255)

            gTTTTAATCACCACCG‑TTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgg                                                                                       >  2:950784/1‑100 (MQ=255)

             ttttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgg                                                                                       <  1:93082/100‑1 (MQ=255)

              tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgaa                                                                                      <  1:2451422/100‑3 (MQ=255)

              tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGa                                                                                      <  1:432570/100‑1 (MQ=255)

              tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTactc                                                                                      <  1:843686/100‑4 (MQ=255)

               ttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaa                                                                                     <  1:258283/100‑1 (MQ=255)

               ttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaa                                                                                     <  2:1036139/100‑1 (MQ=255)

                tAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCACGGTGAGTAATTTTGCCTCGGCAAAATTAGGaaa                                                                                    >  2:2679640/1‑100 (MQ=255)

                 aaTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                   <  1:1145315/100‑1 (MQ=255)

                 aaTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                   >  1:3079993/1‑100 (MQ=255)

                  aTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACa                                                                                  <  1:732997/100‑1 (MQ=255)

                   tCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAg                                                                                 >  1:509418/1‑100 (MQ=255)

                   tCACCACCG‑TTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGa                                                                                <  1:950791/100‑1 (MQ=255)

                    caccacCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGa                                                                                >  1:1055283/1‑100 (MQ=255)

                    caccacCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTAAATCCATTGTCTCTTTCCCTGTGAGTAATTTTAACAAAACAAAATTAGGAAACa                                                                                  <  2:728095/98‑1 (MQ=255)

                          cGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCg                                                                          >  1:369169/1‑100 (MQ=255)

                             ttttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACa                                                                       <  1:1516621/100‑1 (MQ=255)

                              tttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACAt                                                                      >  2:1549484/1‑100 (MQ=255)

                                  cgccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACt                                                                  >  1:1434455/1‑100 (MQ=255)

                                  cgccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACt                                                                  >  1:974223/1‑100 (MQ=255)

                                     cgcacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCaa                                                               >  2:2573602/1‑100 (MQ=255)

                                             aGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGc                                                       <  1:2969008/100‑1 (MQ=255)

                                             aGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGc                                                       <  1:2688887/100‑1 (MQ=255)

                                             aGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGc                                                       <  1:1994246/100‑1 (MQ=255)

                                              gAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCa                                                      >  1:483943/1‑100 (MQ=255)

                                               aaaGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaa                                                     <  1:199904/100‑1 (MQ=255)

                                               aaaGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaa                                                     <  1:3215410/100‑1 (MQ=255)

                                                 aGCAGTCAGCAATTTCATCCAGTGTCGCTTTCCCGGTGATTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAat                                                   <  2:324933/100‑1 (MQ=255)

                                                  gCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAata                                                  <  1:1731633/100‑1 (MQ=255)

                                                       naGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcat                                              >  2:1749393/2‑99 (MQ=255)

                                                       cAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                             <  2:434762/100‑1 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            <  2:2999885/100‑1 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            <  2:2769756/100‑1 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            <  2:1434447/100‑1 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            <  1:540120/100‑1 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            >  1:2080724/1‑100 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            <  2:974216/100‑1 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            >  1:2503764/1‑100 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                            >  1:2678774/1‑100 (MQ=255)

                                                         gCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCCAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCa                                           >  2:496267/1‑100 (MQ=255)

                                                            aTTTCATCCAGTGTCTCTTTCCCGGTGGGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCg                                        >  1:1242415/1‑100 (MQ=255)

                                                                cATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttc                                    <  2:1204104/100‑1 (MQ=255)

                                                                 aTCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttca                                   >  1:2038403/1‑100 (MQ=255)

                                                                    cAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattc                                >  2:2878383/1‑100 (MQ=255)

                                                                       tatctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATACCTCCATGCTTCGCAAATATCATCCAGCGTTTCATTCAAc                             <  2:3033490/98‑1 (MQ=255)

                                                                          ctctTTCCCGTTGAGTAATTTTGCCTCGTCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                          <  2:744278/100‑1 (MQ=255)

                                                                          ctctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                          >  2:1161534/1‑100 (MQ=255)

                                                                          ctctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgc                          <  2:2721707/100‑1 (MQ=255)

                                                                           tctTTCCCGGAGAGTAATTTTGCCTCGGCAAAATTAGGCACCAGATAATAGACATAACTCAATGCTTCGTGAATAACAACCAGCATCTCATTCAAcaggg                         >  2:2040741/1‑96 (MQ=255)

                                                                             ttgcccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGt                       <  2:695844/97‑1 (MQ=255)

                                                                              ttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGtt                      >  2:3059321/1‑100 (MQ=255)

                                                                              ttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGtt                      <  1:3515788/100‑1 (MQ=255)

                                                                               tCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGttt                     >  2:898432/1‑100 (MQ=255)

                                                                                cccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTg                    >  1:2686420/1‑100 (MQ=255)

                                                                                 ccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                   >  2:223070/1‑100 (MQ=255)

                                                                                  cgatctctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                  <  1:2427672/93‑1 (MQ=255)

                                                                                  cGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                  <  1:3309907/100‑1 (MQ=255)

                                                                                  cGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                  <  1:1713052/100‑1 (MQ=255)

                                                                                   ggTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                 >  1:202227/1‑100 (MQ=255)

                                                                                    tttcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCa                >  2:1475291/6‑100 (MQ=39)

                                                                                    tttcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCa                >  1:1083428/6‑100 (MQ=39)

                                                                                      tcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCata              >  1:3509374/4‑100 (MQ=255)

                                                                                      gAGTAATTTTGCCTCGGAAAAATTAGGAAACAGATAAGCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCat               <  2:2649823/99‑1 (MQ=255)

                                                                                        gTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  1:154091/1‑100 (MQ=255)

                                                                                        gTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  2:3563941/1‑100 (MQ=255)

                                                                                        gTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  1:2591078/1‑100 (MQ=255)

                                                                                        ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAtagatcg            >  2:2427656/2‑94 (MQ=255)

                                                                                        ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  1:2884704/2‑100 (MQ=255)

                                                                                        ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  2:642144/2‑100 (MQ=255)

                                                                                        ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  2:591851/2‑100 (MQ=255)

                                                                                        ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc            >  2:2355491/2‑100 (MQ=255)

                                                                                         tAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCg           <  1:3006556/100‑1 (MQ=255)

                                                                                          aaTTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg          >  1:3157524/1‑100 (MQ=255)

                                                                                          aaTTTTGCCGAGGAAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg          >  1:2478164/1‑100 (MQ=255)

                                                                                           aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGc         <  2:647510/100‑1 (MQ=255)

                                                                                           aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGc         <  1:865349/100‑1 (MQ=255)

                                                                                           aTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGc         >  1:172642/1‑100 (MQ=255)

                                                                                            ttttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGca        <  2:3232689/100‑1 (MQ=255)

                                                                                                gCCGCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCCTCCAGCGTTTCATTCAACCGCGGTTTGATCATATCTGca        >  2:3404548/1‑96 (MQ=255)

                                                                                                 ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt   >  2:3581174/1‑100 (MQ=255)

                                                                                                 ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt   >  1:974749/1‑100 (MQ=255)

                                                                                                 ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt   >  2:3231714/1‑100 (MQ=255)

                                                                                                  cTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATa  >  1:760415/1‑100 (MQ=255)

                                                                                                  cTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATa  >  1:3172835/1‑100 (MQ=255)

CCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATA  >  NC_000913/1856234‑1856431

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