BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,856,333	T→G	15.9%	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	15.9%	434.8 / 72.2	182	E199A (GAG→GCG)	ydjH	putative kinase
Reads supporting (aligned to +/- strand): ref base T (77/74); new base G (10/19); total (89/93)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 186 total aligned reads displayed.

CCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAA  >  NC_000913/1856234‑1856432

ccGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCt                                                                                                     <  1:2345411/100‑1 (MQ=255)

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

   tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTc                                                                                                    >  1:1572655/1‑98 (MQ=255)

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

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

    cTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCggn                                                                                                 <  2:315353/100‑2 (MQ=255)

       tACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaa                                                                                              >  2:2740066/1‑100 (MQ=255)

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

         cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                            <  2:1772980/100‑1 (MQ=255)

         cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                            <  2:3322121/100‑1 (MQ=255)

         cacGTTTTAATCACCACGGTGTTTACGCCGCACCCAAGAAACCAGTAAGCAATTTCCTCCAGTGTCTCTTTCCCGGCGATTAATTTGACCTCGGc                                                                                                 >  1:1071670/3‑95 (MQ=255)

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

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

            gTTTTCATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                         <  2:4065885/100‑2 (MQ=255)

            gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAATCAGTCAGCAATTTCATCCATTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                         <  1:2332511/100‑2 (MQ=255)

            gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAg                                                                                         <  2:1974514/100‑1 (MQ=255)

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

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

              tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGa                                                                                       >  2:3738900/1‑100 (MQ=255)

              tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTactc                                                                                       <  2:1486659/100‑4 (MQ=255)

               ttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAgg                                                                                        <  2:586766/98‑1 (MQ=255)

                 aaTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTactcacc                                                                                    <  1:173456/100‑7 (MQ=38)

                   tCACCACCGTTTGCGCGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                    <  2:3884539/98‑1 (MQ=255)

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

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

                     accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGAt                                                                                <  1:136443/100‑1 (MQ=255)

                      ccaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATa                                                                               >  1:2727115/1‑100 (MQ=255)

                        accGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAACATTAGGAAACAGATaa                                                                              >  2:774788/1‑99 (MQ=255)

                          cGTTTTTACGCCGCACGCAAGAAAGCAGGCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGAt                                                                                <  1:3868319/95‑1 (MQ=255)

                           gTTTTTACTCCGCACGCAAGAAAGCAGTTAGCAATTTCATCCAGTTTCTCTTTCCCAGTGAGTAATTTTGCCTCCGCAAAATTAGGAAACAGATACTCg                                                                           <  1:794843/99‑1 (MQ=255)

                           gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                          <  1:1023093/100‑1 (MQ=255)

                           gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGAAAAATTAGGAAACAGATAATCGa                                                                          <  1:628512/100‑1 (MQ=255)

                            tttttACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGAc                                                                         <  2:1269784/100‑1 (MQ=255)

                                   gccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                  <  1:2581404/100‑1 (MQ=255)

                                       cacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                               >  1:3791520/1‑99 (MQ=255)

                                        acgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTTCCTCGGCAAAATTACGAAACAGATAATCGACATAACTCAATGc                                                             >  1:410651/1‑100 (MQ=255)

                                         cgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCGAAATTAGGAAACAGATAATCGACATAACTCAATGCt                                                            >  1:4200401/1‑100 (MQ=255)

                                         cgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCt                                                            >  1:1015951/1‑100 (MQ=255)

                                           caAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                          >  2:2093722/1‑100 (MQ=255)

                                           caAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTc                                                          <  1:2185896/100‑1 (MQ=255)

                                            aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCg                                                         <  2:2163035/100‑1 (MQ=255)

                                             aGAAAGCAGTCAGCAATATCAACCAGTGTCCCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAAGCGACATAACTCAATGCTTCGc                                                        <  1:741181/100‑1 (MQ=255)

                                                 cgCAGTCAGCAATTTCATCAAGTGTCTCTTTCCCGGAGAGTAATTTTTCCTCTACAAAATTAGGAAACAGCTCATCGACATAACTCAGTGCTTCGCaact                                                    >  1:1177146/2‑98 (MQ=255)

                                                   cAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAatat                                                  <  1:1200766/100‑1 (MQ=255)

                                                     gTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatca                                                >  2:1029109/1‑100 (MQ=255)

                                                     gTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatca                                                >  2:2627378/1‑100 (MQ=255)

                                                     gTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatca                                                >  1:502026/1‑100 (MQ=255)

                                                       cAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                              >  2:3464767/1‑100 (MQ=255)

                                                       cAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcat                                               >  2:2790580/1‑99 (MQ=255)

                                                        aGCAATTTCATCCAGTGTCTCTTTCCCAGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                             <  1:498117/100‑1 (MQ=255)

                                                          cAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAg                                           >  2:3727503/1‑100 (MQ=255)

                                                               tCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGttt                                      <  2:50814/100‑1 (MQ=255)

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

                                                                cATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttc                                     >  2:436221/1‑100 (MQ=255)

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

                                                                cATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCACGGCAAAATTAGGAAACAGATAGTCGACATAACTCAAGGCTTCGCGAATATCATCCAGCGta                                       >  1:4170959/1‑97 (MQ=255)

                                                                  tCCTGTGTCTCCTTACCGGTGATTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcat                                   <  1:3837979/100‑1 (MQ=255)

                                                                  tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcat                                   <  1:2390448/100‑1 (MQ=255)

                                                                  tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcat                                   <  2:2452877/100‑1 (MQ=255)

                                                                  tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGGttca                                    >  1:996967/1‑99 (MQ=255)

                                                                    cAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattc                                 <  2:3943905/100‑1 (MQ=255)

                                                                     aGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattca                                <  2:707186/100‑1 (MQ=255)

                                                                     aGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattca                                <  2:2498042/100‑1 (MQ=255)

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

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

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

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

                                                                            ctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCgg                         <  1:2958044/100‑1 (MQ=255)

                                                                                cccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTg                     <  1:1481408/100‑1 (MQ=255)

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

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

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

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

                                                                                  tgtttcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGTAATCACTCAATGCTTCGCAAACATCATCCAGCGTTTCATTCAACTGCGGTTTGa                    >  1:3689081/8‑99 (MQ=25)

                                                                                  cGGTGAGTGATGTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTGATTCAACCGCGGTTTGAt                   >  1:4029100/1‑100 (MQ=255)

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

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

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

                                                                                   agTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                  <  2:2640820/99‑1 (MQ=255)

                                                                                     ttcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCat                >  2:619079/5‑100 (MQ=39)

                                                                                     tGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCat                >  2:3665588/1‑100 (MQ=255)

                                                                                      gAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCata               <  2:2830268/100‑1 (MQ=255)

                                                                                       aGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCatat              >  2:901508/1‑100 (MQ=255)

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

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

                                                                                        ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCatat              >  1:3822833/2‑99 (MQ=255)

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

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

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

                                                                                          aaTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg           <  1:2098052/100‑1 (MQ=255)

                                                                                          aaTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg           >  1:1337287/1‑100 (MQ=255)

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

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

                                                                                             tttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGcac        <  2:4240704/100‑1 (MQ=255)

                                                                                             tttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGca         >  1:2451978/1‑99 (MQ=255)

                                                                                               tGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAg      <  2:1206545/100‑1 (MQ=255)

                                                                                                gCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGa     >  2:2760087/1‑100 (MQ=255)

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

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

                                                                                                   tCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATaa  >  1:1746106/1‑100 (MQ=255)

CCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAA  >  NC_000913/1856234‑1856432

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: `‑`