BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	339,557	T→C	G411G (GGT→GGC)	yahG →	DUF1116 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	339,557	0	T	C	96.5%	390.0 / 1.3	115	G411G (GGT→GGC)	yahG	DUF1116 family protein
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (57/54); total (59/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 117 total aligned reads displayed.

TGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGG  >  NC_000913/339460‑339652

tggtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGctt                                                                                               <  2:3216338/100‑1 (MQ=255)

 ggtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGcttt                                                                                              >  1:3224343/1‑100 (MQ=255)

    gggCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATg                                                                                           >  2:1078697/1‑100 (MQ=255)

    gggCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATg                                                                                           <  2:1577371/100‑1 (MQ=255)

     ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATgg                                                                                          >  2:1457897/1‑100 (MQ=255)

     ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATgg                                                                                          >  2:1457898/1‑100 (MQ=255)

     ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATgg                                                                                          >  1:2490989/1‑100 (MQ=255)

     ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATgg                                                                                          <  1:1943041/100‑1 (MQ=255)

      gcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATggg                                                                                         >  2:1243428/1‑100 (MQ=255)

        ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCg                                                                                       <  1:3029932/100‑1 (MQ=255)

        ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCg                                                                                       <  2:2535120/100‑1 (MQ=255)

          cACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTg                                                                                     <  2:2927221/100‑1 (MQ=255)

            cGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGcc                                                                                   >  1:573164/1‑100 (MQ=255)

             ggTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCATCGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCg                                                                                  <  2:1741931/100‑1 (MQ=255)

                ggAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCg                                                                               <  2:471814/100‑1 (MQ=255)

                  aagaagCTATTGCTTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGc                                                                             <  1:3186517/100‑1 (MQ=255)

                  aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGc                                                                             <  1:1921567/100‑1 (MQ=255)

                  aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGc                                                                             <  1:2149835/100‑1 (MQ=255)

                    gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCaa                                                                           >  1:1023255/1‑100 (MQ=255)

                      agCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATc                                                                         <  1:991500/100‑1 (MQ=255)

                      agCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATc                                                                         >  1:3052249/1‑100 (MQ=255)

                        cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGa                                                                       <  2:2098064/100‑1 (MQ=255)

                        cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGa                                                                       >  1:744574/1‑100 (MQ=255)

                        cTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGa                                                                       <  2:2098381/100‑1 (MQ=255)

                          aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCTAc                                                                      >  1:2768102/1‑99 (MQ=255)

                          aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACa                                                                     <  1:3006463/100‑1 (MQ=255)

                           ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACAt                                                                    >  1:363336/1‑100 (MQ=255)

                           ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACAt                                                                    <  2:1606641/100‑1 (MQ=255)

                            tGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                   >  2:510504/1‑100 (MQ=255)

                             gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATc                                                                   >  1:1933614/1‑99 (MQ=255)

                             gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCa                                                                  >  2:916060/1‑100 (MQ=255)

                             gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCa                                                                  <  1:2993914/100‑1 (MQ=255)

                             gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCa                                                                  <  1:2993940/100‑1 (MQ=255)

                               tttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                                                <  2:2455285/100‑1 (MQ=255)

                                 tctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                                              >  2:2660833/1‑100 (MQ=255)

                                 tctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                                              >  2:1188499/1‑100 (MQ=255)

                                  ctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCcgc                                                             >  1:3242962/1‑100 (MQ=255)

                                    cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGTCTTTATAGGCGTGCCGTCGGCAATCGACATCACCcgcg                                                            >  2:2292081/1‑99 (MQ=255)

                                    cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                           >  2:1499046/1‑100 (MQ=255)

                                    cccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGt                                                           >  1:2690422/1‑100 (MQ=255)

                                     ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTg                                                          >  2:2849592/1‑100 (MQ=255)

                                     ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTg                                                          <  1:1342691/100‑1 (MQ=255)

                                      cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTgg                                                         >  2:975355/1‑100 (MQ=255)

                                       gTCAGATGCGCGAAATCACCCTCGGTGATAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                        <  1:670045/100‑1 (MQ=255)

                                       gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                        <  1:1888292/100‑1 (MQ=255)

                                       gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTggg                                                        <  2:1942844/100‑1 (MQ=255)

                                        tCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGt                                                       <  1:1034202/100‑1 (MQ=255)

                                         cAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTa                                                      <  2:3025444/100‑1 (MQ=255)

                                           gATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagc                                                    >  1:57056/1‑100 (MQ=255)

                                             tGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                                                  >  1:3197763/1‑100 (MQ=255)

                                              gcgcGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                 <  1:178609/100‑1 (MQ=255)

                                              gcgcGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                                                 <  1:1304827/100‑1 (MQ=255)

                                               cgcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCg                                                <  1:902782/100‑1 (MQ=255)

                                               cgcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCg                                                >  1:3183340/1‑100 (MQ=255)

                                                gcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCgg                                               >  2:1331629/1‑100 (MQ=255)

                                                gcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCgg                                               >  2:2903583/1‑100 (MQ=255)

                                                gcgAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCgg                                               <  1:205456/100‑1 (MQ=255)

                                                  gAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                                             >  1:1565192/1‑100 (MQ=255)

                                                     aTCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTc                                          <  1:1170760/100‑1 (MQ=255)

                                                       cACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTg                                        >  2:2623296/1‑100 (MQ=255)

                                                        aCCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGc                                       <  2:2464623/100‑1 (MQ=255)

                                                        aCCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGc                                       <  1:1071205/100‑1 (MQ=255)

                                                         cccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                      <  2:228076/100‑1 (MQ=255)

                                                         cccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGcc                                      >  1:1043234/1‑100 (MQ=255)

                                                          ccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                     >  2:1726686/1‑100 (MQ=255)

                                                          ccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg                                     >  1:2431362/1‑100 (MQ=255)

                                                           cTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCgg                                    <  1:1457911/100‑1 (MQ=255)

                                                           cTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCgg                                    <  1:1457912/100‑1 (MQ=255)

                                                            tCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGt                                   <  2:1931528/100‑1 (MQ=255)

                                                             cGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTg                                  <  1:1109908/100‑1 (MQ=255)

                                                              ggTGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGa                                 >  2:3240160/1‑100 (MQ=255)

                                                                tGAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc                               <  2:1107862/100‑1 (MQ=255)

                                                                 gAAAACCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCa                              >  2:2890638/1‑100 (MQ=255)

                                                                  aaaaCCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa                             >  2:297223/1‑100 (MQ=255)

                                                                   aaaCCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAac                            >  2:1661184/1‑100 (MQ=255)

                                                                    aaCCCCAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca                           >  2:1483437/1‑100 (MQ=255)

                                                                       cccAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                        <  1:2849627/100‑1 (MQ=255)

                                                                        ccAACGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGc                       >  2:2683290/1‑100 (MQ=255)

                                                                        ccAACGTCACCATTCCGCTACTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCg                        >  2:365345/1‑99 (MQ=255)

                                                                          aaCGTCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCa                     <  2:2339446/100‑1 (MQ=255)

                                                                              tCACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGc                 <  2:594646/100‑1 (MQ=255)

                                                                               cACCCTTCCGCTGCTCGGCTTTATGGGCGTGCCGTCCGCAATCGACATCACCCGCGTGGGTAGCAGCGGCCTTCTGCCGGTGATCAACACCGCCATCGcc                <  1:3095430/100‑1 (MQ=255)

                                                                               cACCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGcc                <  1:1499061/100‑1 (MQ=255)

                                                                                aCCATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCg                  <  1:2660902/97‑1 (MQ=255)

                                                                                 ccATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGTAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa              >  1:2293480/1‑100 (MQ=255)

                                                                                 ccATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCa              >  1:2822061/1‑100 (MQ=255)

                                                                                  cATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCAt             >  1:911317/1‑100 (MQ=255)

                                                                                  cATTCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCAt             >  1:911332/1‑100 (MQ=255)

                                                                                    ttCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaa           >  2:2795000/1‑100 (MQ=255)

                                                                                     tCCGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATaaa          <  2:573155/100‑1 (MQ=255)

                                                                                      ccGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg         >  2:1949186/1‑100 (MQ=255)

                                                                                      ccGCTGCTCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAg         <  1:1848900/100‑1 (MQ=255)

                                                                                         ctgctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg      >  1:224411/1‑100 (MQ=255)

                                                                                         ctgctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATg      <  2:398545/100‑1 (MQ=255)

                                                                                           gctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCg    >  1:698093/1‑100 (MQ=255)

                                                                                           gctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCg    >  1:2653459/1‑100 (MQ=255)

                                                                                            ctCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCgg   >  2:2504731/1‑100 (MQ=255)

                                                                                             tCGGCTTTATGGGCGTGCGGTCGGCAATCGACATCACCCGCGTGGGTAGCGGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg  <  2:2249879/100‑1 (MQ=255)

                                                                                             tCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg  <  2:3224306/100‑1 (MQ=255)

                                                                                             tCGGCTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCggg  <  2:2003031/100‑1 (MQ=255)

TGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACACCGCCATCGCCCATAAAGATGCGGG  >  NC_000913/339460‑339652

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