BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,134,509	A→G	E274E (GAA→GAG)	b3943 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,134,509	0	A	G	99.1%	1005.4 / ‑1.8	321	E274E (GAA→GAG)	b3943	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (138/180); total (138/183)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 334 total aligned reads displayed.

CAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTT  >  NC_000913/4134419‑4134598

  gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTa                                                                                       <  1:7128906/93‑1 (MQ=255)

  gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTa                                                                                       >  1:3231323/1‑93 (MQ=255)

  gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTa                                                                                       <  2:6891625/93‑1 (MQ=255)

  gTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTa                                                                                       >  1:822425/1‑93 (MQ=255)

    aTTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGc                                                                                     <  1:7410685/93‑1 (MQ=255)

        ccAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTa                                                                                 >  1:7282872/1‑93 (MQ=255)

        ccAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTa                                                                                 <  1:6234469/93‑1 (MQ=255)

        ccAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGCAGCggga                                                                                 >  2:2237051/1‑91 (MQ=255)

          aGCTTAAGCACACTGGCCGTTCCCGTTTGCGTCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGCGGGCCGGTGCGGGTAGCGGCatt                                                                               >  1:5963819/1‑93 (MQ=255)

           gCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattg                                                                              >  1:1410341/1‑93 (MQ=255)

            cTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgt                                                                             <  2:2843184/93‑1 (MQ=255)

             ttAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtg                                                                            >  2:14059/1‑93 (MQ=255)

             ttAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtg                                                                            >  1:6853248/1‑93 (MQ=255)

              tAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgc                                                                           <  2:6460531/93‑1 (MQ=255)

              tAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgc                                                                           <  1:2885138/93‑1 (MQ=255)

                 gCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctga                                                                        <  2:6289958/93‑1 (MQ=255)

                 gCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctga                                                                        <  1:6178953/93‑1 (MQ=255)

                 gCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctga                                                                        >  2:4664358/1‑93 (MQ=255)

                  cacaCTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgat                                                                       <  2:6837531/93‑1 (MQ=255)

                    cacTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattg                                                                     >  1:1771816/1‑93 (MQ=255)

                    cacTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattg                                                                     <  2:2708183/93‑1 (MQ=255)

                    cacTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattg                                                                     >  1:4854322/1‑93 (MQ=255)

                     ccTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgt                                                                    <  1:5556269/92‑1 (MQ=255)

                       tGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgc                                                                  <  2:5673768/93‑1 (MQ=255)

                        ggCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgct                                                                 <  1:1190235/93‑1 (MQ=255)

                         gCCGTTCCCGTTTGCGGAGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTattgtgctgattgtgctg                                                                <  2:2512524/93‑1 (MQ=255)

                           ctttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGc                                                              <  1:1621187/91‑1 (MQ=255)

                           cGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGc                                                              >  1:2596886/1‑93 (MQ=255)

                           cGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGc                                                              <  1:5755107/93‑1 (MQ=255)

                           cGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATACGGGGGCCGTTGAGGTCAGCGGGATTGTGCTCATTGTGCTGGc                                                              >  2:2191958/1‑93 (MQ=255)

                            gTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCa                                                             >  2:4746462/1‑93 (MQ=255)

                            gTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCa                                                             >  2:845957/1‑93 (MQ=255)

                            gTTCCCGTGTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCa                                                             <  1:1043019/93‑1 (MQ=255)

                             ttCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCAc                                                            >  2:5109743/1‑93 (MQ=255)

                              tCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCAcg                                                           >  2:5528856/1‑92 (MQ=255)

                              tCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACt                                                           <  2:5985634/93‑1 (MQ=255)

                              tCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACt                                                           <  1:4043622/93‑1 (MQ=255)

                                 cGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGgc                                                        >  1:5861426/1‑93 (MQ=255)

                                  gTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGgcg                                                       >  2:3499452/1‑93 (MQ=255)

                                   tttGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGgcgc                                                      >  2:1900087/1‑93 (MQ=255)

                                   tttGCGGCGTACTCTTTTCCTGGGGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGAGCTGGCACTGgcgt                                                      >  2:5440101/1‑92 (MQ=255)

                                     tGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTg                                                    >  2:3202993/1‑93 (MQ=255)

                                       cggcgTACTCTTTTCCTGGAGGCTGCTCGGCGAGAATCCGGGGGCCGATGAGGGTAGCGGTATTGTGCTGATTGAGCTGGCACTGGCGCTGGt                                                  >  2:3553687/1‑93 (MQ=255)

                                        ggcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTg                                                 <  2:3334697/93‑1 (MQ=255)

                                         gcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGa                                                <  2:5706530/93‑1 (MQ=255)

                                         gcgTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGa                                                <  2:2491220/93‑1 (MQ=255)

                                         gcgTACTCTTTTCCTGATGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGa                                                <  2:443616/93‑1 (MQ=255)

                                          ggttctctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAg                                               <  1:1990939/89‑1 (MQ=255)

                                           gTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGc                                              <  1:4569371/93‑1 (MQ=255)

                                           gTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGc                                              >  2:1667346/1‑93 (MQ=255)

                                            tACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGcc                                             >  2:2347943/1‑93 (MQ=255)

                                            tACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGcc                                             >  2:7769209/1‑93 (MQ=255)

                                             aCTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCg                                            >  1:110607/1‑93 (MQ=255)

                                               tctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTa                                          >  1:2773208/1‑93 (MQ=255)

                                               tctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTa                                          >  1:2704301/1‑93 (MQ=255)

                                                ctTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTaa                                         <  2:7312589/93‑1 (MQ=255)

                                                ctTTGCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTaa                                         <  2:2116693/93‑1 (MQ=255)

                                                 ttttCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAg                                        <  2:7158547/93‑1 (MQ=255)

                                                 ttttCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAg                                        >  2:3835107/1‑93 (MQ=255)

                                                  tttCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCTTAAGa                                       <  1:550214/93‑1 (MQ=255)

                                                     ccTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaa                                    <  1:3818599/93‑1 (MQ=255)

                                                     ccTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaa                                    <  1:4709101/93‑1 (MQ=255)

                                                      cTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaaa                                   >  1:4948401/1‑93 (MQ=255)

                                                       tggtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGaaaaaa                                  <  2:1021316/93‑1 (MQ=255)

                                                        ggtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaag                                 <  2:2654633/93‑1 (MQ=255)

                                                         gtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaaga                                >  1:2121359/1‑93 (MQ=255)

                                                         gtggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaaga                                >  1:4375202/1‑93 (MQ=255)

                                                           ggCTGCTCGGCGAGAATCCGGGGGCTGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCCCTGGCGCTGGTGAGCCGTAAGAAAAaagaag                              <  2:1951425/93‑1 (MQ=255)

                                                           ggCTGCTCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAaagaag                              <  2:4044636/93‑1 (MQ=255)

                                                            gctgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGc                             >  2:3047201/1‑93 (MQ=255)

                                                             ctgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGcc                            <  2:3661135/93‑1 (MQ=255)

                                                              tgctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCg                           >  1:662793/1‑93 (MQ=255)

                                                               gctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGt                          <  1:3486890/93‑1 (MQ=255)

                                                                ctCGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTc                         <  2:4943961/93‑1 (MQ=255)

                                                                  cGGCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAg                       <  1:523404/93‑1 (MQ=255)

                                                                   ggCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAgt                      <  1:2499776/93‑1 (MQ=255)

                                                                    gCGAGAATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAgtg                     >  2:895082/1‑93 (MQ=255)

                                                                      gagaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTa                   <  2:7846723/93‑1 (MQ=255)

                                                                      gagaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTa                   >  1:6933862/1‑93 (MQ=255)

                                                                       agaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaa                  >  1:3913836/1‑93 (MQ=255)

                                                                       agaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaa                  <  2:7614860/93‑1 (MQ=255)

                                                                        gaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaaa                 >  1:354846/1‑93 (MQ=255)

                                                                        gaATCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaaa                 <  1:3025021/93‑1 (MQ=255)

                                                                          aTCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaaaaa               <  1:2264888/93‑1 (MQ=255)

                                                                          aTCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTaaaaa               <  2:6711498/93‑1 (MQ=255)

                                                                           tCCGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAg              <  1:7788976/93‑1 (MQ=255)

                                                                             cGGGGGCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGa            <  2:4610989/93‑1 (MQ=255)

                                                                               ggggCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATc          <  2:3901291/93‑1 (MQ=255)

                                                                               ggggCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATc          <  1:3158986/93‑1 (MQ=255)

                                                                                 ggCCGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTg        <  2:3140083/93‑1 (MQ=255)

                                                                                   ccGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGaa      <  1:341846/93‑1 (MQ=255)

                                                                                   ccGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGaa      <  1:6339342/93‑1 (MQ=255)

                                                                                   ccGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGaa      <  2:4889001/93‑1 (MQ=255)

                                                                                    cGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAt     <  2:707234/93‑1 (MQ=255)

                                                                                    cGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAt     <  2:4455760/93‑1 (MQ=255)

                                                                                    cGTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAt     <  2:5254788/93‑1 (MQ=255)

                                                                                     gTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAtt    <  2:6130876/93‑1 (MQ=255)

                                                                                     gTTGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAtt    >  1:4437807/1‑93 (MQ=255)

                                                                                      ttGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAATAAGAAGCCGTCAGTGTAAAAAGGATCTGAAttt   >  1:7653920/1‑93 (MQ=255)

                                                                                      ttGAGGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAAttt   >  1:2063253/1‑93 (MQ=255)

CAGTATTGCCAGCTTAAGCACACTGGCCGTTCCCGTTTGCGGCGTACTCTTTTCCTGGTGGCTGCTCGGCGAGAATCCGGGGGCCGTTGAAGGTAGCGGTATTGTGCTGATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTT  >  NC_000913/4134419‑4134598

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