BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,123,585	+G	coding (953/2286 nt)	glcF ←	glycolate oxidase iron‑sulfur subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,123,582	1	.	G	99.0%	957.9 / 1.2	287	coding (956/2286 nt)	glcF	glycolate oxidase iron‑sulfur subunit
Reads supporting (aligned to +/- strand): ref base . (1/2); new base G (147/136); total (149/138)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 324 total aligned reads displayed.

AGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑GG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCGGCTG  >  NC_000913/3123490‑3123669

aGTTCCGCGTACAGGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGc                                                                                            <  1:7985600/93‑1 (MQ=255)

  ttCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑gg                                                                                         <  1:6643715/93‑1 (MQ=255)

   tCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑gg‑g                                                                                       >  1:1000267/1‑93 (MQ=255)

      gcgTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGACTAAAGGAGCCGATAGC‑gg‑ggca                                                                                    >  1:81632/1‑90 (MQ=255)

      gcgTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑gg‑ggca                                                                                    >  2:5097518/1‑90 (MQ=255)

      gcgTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑gg‑ggca                                                                                    <  1:1322884/93‑4 (MQ=255)

         tACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAg                                                                                 >  2:6913946/1‑93 (MQ=255)

          aCATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGc                                                                                >  1:5222139/1‑93 (MQ=255)

           cATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGcc                                                                               >  1:7082246/1‑93 (MQ=255)

           cATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGcc                                                                               >  2:8484607/1‑93 (MQ=255)

                gACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCAt                                                                          >  2:2964739/1‑93 (MQ=255)

                 aCCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCAtc                                                                         <  2:4565627/93‑1 (MQ=255)

                 aCCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑GGCG‑CAAAGACACCAtc                                                                         >  1:4623199/1‑93 (MQ=255)

                  ccGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCAtct                                                                        >  2:7391669/1‑93 (MQ=255)

                  ccGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCAtct                                                                        >  2:1714230/1‑93 (MQ=255)

                   cGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCAtctc                                                                       >  1:3724427/1‑93 (MQ=255)

                    ggggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTcc                                                                      <  2:1416279/93‑1 (MQ=255)

                     gggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCg                                                                     >  1:5428283/1‑93 (MQ=255)

                     gggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCg                                                                     <  1:1701965/93‑1 (MQ=255)

                      ggTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCCGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCgg                                                                    <  2:7018057/93‑1 (MQ=255)

                       gTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGc                                                                   <  2:3921109/93‑1 (MQ=255)

                          aaaCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACt                                                                >  1:2525038/1‑93 (MQ=255)

                          aaaCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACt                                                                <  2:1606633/93‑1 (MQ=255)

                           aaCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTa                                                               <  2:6414196/93‑1 (MQ=255)

                           aaCACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTa                                                               <  1:2067225/93‑1 (MQ=255)

                             cacGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTaaa                                                             <  1:6479284/93‑1 (MQ=255)

                              acGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAg                                                            >  1:7395939/1‑93 (MQ=255)

                               cgccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGc                                                           <  1:7126347/93‑1 (MQ=255)

                                 ccgcAAGGGTCGAGCTGCTTTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCgg                                                         <  1:6944700/93‑1 (MQ=255)

                                 ccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCgg                                                         <  1:3244760/93‑1 (MQ=255)

                                 ccgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCgg                                                         >  2:2561209/1‑93 (MQ=255)

                                  cgcAAGGGTGGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCggg                                                        <  1:5286242/93‑1 (MQ=255)

                                  cgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCggg                                                        >  1:4216767/1‑93 (MQ=255)

                                  cgcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCggg                                                        <  1:924991/93‑1 (MQ=255)

                                   gcAAGGTGCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGt                                                       <  1:4697739/93‑1 (MQ=255)

                                   gcAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGt                                                       >  2:1069448/1‑93 (MQ=255)

                                    cAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTc                                                      >  1:42924/1‑93 (MQ=255)

                                    cAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTc                                                      <  1:8300110/93‑1 (MQ=255)

                                    cAAGGGTCGAGCTGCGGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑GG‑GACAAAGCCACCATCTCCGGCACTAAAGCGGGTc                                                      <  2:8429447/93‑1 (MQ=255)

                                     aaGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCg                                                     >  1:8053381/1‑93 (MQ=255)

                                     aaGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCg                                                     >  1:3935895/1‑93 (MQ=255)

                                      aGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGc                                                    >  2:7749818/1‑93 (MQ=255)

                                      aGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGc                                                    >  2:6821564/1‑93 (MQ=255)

                                      aGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGc                                                    >  2:1546322/1‑93 (MQ=255)

                                       gggTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGGCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCa                                                   <  1:4503542/93‑1 (MQ=255)

                                        ggTCGGGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                                  <  2:210514/93‑1 (MQ=255)

                                        ggTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCAt                                                  >  2:6583192/1‑93 (MQ=255)

                                         gTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATg                                                 >  2:6274111/1‑93 (MQ=255)

                                          gcGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGa                                                <  2:2146179/92‑1 (MQ=255)

                                           cGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGAc                                               >  1:1486342/1‑93 (MQ=255)

                                            gAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGAcc                                              >  2:7933980/1‑93 (MQ=255)

                                             aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                                             >  1:5004150/1‑93 (MQ=255)

                                             aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                                             >  1:8591437/1‑93 (MQ=255)

                                             aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                                             >  2:3087783/1‑93 (MQ=255)

                                             aGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCg                                             >  2:7709680/1‑93 (MQ=255)

                                              gctgctGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGc                                            >  1:7485728/1‑93 (MQ=255)

                                               ctgctgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGcc                                           >  2:8230240/1‑93 (MQ=255)

                                               ctgctgTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGcc                                           >  2:3699245/1‑93 (MQ=255)

                                                    gTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGt                                      >  1:5508877/1‑93 (MQ=255)

                                                    gTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGt                                      >  1:2228700/1‑93 (MQ=255)

                                                     tttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGtt                                     <  2:2722301/93‑1 (MQ=255)

                                                     tttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGtt                                     <  1:8226146/93‑1 (MQ=255)

                                                     tttAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGtt                                     <  1:3095750/93‑1 (MQ=255)

                                                       tAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCTGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGc                                   >  1:4410763/1‑93 (MQ=255)

                                                       tAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGc                                   <  2:3568171/93‑1 (MQ=255)

                                                        aaGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCg                                  >  1:567342/1‑93 (MQ=255)

                                                         aGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCgg                                 >  2:3396779/1‑93 (MQ=255)

                                                          gctgctGGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGCCCGCCAGCGTTGCggg                                <  1:3821466/93‑1 (MQ=255)

                                                           ctgctgGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgc                               <  2:5946630/93‑1 (MQ=255)

                                                            tgctgGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcg                              >  2:1493631/1‑93 (MQ=255)

                                                             gctgGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcga                             <  2:4170587/93‑1 (MQ=255)

                                                               tgGTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatg                           >  2:6651793/1‑93 (MQ=255)

                                                               tgGTGATAGCGGAAAAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatg                           <  2:1686650/93‑1 (MQ=255)

                                                                ggTGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgc                          <  1:5583356/93‑1 (MQ=255)

                                                                  tGATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcga                        >  1:593880/1‑93 (MQ=255)

                                                                   gATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcgat                       <  1:1244586/93‑1 (MQ=255)

                                                                   gATAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGgcgatgcgat                       >  2:6476065/1‑93 (MQ=255)

                                                                    aTAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGcgatgcgatg                      >  1:2374575/1‑93 (MQ=255)

                                                                    aTAGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGcgatgcgatg                      >  1:8574399/1‑93 (MQ=255)

                                                                    aTAGCGGAATAAAGGAGCCGATAGC‑GG‑GACAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGcgatgcgatg                      <  2:4922912/93‑1 (MQ=255)

                                                                     aggCGGAAAAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACGATCTCCGGCACTAATGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATgg                     <  1:6807749/91‑1 (MQ=255)

                                                                      aGCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGa                    <  2:1303812/93‑1 (MQ=255)

                                                                       gCGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAt                   >  1:2763556/1‑93 (MQ=255)

                                                                       gCGGAAAAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAt                   <  2:6220820/93‑1 (MQ=255)

                                                                        cGGAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAtt                  <  1:2843406/93‑1 (MQ=255)

                                                                         ggAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGAttt                 >  1:5074738/1‑93 (MQ=255)

                                                                          gAATAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCGGCGTTGCGGGCGATGCGATGGATTTg                >  1:1526788/1‑93 (MQ=255)

                                                                           aaTAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGa               <  2:4691056/93‑1 (MQ=255)

                                                                             tAAAGGAGGCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGAGGGATTTGAtt             >  2:5482455/1‑93 (MQ=255)

                                                                             tAAAGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGAtt             >  1:3445961/1‑93 (MQ=255)

                                                                             aaaaGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGAtt             <  1:1613679/92‑1 (MQ=255)

                                                                                aGGAGCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGAGATGGATTTGATTGTc          >  1:6118139/1‑93 (MQ=255)

                                                                                    gCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCg      <  1:775745/93‑1 (MQ=255)

                                                                                    gCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCg      >  2:5680688/1‑93 (MQ=255)

                                                                                    gCCGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCg      <  1:1433351/93‑1 (MQ=255)

                                                                                     ccGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCgg     >  2:86268/1‑93 (MQ=255)

                                                                                      cGATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCGGc    >  1:1868426/1‑93 (MQ=255)

                                                                                      cGATAGCGGG‑G‑CAAAGCCACCAGCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCGGc    <  1:6841571/93‑1 (MQ=255)

                                                                                       gATAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCGGCt   >  1:4831957/1‑93 (MQ=255)

                                                                                        aTAGCGGG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATATGATTGTCCTCGGCTg  <  2:7836380/93‑1 (MQ=255)

AGTTCCGCGTACATGCGACCGGGGTTAAACACGCCGCAAGGGTCGAGCTGCTGTTTAAGCTGCTGGTGATAGCGGAATAAAGGAGCCGATAGC‑GG‑G‑CAAAGCCACCATCTCCGGCACTAAAGCGGGTCGCATGACCGCCAGCGTTGCGGGCGATGCGATGGATTTGATTGTCCTCGGCTG  >  NC_000913/3123490‑3123669

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