BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	260,537	G→C	15.2%	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.2%	284.5 / 39.8	108	S50S (TCG→TCC)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base G (46/42); new base C (6/10); total (53/55)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 110 total aligned reads displayed.

CGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTG  >  NC_000913/260441‑260636

cGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCgggn                                                                                                   <  2:2019201/100‑2 (MQ=255)

 ggATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGgc‑g                                                                                                 >  2:683803/1‑100 (MQ=255)

 ggATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGgc‑g                                                                                                 >  2:559882/1‑100 (MQ=255)

      gcgcCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATc                                                                                            <  2:1245567/100‑1 (MQ=255)

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

       cgcCGTCTGAATCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcg                                                                                           >  2:1496175/1‑100 (MQ=255)

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

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

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

         ccGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgcc                                                                                         <  2:1226682/100‑1 (MQ=255)

         ccGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgcc                                                                                         <  2:265099/100‑1 (MQ=255)

           gTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgccgc                                                                                       <  2:1768786/100‑1 (MQ=255)

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

              tGAACCGTGCCCATATCGTTGAAATTGATCGCCAGTGCGCACCGTTACAAGCCGCCGGACCACGGATTGTTATTGTACCGTCGGGA‑GCGATCACCGCCgg                                                                                    >  2:1509488/1‑100 (MQ=255)

               gAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGa                                                                                   >  1:35458/1‑100 (MQ=255)

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

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

                 aCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACg                                                                                 <  1:2505443/100‑1 (MQ=255)

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

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

                 aCCGTGCCAATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACg                                                                                 <  1:1087770/100‑1 (MQ=38)

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

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

                   cGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtc                                                                               <  2:2315138/100‑2 (MQ=38)

                   cGTGCCCAGATCGTTGAACGTTTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtc                                                                               <  2:1771814/100‑2 (MQ=25)

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

                     tGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCAGGC‑GCGATCGCCGCCGGAc                                                                                  <  1:1033277/95‑1 (MQ=255)

                       cccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCa                                                                           <  2:1876492/100‑1 (MQ=255)

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

                        ccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtcacaat                                                                          <  2:2590894/100‑7 (MQ=25)

                        ccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtcacaat                                                                          <  1:2557639/100‑7 (MQ=25)

                         cctatCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGCCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtcacaaga                                                                         <  1:1772509/98‑8 (MQ=25)

                          atatCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCt                                                                        >  2:2211455/1‑100 (MQ=255)

                           tatCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTg                                                                       >  2:2339627/1‑100 (MQ=255)

                            atCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                                                      <  2:309914/100‑1 (MQ=255)

                            atCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTgg                                                                      >  2:1038621/1‑100 (MQ=255)

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

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

                              cGTTGAACTTGTTGGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGt                                                                    >  2:1758012/1‑100 (MQ=255)

                                   aaCTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTAccc                                                               >  2:1090609/1‑100 (MQ=255)

                                    aCTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCg                                                              >  2:402015/1‑100 (MQ=255)

                                     cTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCg                                                              <  1:2124492/99‑1 (MQ=255)

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

                                       tGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGaa                                                           >  2:1200391/1‑100 (MQ=255)

                                       tGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGaa                                                           >  2:888935/1‑100 (MQ=255)

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

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

                                         ttCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACt                                                         >  1:106951/1‑100 (MQ=255)

                                           cGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGc                                                       >  2:791484/1‑100 (MQ=255)

                                           cGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGc                                                       >  2:2323962/1‑100 (MQ=255)

                                           cGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGc                                                       >  1:1874504/1‑100 (MQ=255)

                                           cGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGc                                                       >  2:512810/1‑100 (MQ=255)

                                               aGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGc                                                   <  1:2188745/100‑1 (MQ=255)

                                                gTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCg                                                  <  2:2169223/100‑1 (MQ=255)

                                                  gcgcgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGAc                                                >  2:1592679/1‑100 (MQ=255)

                                                   cgcgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGAcc                                               >  2:900559/1‑100 (MQ=255)

                                                   cgcgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGAcc                                               >  1:2544775/1‑100 (MQ=255)

                                                      gcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATc                                            <  1:1090605/100‑1 (MQ=255)

                                                       cAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGAc                                                >  1:2373943/1‑95 (MQ=255)

                                                       cAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCg                                           >  1:1320495/1‑100 (MQ=255)

                                                         gTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc                                         >  2:770555/1‑100 (MQ=255)

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

                                                            aCATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCg                                      <  1:512001/100‑1 (MQ=255)

                                                               tGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGaaa                                   >  1:1727992/1‑100 (MQ=255)

                                                               tGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGaaa                                   <  1:523022/100‑1 (MQ=255)

                                                                 ccgccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaaca                                 >  2:1540712/1‑100 (MQ=255)

                                                                   gccgGGCATCGGATTGTTATTGTGCCGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaacaac                               <  1:2152809/100‑1 (MQ=255)

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

                                                                   gccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCACCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaacaac                               >  1:2053264/1‑100 (MQ=255)

                                                                      gggCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgc                            <  2:1654678/100‑1 (MQ=255)

                                                                        gCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActgctg                          <  2:2154314/100‑1 (MQ=255)

                                                                          aTCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggc                        <  2:949941/100‑1 (MQ=255)

                                                                          aTCGGATTGTTATTGTGACGTCAGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggc                        <  1:864763/100‑1 (MQ=255)

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

                                                                               aTTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGcggcgg                   <  1:2574593/100‑1 (MQ=255)

                                                                               aTTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCAGCGCCTCGAAACAACTGCTGGcggcgg                   <  2:1580547/100‑1 (MQ=255)

                                                                                  gTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAg                <  2:1147152/100‑1 (MQ=255)

                                                                                  gTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGATACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAg                >  1:1215595/1‑100 (MQ=255)

                                                                                    tATTGTGACGTCGGGC‑GCGCTCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAggg              >  1:1316326/1‑100 (MQ=255)

                                                                                     gtgctcaCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGc             >  1:2379909/7‑100 (MQ=25)

                                                                                     gtgctcaCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGc             >  2:774961/7‑100 (MQ=25)

                                                                                     aTTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGc             <  2:1252293/100‑1 (MQ=255)

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

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

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

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

                                                                                        gtgACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCaga          >  1:1624643/1‑100 (MQ=255)

                                                                                        gtgACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCaga          >  2:2378539/1‑100 (MQ=255)

                                                                                         tcaCGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAggg              >  1:2573143/3‑95 (MQ=37)

                                                                                          gACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGt        <  1:2229844/100‑1 (MQ=255)

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

                                                                                            cGTCTGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTc       <  1:2051484/99‑1 (MQ=255)

                                                                                            cGTCGGGC‑GCGATCGCCGCCGGGCGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCg      >  2:2503712/1‑100 (MQ=255)

                                                                                            cGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGt        >  1:2386833/1‑98 (MQ=255)

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

                                                                                            cGTCCGGCGGCGATCG‑CGCCCGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAgggagatcgga      >  1:2013996/1‑92 (MQ=37)

                                                                                              tCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGAc    >  2:653643/1‑100 (MQ=255)

                                                                                               cGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACt   <  1:1592670/100‑1 (MQ=255)

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

                                                                                                gggC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTg  >  1:1893669/1‑100 (MQ=255)

CGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTG  >  NC_000913/260441‑260636

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