BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,949,410	0	C	T	52.2%	96.0 / inf	270	L196L (CTC→CTT)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (70/59); new base T (74/67); total (144/126)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 276 total aligned reads displayed.

GCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCA  >  NC_000913/1949318‑1949501

gCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTCCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTc                                                                                             <  1:3651054/93‑1 (MQ=255)

  cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                                                                           <  2:3122743/93‑1 (MQ=255)

  cTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCaa                                                                                           <  1:2022903/93‑1 (MQ=255)

   tGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAc                                                                                          <  1:3564751/93‑1 (MQ=255)

    gTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcg                                                                                         >  2:871292/1‑93 (MQ=255)

      aGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgcg                                                                                       >  1:435877/1‑93 (MQ=255)

        tGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGtt                                                                                     <  1:2884651/93‑1 (MQ=255)

        tGCCGCCAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGtt                                                                                     <  1:1858026/93‑1 (MQ=255)

         gccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTa                                                                                    <  1:793503/93‑1 (MQ=255)

            gccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGa                                                                                 <  2:2140816/93‑1 (MQ=255)

            gccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGa                                                                                 <  1:4056774/93‑1 (MQ=255)

             ccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAAATCCTCAACGTGTTATGAt                                                                                <  2:4571536/93‑1 (MQ=255)

             ccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAt                                                                                <  2:2276253/93‑1 (MQ=255)

              cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt                                                                               <  1:1790526/93‑1 (MQ=255)

               aGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGAttt                                                                              <  2:1495691/93‑1 (MQ=255)

                gcgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTa                                                                             >  1:375648/1‑93 (MQ=255)

                  gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa                                                                           <  2:2968561/93‑1 (MQ=255)

                  gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                                                                           <  2:3508615/93‑1 (MQ=255)

                  gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                                                                           >  1:4538763/1‑93 (MQ=255)

                  gcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACa                                                                           <  1:503616/93‑1 (MQ=255)

                   ccGAGCAACACAATAACAGCATTAATCATATTTACCCGCGCATCGCCGGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt                                                                          <  2:2838365/93‑1 (MQ=255)

                    cGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATc                                                                         <  1:3116961/93‑1 (MQ=255)

                     gAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCg                                                                        >  1:2156901/1‑93 (MQ=255)

                       gcagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGt                                                                      <  1:1735948/93‑1 (MQ=255)

                       gcagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGt                                                                      <  2:1316973/93‑1 (MQ=255)

                        cagcaCAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTc                                                                     <  2:2660738/93‑1 (MQ=255)

                        caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                                                                     <  1:1704208/90‑1 (MQ=255)

                          acacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCtg                                                                   >  2:1550099/2‑92 (MQ=255)

                           cacaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAc                                                                  <  2:4544535/93‑1 (MQ=255)

                            acaATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTAcc                                                                 >  2:3425573/1‑93 (MQ=255)

                             caATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCg                                                                >  1:1404891/1‑93 (MQ=255)

                              aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTATCGc                                                               >  1:2063798/1‑93 (MQ=255)

                               aTAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCa                                                              >  1:2313919/1‑93 (MQ=255)

                                tAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCaa                                                             >  1:4581117/1‑93 (MQ=255)

                                tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGTAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                             <  1:3295921/93‑1 (MQ=255)

                                tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                             >  2:3080838/1‑93 (MQ=255)

                                 aaCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAAt                                                            <  2:2504081/93‑1 (MQ=255)

                                  aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATg                                                           >  2:1928487/1‑93 (MQ=255)

                                     gTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTc                                                        >  2:485018/1‑93 (MQ=255)

                                      tATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCg                                                       <  1:342489/93‑1 (MQ=255)

                                      tATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGCTATGATTTACATCGGTCTACCGCAATGGTCg                                                       >  2:3169880/1‑93 (MQ=255)

                                         taatcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCAt                                                    >  1:4375826/5‑93 (MQ=255)

                                          aaCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTATACCGCAATGGTCGCATc                                                   >  2:1444889/1‑93 (MQ=255)

                                           atcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATcc                                                  <  1:1300093/91‑1 (MQ=255)

                                           aCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATcc                                                  <  1:1610068/93‑1 (MQ=255)

                                            tcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCt                                                 <  2:2093367/92‑1 (MQ=255)

                                             cATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTa                                                >  2:409379/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaaa                                              >  2:1358332/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa                                              >  1:4304441/1‑93 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa                                              <  2:2919525/93‑1 (MQ=255)

                                               tatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTaaa                                              <  2:3706590/93‑1 (MQ=255)

                                                 tCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATg                                            >  1:2757067/1‑93 (MQ=255)

                                                  ctccccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATgg                                           <  2:975968/90‑1 (MQ=255)

                                                  cTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATgg                                           >  2:2875700/1‑93 (MQ=255)

                                                   tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                          >  1:2209913/1‑93 (MQ=255)

                                                   tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATggg                                          <  2:631424/93‑1 (MQ=255)

                                                    tcccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAAGGGGt                                         <  2:910239/92‑1 (MQ=255)

                                                     cccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTg                                        >  2:4385550/1‑93 (MQ=255)

                                                        gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCgg                                     >  2:1034765/1‑93 (MQ=255)

                                                         cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGt                                    >  2:1758732/1‑93 (MQ=255)

                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                   <  2:1116407/93‑1 (MQ=255)

                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                   <  2:4263616/93‑1 (MQ=255)

                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                   >  2:2243134/1‑93 (MQ=255)

                                                          gcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGtt                                   >  2:3656793/1‑93 (MQ=255)

                                                           cATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTg                                  >  1:1171915/1‑93 (MQ=255)

                                                            aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                 <  1:4477174/93‑1 (MQ=255)

                                                            aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                 >  2:2577648/1‑93 (MQ=255)

                                                            aTCGCCCGTCTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                 <  2:1604291/93‑1 (MQ=255)

                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                >  2:3098493/1‑93 (MQ=255)

                                                             tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                <  1:4520386/93‑1 (MQ=255)

                                                               gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggta                              <  1:2494259/93‑3 (MQ=255)

                                                                cccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTAt                             >  2:1649340/1‑93 (MQ=255)

                                                                 ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                            >  1:2674326/1‑93 (MQ=255)

                                                                    tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                         >  1:2995539/1‑93 (MQ=255)

                                                                    tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCAcc                         >  2:4355435/1‑93 (MQ=255)

                                                                    tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCAcc                         >  1:2275963/1‑93 (MQ=255)

                                                                        cGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGcc                     <  1:3833181/93‑1 (MQ=255)

                                                                         gTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGcgt                    <  2:1379360/93‑3 (MQ=255)

                                                                          gaGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCtt                   <  2:2139327/92‑1 (MQ=255)

                                                                             cGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGaa                <  2:1826249/93‑1 (MQ=255)

                                                                             cGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGaa                <  2:3769942/93‑1 (MQ=255)

                                                                              gTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAag               <  2:56434/93‑1 (MQ=255)

                                                                              gTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAag               >  2:590877/1‑93 (MQ=255)

                                                                              gTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAag               >  1:1389682/1‑93 (MQ=255)

                                                                               tagAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAaga              >  2:2446337/3‑93 (MQ=255)

                                                                               tGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAaga              <  1:883556/93‑1 (MQ=255)

                                                                                ggAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagag             <  1:1979535/93‑1 (MQ=255)

                                                                                 gAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGt            >  1:3500956/1‑93 (MQ=255)

                                                                                 gAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGt            >  1:3862481/1‑93 (MQ=255)

                                                                                 gAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGt            >  1:3224698/1‑93 (MQ=255)

                                                                                 gAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGt            <  2:3232730/93‑1 (MQ=255)

                                                                                 gAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGt            >  2:1999462/1‑93 (MQ=255)

                                                                                   agagATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTAt          >  2:3928069/1‑93 (MQ=255)

                                                                                   agagATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTAt          >  1:1032276/1‑93 (MQ=255)

                                                                                      gaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGcc       >  2:3527746/1‑93 (MQ=255)

                                                                                      gaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGcc       >  2:3309344/1‑93 (MQ=255)

                                                                                       aTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGccc      >  1:3470206/1‑93 (MQ=255)

                                                                                        tCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCg     >  1:1026161/1‑93 (MQ=255)

                                                                                         ccTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGc    >  1:3878932/1‑93 (MQ=255)

                                                                                           tCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGccc  >  1:917902/1‑92 (MQ=255)

GCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCA  >  NC_000913/1949318‑1949501

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