BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,542,521	A→G	15.5%	H162R (CAT→CGT)	fimE →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,521	0	A	G	15.5%	458.1 / 65.0	155	H162R (CAT→CGT)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand): ref base A (70/61); new base G (13/11); total (83/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 156 total aligned reads displayed.

AACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAA  >  NC_000913/4542426‑4542619

 cCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGa                                                                                                <  1:249669‑M2/98‑1 (MQ=255)

 aCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGa                                                                                                <  2:1660441/99‑1 (MQ=255)

 aCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGa                                                                                                >  2:1999317/1‑99 (MQ=255)

 aCGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGa                                                                                                >  1:1821677/1‑99 (MQ=255)

  cGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGaa                                                                                               <  2:2097200/99‑1 (MQ=255)

   gCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGaaa                                                                                              <  1:602714/99‑1 (MQ=255)

     aGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAata                                                                                            >  1:2718119/1‑99 (MQ=255)

      gACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAatat                                                                                           >  2:1185169/1‑99 (MQ=255)

       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATAtt                                                                                          >  1:1665402/1‑99 (MQ=255)

       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAtt                                                                                          <  1:773240/99‑1 (MQ=255)

       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAtt                                                                                          <  2:1067874/99‑1 (MQ=255)

       aCTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATAtt                                                                                          >  2:1601409/1‑99 (MQ=255)

        cTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTc                                                                                         >  1:2341888/1‑99 (MQ=255)

         tCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCg                                                                                        >  2:3193060/1‑99 (MQ=255)

         tCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCg                                                                                        >  2:1486718/1‑99 (MQ=255)

           aTCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGcc                                                                                      >  2:2528440/1‑99 (MQ=255)

           aTCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGcc                                                                                      <  1:386835/99‑1 (MQ=255)

              cTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATa                                                                                   <  1:2859234/99‑1 (MQ=255)

              cTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATa                                                                                   <  1:777851/99‑1 (MQ=255)

               tCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATAc                                                                                  <  1:2113375/99‑1 (MQ=255)

               tCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATAc                                                                                  <  2:1840071/99‑1 (MQ=255)

                cATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACt                                                                                 <  1:3033723/99‑1 (MQ=255)

                    tGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGc                                                                             >  1:2218790/1‑99 (MQ=255)

                    tGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGc                                                                             >  2:2559445/1‑99 (MQ=255)

                    tGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGc                                                                             >  2:2895288/1‑99 (MQ=255)

                      ttAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGt                                                                           <  1:2983314/99‑1 (MQ=255)

                        aaGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTta                                                                         <  2:904503/99‑1 (MQ=255)

                        aaGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTta                                                                         <  2:2299513/99‑1 (MQ=255)

                        aaGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTta                                                                         <  2:2656012/99‑1 (MQ=255)

                        aaGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTta                                                                         >  2:489691/1‑99 (MQ=255)

                         aGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTtat                                                                        >  2:2822461/1‑99 (MQ=255)

                         aGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTtat                                                                        >  2:1285258/1‑99 (MQ=255)

                          ggCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTtata                                                                       <  1:3032232/99‑1 (MQ=255)

                           gCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATAc                                                                      <  1:1285597/99‑1 (MQ=255)

                           gCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATAc                                                                      <  1:545218/99‑1 (MQ=255)

                            cATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATAcc                                                                     <  1:3343895/99‑1 (MQ=255)

                             aTGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCg                                                                    >  2:124661/1‑99 (MQ=255)

                              tGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGc                                                                   >  1:816789/1‑99 (MQ=255)

                              tGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGc                                                                   >  2:2043721/1‑99 (MQ=255)

                                cTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCa                                                                 >  1:1425876/1‑99 (MQ=255)

                                 ttGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAg                                                                <  1:3397158/99‑1 (MQ=255)

                                      gTTCTGAATTGGCGGAGCGTGGAACAGATAATCGCTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAta                                                           >  2:3416939/1‑98 (MQ=255)

                                      gTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtg                                                           <  2:1164991/99‑1 (MQ=255)

                                      gTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtg                                                           <  2:352676/99‑1 (MQ=255)

                                       ttATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgc                                                          >  1:1935291/1‑99 (MQ=255)

                                       ttATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgc                                                          >  2:1621503/1‑99 (MQ=255)

                                          tGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAAtgctgc                                                       >  2:1089041/1‑99 (MQ=255)

                                             aTTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCg                                                    >  1:1022901/1‑99 (MQ=255)

                                             aTTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCg                                                    <  1:1632483/99‑1 (MQ=255)

                                               tGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGtt                                                  <  2:1676749/99‑1 (MQ=255)

                                                ggCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGttt                                                 <  2:682021/99‑1 (MQ=255)

                                                ggCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGttt                                                 <  1:3094923/99‑1 (MQ=255)

                                                   ggAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGc                                              >  1:1818956/1‑99 (MQ=255)

                                                   ggAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGc                                              <  1:2397748/99‑1 (MQ=255)

                                                      gCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCgg                                           >  2:1845912/1‑99 (MQ=255)

                                                      gCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCgg                                           >  2:2700041/1‑99 (MQ=255)

                                                       cGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGa                                          <  1:2980795/99‑1 (MQ=255)

                                                        gtggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGAt                                         <  2:1021906/99‑1 (MQ=255)

                                                        gtggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGAt                                         >  2:3177683/1‑99 (MQ=255)

                                                         tggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGAtt                                        >  2:283850/1‑99 (MQ=255)

                                                         tggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGAtt                                        <  2:2691735/99‑1 (MQ=255)

                                                          ggtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTa                                       <  1:2250837/99‑1 (MQ=255)

                                                           gtgCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTAt                                      >  1:144475/1‑99 (MQ=255)

                                                            tgCAGATACTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATg                                     >  2:41979/1‑99 (MQ=255)

                                                              cAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATggg                                   >  2:2349184/1‑99 (MQ=255)

                                                               aGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGga                                  <  1:283005/99‑1 (MQ=255)

                                                                gATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaa                                 >  1:1371459/1‑99 (MQ=255)

                                                                gATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaa                                 <  2:78831/99‑1 (MQ=255)

                                                                gATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaa                                 >  2:2182492/1‑99 (MQ=255)

                                                                  tACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaag                               <  2:3120886/99‑1 (MQ=255)

                                                                   aCTCGTTTAATTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaaga                              <  2:54350/99‑1 (MQ=255)

                                                                   aCTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaaga                              <  2:984910/99‑1 (MQ=255)

                                                                   aCTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaaga                              <  1:547877/99‑1 (MQ=255)

                                                                   aCTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaaga                              <  1:3469357/99‑1 (MQ=255)

                                                                    cTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaagaa                             <  2:3418329/99‑1 (MQ=255)

                                                                    cTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGgaaagaa                             <  1:1587697/99‑1 (MQ=255)

                                                                       gTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAaata                          >  1:2064676/1‑99 (MQ=255)

                                                                        tttAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAaataa                         >  2:75716/1‑99 (MQ=255)

                                                                        tttAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAaataa                         <  1:627574/99‑1 (MQ=255)

                                                                        tttAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAaataa                         >  2:174220/1‑99 (MQ=255)

                                                                           aaTTCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAAtct                      >  1:1510098/1‑99 (MQ=255)

                                                                           aaTTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAAtct                      <  2:1521860/99‑1 (MQ=255)

                                                                            aTTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAAtctc                     >  2:2575726/1‑99 (MQ=255)

                                                                             ttCAGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCa                    <  1:3212410/99‑1 (MQ=255)

                                                                             ttCAGGATTATCTCGGGCATCGAAATATTCTCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCa                    <  2:140052/99‑1 (MQ=255)

                                                                             ttCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCa                    <  2:164860/99‑1 (MQ=255)

                                                                              tCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCAt                   <  1:937553/99‑1 (MQ=255)

                                                                                aGGATTATCTCGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATaa                 >  2:2445860/1‑99 (MQ=255)

                                                                                 ggATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATaaa                <  1:2257905/99‑1 (MQ=255)

                                                                                 ggATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATaaa                >  1:668175/1‑99 (MQ=255)

                                                                                  gATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAAc               >  1:101586/1‑99 (MQ=255)

                                                                                   aTTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACg              >  1:1971825/1‑99 (MQ=255)

                                                                                       tctcGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGaaaa          >  1:1789752/1‑99 (MQ=255)

                                                                                       tctcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGaaaa          >  2:2115611/1‑99 (MQ=255)

                                                                                       tctcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGAGAAAGAAATAATCTCATAAACGaaaa          >  1:1791957/1‑99 (MQ=255)

                                                                                         tcGGGCGTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAAt        <  2:3072767/99‑1 (MQ=255)

                                                                                         tcGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAAt        <  1:2886033/99‑1 (MQ=255)

                                                                                           gggCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTa      >  2:904377/1‑99 (MQ=255)

                                                                                            ggCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTaa     <  1:328172/99‑1 (MQ=255)

                                                                                               aTCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTaaaaa  >  1:1480601/1‑99 (MQ=255)

AACGCAGACTCATCCTCATATGTTAAGGCATGCTTGCGGTTATGAATTGGCGGAGCGTGGTGCAGATACTCGTTTAATTCAGGATTATCTCGGGCATCGAAATATTCGCCATACTGTGCGTTATACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATGGGAAAGAAATAATCTCATAAACGAAAAATTAAAAA  >  NC_000913/4542426‑4542619

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