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	66.1%	602.0 / inf	451	L196L (CTC→CTT)	yecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (90/63); new base T (183/115); total (273/178)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 461 total aligned reads displayed.

TGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCA  >  NC_000913/1949321‑1949501

  tAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAcgc                                                                                        <  1:3923616/93‑1 (MQ=255)

       ccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                                                                   <  2:1867003/93‑1 (MQ=255)

       ccgccAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTAt                                                                                   <  2:816998/93‑1 (MQ=255)

       ccgccAGCGCCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTAt                                                                                   >  2:4654638/1‑93 (MQ=255)

           cAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGAtt                                                                               >  1:843820/1‑93 (MQ=255)

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

             gcgcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTa                                                                             <  2:4355313/93‑1 (MQ=255)

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

              cgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCGCGCGTTATGATTTAc                                                                            <  2:1592629/93‑1 (MQ=255)

              cgcCGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTAc                                                                            >  1:2008401/1‑93 (MQ=255)

               gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa                                                                           >  2:1456629/1‑93 (MQ=255)

               gcCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACa                                                                           >  1:195139/1‑93 (MQ=255)

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

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

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

                ccGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACAt                                                                          <  2:3375867/93‑1 (MQ=255)

                ccGAGCAACACAATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAt                                                                          >  1:1650041/1‑93 (MQ=255)

                 cGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATc                                                                         <  2:441023/93‑1 (MQ=255)

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

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

                   agcagcACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCgg                                                                       <  1:460286/93‑1 (MQ=255)

                    gcagcaCAATAACAGTATCAACCATATCTACCCGCGCATCCCCCGTATGCGGCGCGCGGACGAGCTCCTCAACGCGTTATGATTTACAtttgt                                                                      >  2:4734241/1‑89 (MQ=255)

                     caacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTc                                                                     >  1:3512033/4‑93 (MQ=255)

                      aacacaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCt                                                                    >  1:1872695/3‑93 (MQ=255)

                         acaATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGcc                                                                 <  1:2637264/93‑1 (MQ=255)

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

                          caATAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCg                                                                >  2:2236536/1‑93 (MQ=255)

                           aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc                                                               >  1:2516727/1‑93 (MQ=255)

                           aaTAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGc                                                               >  1:1007355/1‑93 (MQ=255)

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

                             tAACAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCaa                                                             <  1:432443/93‑1 (MQ=255)

                               aCAGCATTAATCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTAGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATg                                                           <  1:3709821/93‑1 (MQ=255)

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

                                     tCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCa                                                     >  1:2441074/1‑93 (MQ=255)

                                       aaCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAAGGCATTATGATTTACATCGGTCTACCGCAATGGTCGCATc                                                   >  2:1321114/1‑93 (MQ=255)

                                       aaCCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATc                                                   >  1:1288858/1‑93 (MQ=255)

                                        atcATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATcc                                                  >  1:1409267/3‑93 (MQ=255)

                                         ccATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCt                                                 >  1:3385287/1‑93 (MQ=255)

                                           atatCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTaa                                               >  2:3162149/1‑93 (MQ=255)

                                             atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAAt                                             >  2:2899632/1‑93 (MQ=255)

                                             atCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAAt                                             >  1:775359/1‑93 (MQ=255)

                                              tCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCCTCCTAAATg                                            >  2:883006/1‑93 (MQ=255)

                                               cTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATgg                                           >  1:4533199/1‑93 (MQ=255)

                                                tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATggg                                          <  1:1526321/93‑1 (MQ=255)

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

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

                                                tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATggg                                          >  1:2697485/1‑93 (MQ=255)

                                                tACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATggg                                          >  2:15141/1‑93 (MQ=255)

                                                 aCCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGt                                         <  2:3728180/93‑1 (MQ=255)

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

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

                                                   ccGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGc                                       <  1:4466166/93‑1 (MQ=255)

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

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

                                                     gcgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCgg                                     >  1:513725/1‑93 (MQ=255)

                                                      cgcATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGt                                    <  2:1937024/93‑1 (MQ=255)

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

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

                                                         aTCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgg                                 <  2:1739653/93‑1 (MQ=255)

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

                                                          tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                <  2:4250826/93‑1 (MQ=255)

                                                          tCGCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTggg                                <  2:3225245/93‑1 (MQ=255)

                                                            gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggta                              >  2:1365579/1‑91 (MQ=255)

                                                            gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggta                              >  1:4411444/1‑91 (MQ=255)

                                                            gCCCGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTgggta                              >  1:3036112/1‑91 (MQ=255)

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

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

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

                                                              ccGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                            <  2:114906/93‑1 (MQ=255)

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

                                                              cagtgtGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACAGCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATc                            <  1:2223784/91‑1 (MQ=255)

                                                               cGTGTGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCa                           >  1:4261296/1‑93 (MQ=255)

                                                                gtgtGCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATTGGTGCGGTTGGGTATCAc                          >  2:2571965/1‑93 (MQ=255)

                                                                 tgtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCAcc                         <  1:2819714/93‑1 (MQ=255)

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

                                                                 tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATTAcc                         <  2:4138742/93‑1 (MQ=255)

                                                                 tgtgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTAGCATCCTAAATGGGTGCGGTTGGGGATCAcc                         >  2:2629308/1‑93 (MQ=255)

                                                                  gtgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCa                        >  1:2156294/1‑93 (MQ=255)

                                                                   tgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAg                       >  2:564204/1‑93 (MQ=255)

                                                                   tgCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAg                       >  2:4499170/1‑93 (MQ=255)

                                                                   tgCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAg                       <  1:1807415/93‑1 (MQ=255)

                                                                    gCGTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGc                      >  2:3903249/1‑93 (MQ=255)

                                                                      gTAGCGTGGAAGAGATCTTTCACGCGTTATGATTTACATCGGTCTGCCGTAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCATCAGCCt                    >  2:4212845/1‑93 (MQ=255)

                                                                      gTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCt                    >  2:1225855/1‑93 (MQ=255)

                                                                      gTAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCt                    <  2:2269601/93‑1 (MQ=255)

                                                                      gTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCGTAAATGGGTGCGGTTGGGGATCACCAGCCt                    >  1:674687/1‑93 (MQ=255)

                                                                       tAGCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCtt                   <  2:3057390/93‑1 (MQ=255)

                                                                         gCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGa                 <  2:2725300/93‑1 (MQ=255)

                                                                         gCGTGGAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGa                 <  1:4030704/93‑1 (MQ=255)

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

                                                                             ggAAGAGATCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAagag             >  2:2044106/1‑93 (MQ=255)

                                                                             ggAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAagag             >  2:4425892/1‑93 (MQ=255)

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

                                                                                 gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg         >  1:3602797/1‑93 (MQ=255)

                                                                                 gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg         >  2:2067572/1‑93 (MQ=255)

                                                                                 gagaTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATg         <  2:1696919/93‑1 (MQ=255)

                                                                                   gaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGcc       <  2:3463281/93‑1 (MQ=255)

                                                                                   gaTCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGcc       <  1:2312460/93‑1 (MQ=255)

                                                                                    aTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccc      >  2:2856586/1‑93 (MQ=255)

                                                                                    aTCCTTCACGCGTTATGATTTACATCGGTCTGCCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGTATCACCAGCCTTGAAGAGTATGccc      <  2:2794399/93‑1 (MQ=255)

TGTAGTGCCGCCAGCGCCGAGCAGCACAATAACAGTATCAACCATATCTACCCGCGCATCGCCCGTGTGCGTAGCGTGGAAGAGATCCTCAACGCGTTATGATTTACATCGGTCTACCGCAATGGTCGCATCCTAAATGGGTGCGGTTGGGGATCACCAGCCTTGAAGAGTATGCCCGCCA  >  NC_000913/1949321‑1949501

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