BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,071,190	T→C	T65A (ACG→GCG)	b1010 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,071,190	0	T	C	99.2%	385.5 / ‑4.6	125	T65A (ACG→GCG)	b1010	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (67/57); total (67/58)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

Only 100 of 134 total aligned reads displayed.

CGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGATCAAAA  >  NC_000913/1071100‑1071278

cGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGcct                                                                                         <  1:764024/93‑1 (MQ=255)

   cGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGa                                                                                      <  1:3381090/93‑1 (MQ=255)

     tAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATc                                                                                    <  1:4135724/93‑1 (MQ=255)

      aGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCa                                                                                   >  2:381471/1‑93 (MQ=255)

        aTTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCAcc                                                                                 >  2:38937/1‑93 (MQ=255)

         tttCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGACGCCTCGATCACCt                                                                                <  1:1041145/93‑1 (MQ=255)

         ttgcGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCt                                                                                <  1:42491/90‑1 (MQ=255)

            cGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGTAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGc                                                                             <  1:3042209/93‑1 (MQ=255)

            cGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGc                                                                             <  2:2641466/93‑1 (MQ=255)

             gTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCg                                                                            >  1:3852793/1‑93 (MQ=255)

             gTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCg                                                                            >  2:875038/1‑93 (MQ=255)

              ttAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCgg                                                                           >  2:1194117/1‑93 (MQ=255)

              ttAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCgg                                                                           >  2:3775401/1‑93 (MQ=255)

               tAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGa                                                                          <  1:2732449/93‑1 (MQ=255)

                aaTCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGAt                                                                         >  1:2174696/1‑93 (MQ=255)

                 aTCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATa                                                                        <  2:4239605/93‑1 (MQ=255)

                 aTCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATa                                                                        <  2:2318646/93‑1 (MQ=255)

                 aTCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATa                                                                        >  2:4618757/1‑93 (MQ=255)

                  tCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAg                                                                       <  2:2989091/93‑1 (MQ=255)

                  tCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAg                                                                       >  2:4235010/1‑93 (MQ=255)

                    gcggcgTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGtt                                                                     >  2:4016638/1‑93 (MQ=255)

                     cggcgTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGttt                                                                    >  1:1043622/1‑93 (MQ=255)

                       gcgTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTcc                                                                  >  2:4015844/1‑93 (MQ=255)

                        cgTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCa                                                                 <  1:1017540/93‑1 (MQ=255)

                         gTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAg                                                                >  1:3468881/1‑93 (MQ=255)

                          tAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGa                                                               >  2:3982362/1‑93 (MQ=255)

                          tAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCAGCCGCCTCGATCACCTTGCGGATAGTTTCAAGa                                                               >  1:901346/1‑93 (MQ=255)

                           aaTTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGaa                                                              >  1:4505342/1‑93 (MQ=255)

                            aTTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAAc                                                             <  1:3247772/93‑1 (MQ=255)

                             tttttCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACg                                                            >  1:294537/1‑93 (MQ=255)

                             tttttCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACg                                                            >  1:4628478/1‑93 (MQ=255)

                              ttttCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGt                                                           <  1:3828279/93‑1 (MQ=255)

                              ttttCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGt                                                           <  1:60359/93‑1 (MQ=255)

                                   cAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCgg                                                      >  1:1136001/1‑93 (MQ=255)

                                   cAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCgg                                                      <  1:1799807/93‑1 (MQ=255)

                                     gTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGt                                                    >  1:2921438/1‑93 (MQ=255)

                                     gTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGt                                                    >  1:1463977/1‑93 (MQ=255)

                                      tCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCggggt                                                   >  1:89014/1‑91 (MQ=255)

                                        ggTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTg                                                 <  1:22211/93‑1 (MQ=255)

                                        ggTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTg                                                 >  1:2724736/1‑93 (MQ=255)

                                         gTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTgg                                                >  1:828949/1‑92 (MQ=255)

                                          taataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCGCCCGCCGCCTTGATCACCTTGCGGATAGTTTCCAGAACGTGGCCGGTTTgcg                                               <  2:4548882/93‑1 (MQ=255)

                                          taataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTgcg                                               <  1:4243813/93‑1 (MQ=255)

                                          taataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTgcg                                               <  2:698271/93‑1 (MQ=255)

                                          taataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTgcg                                               >  2:3743420/1‑93 (MQ=255)

                                           aataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTgcgc                                              >  1:1356478/1‑93 (MQ=255)

                                            ataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTggggc                                             >  1:3592056/1‑89 (MQ=255)

                                            ataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTgcctc                                             >  1:1574/1‑90 (MQ=255)

                                            ataaAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGcc                                             >  2:610693/1‑93 (MQ=255)

                                              aaaGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCtt                                           <  2:3321709/93‑1 (MQ=255)

                                              aaaGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCtt                                           >  2:2357536/1‑93 (MQ=255)

                                              aaaGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCtt                                           >  2:3582320/1‑93 (MQ=255)

                                               aaGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCttt                                          <  2:4175338/93‑1 (MQ=255)

                                                 gATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGGCttttg                                        >  2:1661129/1‑91 (MQ=255)

                                                 gATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTgg                                        <  2:4665349/93‑1 (MQ=255)

                                                   tGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGGCTTTgggg                                      >  1:1113188/1‑92 (MQ=255)

                                                   tGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGt                                      <  2:2192179/93‑1 (MQ=255)

                                                    gCTGTTGAAGGCCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGTGGCCTTGGGTc                                     >  1:1679057/1‑93 (MQ=255)

                                                     cTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCa                                    >  1:1972634/1‑93 (MQ=255)

                                                      tGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTTGGTCAt                                   >  1:4169771/1‑93 (MQ=255)

                                                       gTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑c                                 <  1:376462/93‑1 (MQ=255)

                                                        ttGAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑Cg                                <  2:2817341/93‑1 (MQ=255)

                                                          gAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGc                              >  1:2882328/1‑93 (MQ=255)

                                                          gAAGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTGGGTCAT‑CGGc                              >  2:2977346/1‑93 (MQ=255)

                                                            aGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCaa                            <  2:3499937/93‑1 (MQ=255)

                                                            aGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCaa                            >  1:2279683/1‑93 (MQ=255)

                                                            aGGTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTTGGTCAT‑CGGCaa                            >  2:3084342/1‑93 (MQ=255)

                                                             ggTCACATCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTTGGCCTTTTGGTCAT‑CGGCaca                           >  1:3986984/1‑91 (MQ=255)

                                                                cacaTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGTGCCTTTTGGTCATCCGGCaaaag                         >  1:4652868/1‑91 (MQ=255)

                                                                cacaTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTGGGTCAT‑CGGCAAACAg                        >  2:3548766/1‑93 (MQ=255)

                                                                cacaTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTGGGGCCT‑CGGCAAAccg                        >  1:3216072/1‑91 (MQ=255)

                                                                 acaTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGc                       <  2:554956/93‑1 (MQ=255)

                                                                 acaTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGACCTTTGGTCAC‑AGGCaaagaac                       >  2:1845162/1‑89 (MQ=255)

                                                                  caTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTGGGGCAT‑CGGGAAACAGCa                      >  2:1999070/1‑93 (MQ=255)

                                                                   aTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCTCCTTTGGGACAT‑CGGCAAACAGCAc                     >  2:1453839/1‑93 (MQ=255)

                                                                   aTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCAc                     <  1:1356564/93‑1 (MQ=255)

                                                                   aTCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGCTTTGGCCTTTGGGTCAT‑CGGCAAACAGCAc                     >  1:2037191/1‑93 (MQ=255)

                                                                    tCCGCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTTCGCCTTTGGGTCAT‑CGGCAAACAGCACg                    >  2:1050909/1‑93 (MQ=255)

                                                                      cGCCATCGTGCCACCCGCCGCCGCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGtt                  <  2:3287090/93‑1 (MQ=255)

                                                                       gCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGtta                 <  2:3613077/93‑1 (MQ=255)

                                                                       gCCATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGtta                 <  1:4440189/93‑1 (MQ=255)

                                                                        ccATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGttat                <  2:110571/93‑1 (MQ=255)

                                                                         cATCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGttatt               <  2:2484771/93‑1 (MQ=255)

                                                                           tCGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTtattat             <  1:2728401/93‑1 (MQ=255)

                                                                            cGTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATg            <  2:1621780/93‑1 (MQ=255)

                                                                             gTGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCTTTGGGGCAT‑CGGCAAACAGCCCGTTATTATGt           >  1:748577/1‑93 (MQ=255)

                                                                              tGCCCCCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGtt          <  2:3506461/93‑1 (MQ=255)

                                                                              tGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGtt          <  1:1125755/93‑1 (MQ=255)

                                                                              tGCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGtt          <  2:976296/93‑1 (MQ=255)

                                                                               gCCACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTg         >  1:2001635/1‑93 (MQ=255)

                                                                                ccACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGa        <  2:2898283/93‑1 (MQ=255)

                                                                                ccACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGa        <  1:38937/93‑1 (MQ=255)

                                                                                 cACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGAt       <  1:4622969/93‑1 (MQ=255)

                                                                                 cACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGAt       <  1:2938883/93‑1 (MQ=255)

                                                                                 cACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGAt       >  1:1521086/1‑93 (MQ=255)

                                                                                 cACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGGCTTTGGGGCAT‑CGGCAAACAGCACGTTATTATGTTTAt       >  1:852407/1‑93 (MQ=255)

                                                                                 cACCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGGTTGCGCCCTTGGGGCAT‑CGGCAAACAGCACGTTATTATGTTTAt       >  1:1675895/1‑93 (MQ=255)

                                                                                  aCCCGCCGCCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGAAAACAGCACGTTATTATGTTGATc      <  1:1376353/93‑1 (MQ=255)

                                                                                      gccgccTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGATCaaaa  <  2:3417130/93‑1 (MQ=255)

                                                                                      gccgccTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGATCaaaa  <  1:352802/93‑1 (MQ=255)

CGGCGTAGATTTCGTTAATCGCGGCGTAATTTTTCCAGTCGGTAATAAAGATGCTGTTGAAGGTCACATCCGCCATCGTGCCACCCGCCGTCTCGATCACCTTGCGGATAGTTTCCAGAACGTGGCGGGTTTGCGCCTTTGGGTCAT‑CGGCAAACAGCACGTTATTATGTTGATCAAAA  >  NC_000913/1071100‑1071278

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