BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	260,537	G→C	15.1%	S50S (TCG→TCC)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	260,537	0	G	C	15.1%	366.9 / 47.3	132	S50S (TCG→TCC)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base G (53/58); new base C (5/15); total (58/74)
Fisher's exact test for biased strand distribution p-value = 8.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 135 total aligned reads displayed.

GGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACT  >  NC_000913/260439‑260635

ggCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCgg                                                                                                    <  1:2895538/100‑1 (MQ=255)

 gCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCggg                                                                                                   <  1:14210/100‑1 (MQ=255)

        gcgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATc                                                                                           >  2:2068448/1‑100 (MQ=255)

        gcgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATc                                                                                           >  2:279814/1‑100 (MQ=255)

         cgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCATTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcg                                                                                          <  2:1585151/100‑1 (MQ=255)

         cgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcg                                                                                          <  1:773056/100‑1 (MQ=255)

           ccGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGCTTGTTATTGTGACGTCGGGC‑GCCATcgcc                                                                                        >  2:3593274/1‑100 (MQ=255)

            cGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgccg                                                                                       >  1:1876468/1‑100 (MQ=255)

            cGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGGGACGTCGGGC‑GCGATcgccg                                                                                       >  2:2358144/1‑100 (MQ=255)

              tCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCgccgcc                                                                                     <  1:3419231/100‑1 (MQ=255)

              tCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCgccgcc                                                                                     <  1:418350/100‑1 (MQ=255)

               cTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGccgccg                                                                                    <  2:945044/100‑1 (MQ=255)

               cTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGccgccg                                                                                    >  1:424942/1‑100 (MQ=255)

               cTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGccgccg                                                                                    <  2:3005644/100‑1 (MQ=255)

               cTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCACCGGATAGTTATTGTGACGTCGGGC‑GCGATCGccgccg                                                                                    >  1:1981977/1‑100 (MQ=255)

                tGAACCGTGCCCATATCGTTGAACTTGTCCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCgg                                                                                   <  1:3448763/100‑1 (MQ=255)

                  aaCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGAc                                                                                 >  1:1792113/1‑100 (MQ=255)

                  aaCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccgac                                                                                  <  2:3624548/100‑8 (MQ=255)

                  aaCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccgac                                                                                  <  1:2888708/100‑8 (MQ=255)

                   aCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACg                                                                                >  2:1308758/1‑100 (MQ=255)

                   aCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACg                                                                                <  1:1344399/100‑1 (MQ=38)

                    ccGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACg                                                                                <  1:3434165/99‑1 (MQ=38)

                    ccGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGt                                                                               <  2:2149116/100‑1 (MQ=38)

                    ccGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGt                                                                               <  2:3313966/100‑1 (MQ=38)

                     cGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTg                                                                              <  1:920308/100‑1 (MQ=255)

                     cGTGCCCATATCGTTGAACTTGTTCGCCAGGGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTg                                                                              <  2:1986186/100‑1 (MQ=255)

                      gTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGa                                                                             <  2:2444536/100‑1 (MQ=255)

                      gTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGa                                                                             >  1:3105581/1‑100 (MQ=255)

                      gTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtca                                                                             <  1:2475068/100‑3 (MQ=37)

                       tGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAg                                                                            <  1:2931207/100‑1 (MQ=255)

                       tGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAg                                                                            <  1:1797728/100‑1 (MQ=255)

                       tGCCCATATCGTTGAACTTGTTCGCCAGGGCGCGAAGTTACAGGCCGGCGGGCATCGGATTGGTATTGTGACGTCCGGCGGCGACCGCGGCCggcgtcaa                                                                             <  2:829612/100‑7 (MQ=255)

                        gCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtcaca                                                                           <  2:3262729/100‑5 (MQ=37)

                         cccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCa                                                                          >  1:556825/1‑100 (MQ=255)

                          caatatCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtcacaat                                                                         <  2:3369264/98‑7 (MQ=25)

                           cATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCAcc                                                                        >  1:3392538/1‑100 (MQ=255)

                               tCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTggg                                                                    <  1:2867867/100‑1 (MQ=255)

                               tCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTggg                                                                    <  1:2552373/100‑1 (MQ=255)

                               tCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTggg                                                                    <  2:556823/100‑1 (MQ=255)

                                   tGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTAc                                                                >  1:1057942/1‑100 (MQ=255)

                                      aCTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCg                                                             <  1:1738074/100‑1 (MQ=255)

                                        ttGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGa                                                           <  1:1037655/100‑1 (MQ=255)

                                          gTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAAc                                                         >  2:252458/1‑100 (MQ=255)

                                          gTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAAc                                                         <  1:1666269/100‑1 (MQ=255)

                                          gTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAAc                                                         >  1:3047372/1‑100 (MQ=255)

                                           ttCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACt                                                        <  1:358381/100‑1 (MQ=255)

                                            tCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTg                                                       <  2:2369443/100‑1 (MQ=255)

                                              gCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGcc                                                     <  2:3047349/100‑1 (MQ=255)

                                              gCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGcc                                                     <  1:516185/100‑1 (MQ=255)

                                               ccAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGcc                                                     >  2:1771674/1‑99 (MQ=255)

                                                    gcgcgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGAc                                               <  2:2643816/100‑1 (MQ=255)

                                                     cgcgcAGTTACATCCCGCCGGGCATCGGATTGTTATTGTGACGGCGGGC‑GCGATCGCCGCCGGACGCGAGCACCTGCGTTCCCCGGCGCTCCCAGCGa                                                >  2:2800919/1‑98 (MQ=255)

                                                          aGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGc                                         <  1:2802617/100‑1 (MQ=255)

                                                           gTTACATGACGCCGGGCATCGGATTGTTATTGTGACGTCGTGC‑GCGATCGCCGCCGGACGTGAGCACCTATGTTACCCGGAACTGCCAGCGCCCATCGcc                                        <  2:1257584/100‑1 (MQ=255)

                                                            ttACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc                                        <  1:1238890/99‑1 (MQ=255)

                                                             tACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc                                      >  1:578951/1‑100 (MQ=255)

                                                             tACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc                                      <  2:2070744/100‑1 (MQ=255)

                                                             tACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc                                      >  1:1684933/1‑100 (MQ=255)

                                                             tACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc                                      >  1:1088965/1‑100 (MQ=255)

                                                             tACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGAGCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc                                      <  2:1008278/100‑1 (MQ=255)

                                                              tcttGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCg                                     <  1:2342321/97‑1 (MQ=255)

                                                              aCATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCg                                     >  1:1293086/1‑100 (MQ=255)

                                                                aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGaa                                   <  1:875909/100‑1 (MQ=255)

                                                                aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGaa                                   <  2:2356396/100‑1 (MQ=255)

                                                                aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGGTACCCGGAACTGCCAGCGACCATCGCCTCGaa                                   >  1:2681730/1‑100 (MQ=255)

                                                                 tGCCGCCGGGCATCGGACTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCTGAACGGCCAGCGACCATCGCCTCg                                     >  2:2173540/1‑97 (MQ=255)

                                                                  gccgccGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaac                                 <  1:3355089/100‑1 (MQ=255)

                                                                  gccgccGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaac                                 >  1:2807719/1‑100 (MQ=255)

                                                                   ccgccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaaca                                <  2:1909890/100‑1 (MQ=255)

                                                                   ccgccgGACATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaaca                                >  1:1909900/1‑100 (MQ=255)

                                                                     gccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaacaac                              <  2:688832/100‑1 (MQ=255)

                                                                     gccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaacaac                              <  1:1438268/100‑1 (MQ=255)

                                                                         ggCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgct                          <  1:1666247/100‑1 (MQ=255)

                                                                         ggCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgct                          >  1:542386/1‑100 (MQ=255)

                                                                          gCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgctg                         >  2:3101899/1‑100 (MQ=255)

                                                                           cATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTgg                        >  2:1700915/1‑100 (MQ=255)

                                                                            aacGGAATGGTATTGGGGCGTTGGGG‑GGGATCGGCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggc                       <  1:480324/98‑1 (MQ=255)

                                                                             tCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcg                      <  2:1859744/100‑1 (MQ=255)

                                                                             tCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggc                       >  2:3633108/1‑99 (MQ=255)

                                                                              cGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcgg                     >  1:463117/1‑100 (MQ=255)

                                                                               ggATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcggc                    >  1:354484/1‑100 (MQ=255)

                                                                               ggATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcggc                    >  1:354321/1‑100 (MQ=255)

                                                                               ggATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcggc                    <  2:1070623/100‑1 (MQ=255)

                                                                                  ttGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGcggcgg                  <  2:1301798/99‑1 (MQ=255)

                                                                                    gTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAg               >  2:465346/1‑100 (MQ=255)

                                                                                      tATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAggg             <  1:1530407/100‑1 (MQ=255)

                                                                                      ggtgctcaCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCGCCTGGGGTACCCGGAACTGCCAGCGACCATCGCCTCGAAACACCTGCTGGCGGCGATAg               >  2:1339664/8‑98 (MQ=16)

                                                                                        ttGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGc            <  2:141676/99‑1 (MQ=255)

                                                                                        ttGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCa           <  2:2055072/100‑1 (MQ=255)

                                                                                         tgtgACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag          <  1:390207/100‑1 (MQ=255)

                                                                                         gctcaCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag          >  2:2229268/5‑100 (MQ=37)

                                                                                         gctcaCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag          >  1:3171417/5‑100 (MQ=37)

                                                                                          gtgACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag          >  1:1565934/1‑99 (MQ=255)

                                                                                            gACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGACCTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCagag        >  1:391257/1‑99 (MQ=255)

                                                                                            caCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGt       >  1:203208/2‑100 (MQ=38)

                                                                                             aCGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTc      <  1:1536654/100‑1 (MQ=255)

                                                                                             aCGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTc      <  2:1104224/100‑1 (MQ=255)

                                                                                              cGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCg     >  1:1433639/1‑100 (MQ=255)

                                                                                              cGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCg     <  1:821978/100‑1 (MQ=255)

                                                                                                 cGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACt  >  2:1370115/1‑100 (MQ=255)

GGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACT  >  NC_000913/260439‑260635

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