BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	260,541:1	+G	15.3%	coding (154/1104 nt)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	260,540	1	.	G	15.3%	452.9 / 85.8	131	coding (153/1104 nt)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base . (53/58); new base G (5/15); total (58/73)
Fisher's exact test for biased strand distribution p-value = 8.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

Only 100 of 137 total aligned reads displayed.

GGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATT  >  NC_000913/260442‑260639

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

     gcgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATc                                                                                               >  2:1106626/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)

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

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

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

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

             tGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCgg                                                                                       <  1:1457964/100‑1 (MQ=255)

              gAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCgg                                                                                       >  2:3545014/1‑99 (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:2765162/100‑1 (MQ=38)

                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:3149868/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                                                                                 >  1:3105581/1‑100 (MQ=255)

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

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

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

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

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

                      cccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGGCGCCGGACGTGAGCa                                                                              >  1:3590056/1‑100 (MQ=255)

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

                       ccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCAc                                                                             >  2:2991233/1‑100 (MQ=255)

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

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

                        cATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCAcc                                                                            <  1:3495686/100‑1 (MQ=255)

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

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

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

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

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

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

                                       gTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGCTTGTTATTGTGACGTCGGGC‑GCGATCGCCGCAGGACGTGAGCAACTGGGGTACCCGGAAc                                                             >  2:2076002/1‑100 (MQ=255)

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

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

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

                                         tCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTg                                                           <  2:2369443/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)

                                                    cgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCAt                                                >  1:42111/1‑100 (MQ=255)

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

                                                        gTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc                                            <  1:2877192/100‑1 (MQ=255)

                                                         ttACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc                                            <  1:1238890/99‑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)

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

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

                                                              tGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGaga                                      >  2:2342305/1‑98 (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                                  <  1:1438268/100‑1 (MQ=255)

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

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

                                                                      ggCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgct                              <  1:1666247/100‑1 (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‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcg                          >  1:510948/1‑100 (MQ=255)

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

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

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

                                                                             naTTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGGCGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGgcggcg                       >  1:34293/2‑100 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                         gACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGt           >  1:3033015/1‑100 (MQ=255)

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

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

                                                                                          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)

                                                                                           cGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTCGCGGCGGTAGGGCAGAGTCg         >  2:1332192/1‑100 (MQ=255)

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

                                                                                                 gc‑gcGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGAt   <  2:313195/100‑1 (MQ=255)

                                                                                                  c‑gcgATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGAtt  <  2:2675085/100‑1 (MQ=255)

                                                                                                  c‑gcgATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGAtt  <  1:2027243/100‑1 (MQ=255)

GGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATT  >  NC_000913/260442‑260639

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